CONTENTS 


ONE INTRODUCTION I 

A Brief History 1 

A Brief Introduction to the Technology of Aeronautics 9 

TWO FLUID MECHANICS 16 

Fluid Statics and the Atmosphere 16 

Fluid Dynamics 22 

Conservation of Mass 26 

The Momentum Theorem 29 

Euler’s Equation of Motion 31 

Bernoulli’s Equation 33 

Determination of Free-Stream Velocity 35 

Determination of True Airspeed 36 

Potential Flow 38 

Velocity Potential and Stream Function 39 

Elementary Flow Functions 42 

Vortex 43 

Source 44 

Biot-Savart Law 45 

The Calculation of Flows for Well-Defined Body Shapes 47 

The Circular Cylinder 51 


The Numerical Calculation of Potential Flow Around Arbitrary Body 
Shapes 


THREE THE GENERA TION OF LIFT 
Wing Geometry 
Airfoils 

Airfoil Families 

NACA Four-Digit Series 
NACA Five-Digit Series 
NACA 1-Series (Series 16) 

NACA 6-Series 
Modern Airfoil Developments 
Prediction of Airfoil Behavior 
Maximum Lift 
Flaps 

Plain Flaps 
Split Flaps 
Slotted Flaps 


61 

61 

63 

72 

72 

72 

74 

74 

76 

82 

93 

95 


8SS 



CONTENTS 


ONE INTRODUCTION 1 

A Brief History 1 

A Brief Introduction to the Technology of Aeronautics 9 

TWO FLUID MECHANICS 16 

Fluid Statics and the Atmosphere 16 

Fluid Dynamics 22 

Conservation of Mass 26 

The Momentum Theorem 29 

Euler’s Equation of Motion 31 

Bernoulli’s Equation 33 

Determination of Free-Stream Velocity 35 

Determination of True Airspeed 36 

Potential Flow 38 

Velocity Potential and Stream Function 39 

Elementary Row Functions 42 

Vortex 43 

Source 44 

Biot-Savart Law 45 

The Calculation of Rows for Well-Defined Body Shapes 47 

The Circular Cylinder 51 

The Numerical Calculation of Potential Flow Around Arbitrary Body 

Shapes ^4 

THREE THE GENERATION OF LIFT 61 

Wing Geometry 

Airfoils **3 

Airfoil Families 72 

NACA Four-Digit Series 72 

NACA Five-Digit Series 72 

NACA 1-Series (Series 16) 74 

NACA 6-Series 74 

Modern Airfoil Developments 7< * 

Prediction of Airfoil Behavior * 2 

Maximum Lift 93 

Raps 95 

Plain Raps 99 

Split Raps 99 

Slotted Raps ^ 


Xii CONTENTS 


Flap Effectiveness in the Linear Range 
Leading Edge Devices 
The Optimum Airfoil for High Lift 
Powered-Lift Systems 
The Lifting Characteristics of a Finite Wing 
The Vortex System for a Wing 
The Maximum Lift of a Finite Wing 
Effect of Fuselage on C^ 

Effect of Trim on C f , ia> 

Estimation of C /nax for a Complete Airplane Configuration 
Airfoil Characteristics at Low Reynolds Numbers 

FOUR DRAG 

Skin Friction Drag 
Form Drag 

Drag of Streamlined Shapes 
Interference Drag 
Induced Drag 

Calculation of Induced Drag 
Effective Aspect Ratio 

Drag Breakdown and Equivalent Flat-Plate Area 
Drag Counts 

Average Skin Friction Coefficients 
Example Estimates of Drag Breakdown 
Trim Drag 
Cooling Drag 
Drag Reduction 
Winglets 

Reduction of Skin Friction Drag 
Drag Cleanup 
Total Airplane Drag 

FIVE LIFT AND DRAG AT HIGH MACH NUMBERS 

Qualitative Behavior of Airfoils as a Function of Mach Number 
Subsonic Flow at High Mach Numbers 
Fundamentals of Gas Dynamics 
One-Dimensional Isentropic Flow 
Normal Shock Waves 
Oblique Shock Waves 
Expansion Waves 
Transonic Airfoils 
Supersonic Airfoils 

Linearized Compressible Potential Flow 
Subsonic Flow 

Supersonic Flow (Ackeret Theory) 

Three-Dimensional Wings 
Characteristics of Sweptback Wings 
Delta Wings 


CONTENTS 


xiii 


Supersonic Wings 303 

Subsonic Leading Edges 304 

Supersonic Leading Edges 309 

Effect of Mach Number on the Zero Lift Drag of Two- and Three- 

Dimensional Shapes 315 

Area Rule for Transonic Flow 321 


SIX THE PRODUCTION OF THRUST 
A Brief History of the Piston Engine 
Piston Engine Characteristics 
Supercharged Engines 
Propeller Analysis 
Momentum Theory 
Blade Element Theories 
Momentum-Blade Element Theory 
Vortex Theory 

Practical Use of Propeller Charts 
Approximate Useful Relationships for Propellers 
Propeller Selection 
Design of a New Propeller 
A Brief History of the Turbojet 
Description of the Gas Turbine Engine 
Engine Ratings 
Flat Rating 

Some Considerations Relating to Gas Turbine Performance 
Qualitative Comparison of the Performance of Turbojet, Turbofan, and 
Turboprop Engines 

Specific Engine Characteristics and Performance 
Turbojet 
Turbofan 
Turboprop 
Installation Losses 
Trends in Aircraft Propulsion 


332 

332 

333 
339 
343 
343 
347 
349 
351 
359 
364 
366 
368 

370 

371 

377 

378 


382 

385 

385 

398 

405 

411 

411 


SEVEN AIRPLANE PERFORMANCE 

Takeoff 
Ground Roll 
Effect of Wind 
Airborne Distance 
Balanced Field Length 
Rate of Climb, Time to Climb, and Ceilings 
Generalized Power-Required Curve 
Time to Climb 
Range 

Maximum Endurance 

Descent 

Landing 


418 

418 

419 
425 
427 

431 

432 

438 

439 

440 

446 

447 

448 



XiV CONTENTS 


Airborne Distance 
Ground Roll 
Range Payload 
Operating Limitations 
Flight Envelope 

Maneuvering Envelope ( V-n Diagram) 
Gust Load Factors 

Energy Methods for Optimal Trajectories 
The Art of Estimating and Scaling 


450 

451 

452 
460 
460 
463 
466 
468 
471 


EIGHT STATIC STABILITY AND CONTROL 477 

Introduction 477 

Coordinate System—Forces, Moments, and Velocities 478 

Longitudinal Static Stability 479 

Stick-Fixed Stability 479 

Stick-Fixed Neutral Point and Static Margin 481 

Cwi, and Aerodynamic Center Location for a Finite Wing 484 

Downwash Angle 486 

Longitudinal Control 488 

Control Position as a Function of Lift Coefficient 488 

All-Movable Tail 489 

Stabilizer-Elevator 489 

Stabilator 491 

Control Forces 491 

Gearing 492 

Stick Force for A Stabilator 493 

Stick Force for a Horizontal Stabilizer-Elevator Combination 494 

Estimation of Aerodynamic Hinge Moments 495 

Example Calculation of Stick Force 502 

Stick-Free Longitudinal Static Stability 505 

Elevator-Stabilizer Configuration 505 

Stabilator Configuration 506 

Stick-Free Static Margin 508 

Steady Maneuvering 508 

Horizontal Stabilizer-Elevator Configuration: Elevator Angle per g 509 

Stabilator Angle per g 512 

Stick Force per g 512 

Stabilizer-Elevator Configuration 512 

Stabilator 513 

Effect of Fuselage and Nacelles 513 

Effects of Propulsion System 515 

Propellers 515 

Jets 52i 

Ground Effect 523 

Lateral and Directional Static Stability and Control 524 


CONTENTS 


XV 


Directional Static Stability 
Directional Control 
Lateral Control 
Adverse Yaw 

Roll Control by the Use of Spoilers 
Aileron Reversal 
Steady Rolling Motion 
Coupling Effects 
Rolling Moment with Rudder 
Rolling Moment with Yaw Rate 
Yawing Moment with Roll Rate 
Rolling Moment with Sideslip Angle—Dihedral Effect 

NINE LONGITUDINAL DYNAMIC STABILITY AND CONTROL 

Equations of Motion 
Linearization of the Equations 

A Summary Look at the Stability Derivatives and Other Parameters 
Affecting Longitudinal Dynamic Motion 
X Derivatives and Parameters 
Z Derivatives and Parameters 
M Derivatives and Parameters 

Examination and Reduction of Equations of Longitudinal Motion 
Solution of nth-Order Linear Differential Equations with Constant 
Coefficients 
Mode Shapes 

Phugoid (Long-Period Mode) 

Short-Period Mode 

Solution of the Longitudinal Equations of Motion Using an Analog 
Computer 

Longitudinal Flying Qualities 
Phugoid Mode 
Flight Path Stability 
Short Period Mode 

TEN LA TER A L -DIR EC TIONA L DYNAMIC STABILITY AND 
CONTROL 

Equations of Motion 

Euler Angles 

Reduction of the Lateral-Directional Equations of Motion 

A Summary Look at the Stability Derivatives and Other Parameters 
Affecting Lateral-Directional Dynamic Motion 
Y Derivatives 

C Yfl 

C Yp 

C Yt 

Cy s 


524 

526 

528 

534 

534 

537 

538 
541 

541 

542 

543 
543 

552 

552 


561 

561 

563 

564 

565 

565 

569 

574 


576 

585 

586 

587 
587 


591 

594 

595 

5% 

596 

597 

597 

598 
600 



CONTENTS 


l Derivatives 600 

C h 600 

0„ 601 

O f 602 

0* a 603 

O s ” 603 

N Derivatives 604 

Cn 0 604 

C Nf 604 

C Nf 604 

Cn 5j 605 

C Ntf 605 

Lateral-Directional Equations for the Cherokee 180 605 

Mode Shapes 606 

Roll Mode 606 

Spiral Mode 607 

Oscillatory or “Dutch Roll” Mode 609 

Lateral Directional Flying Qualities 610 

Roll Mode 611 

Dutch Roll Mode 612 

Spiral Stability 612 

Spinning 613 

APPENDIX A.1 THE SI SYSTEM 621 

APPENDIX A.2 STANDARD ATMOSPHERE 624 

APPENDIX A.3 AIRPLANE DATA 629 

APPENDIX A.4 NOMENCLATURE AND ABBREVIATIONS 640 


INDEX 


649 


AERODYNAMICS, AERONAUTICS, AND FLIGHT MECHANICS 



ONE 

INTRODUCTION 


Aeronautics is defined as “the science that treats of the operation of 
aircraft; also, the art or science of operating aircraft.” Basically, with 
aeronautics, one is concerned with predicting and controlling the forces and 
moments on an aircraft that is traveling through the atmosphere. 


A BRIEF HISTORY 


Thursday, December 17, 1903 

“When we got up a wind of between 20 and 25 miles was blowing from the 
north. We got the machine out early and put out the signal for the men at the 
station. Before we were quite ready, John T. Daniels, W. S. Dough, A. D. 
Etheridge, W. C. Brinkly of Manteo, and Johnny Moore of Nags Head arrived. 
After running the engine and propellers a few minutes to get them in working 
order, I got on the machine at 10:35 for the first trial. The wind, according to 
our anemometers at this time, was blowing a little over 20 miles (corrected) 27 
miles according to the government anemometer at Kitty Hawk. On slipping the 
rope the machine started off increasing in speed to probably 7 or 8 miles. The 
machine lifted from the truck just as it was entering the fourth rail. Mr. 
Daniels took a picture just as it left the tracks. I found the control of the front 
rudder quite difficult on account of its being balanced too near the center and 
thus had a tendency to turn itself when started so that the rudder was turned 
too far on one side and then too far on the other. As a result the machine 
would rise suddenly to about 10ft. and then as suddenly, on turning the 
rudder, dart for the ground. A sudden dart when out about 100 feet from the 
end of the tracks ended the flight. Time about 12 seconds (not known exactly 
as watch was not promptly stopped). The level for throwing off the engine was 
broken, and the skid under the rudder cracked. After repairs, at 20 min. after 
11 o'clock Will made the second trial." 

The above, taken from Orville Wright’s diary, as reported in Reference 
1.1, describes mankind’s first sustained, controlled, powered flight in a 




Figure 1.1 The first flight, December t7. 1903 (Courtesy of the National Air and 
Space Museum Smithsonian Institution) 


heavier-than-air machine. The photograph, mentioned by Orville Wright, is 
shown here as Figure 1.1. Three more flights were made that morning. The 
last one, by Wilbur Wright, began just at 12 o’clock and covered 260 m in 59 s. 
Shortly after this flight, a strong gust of wind struck the airplane, turning it 
over and over. Although the machine was severely damaged and never flew 
again, the Wright Brothers achieved their goal, begun approximately 4yr 
earlier. 

Their success was no stroke of luck. The Wright Brothers were painstak¬ 
ing in their research and confident of their own results. They built their own 
wind tunnel and tested, in a methodical manner, hundreds of different airfoil 
and wing planform shapes. They were anything but a “couple of bicycle 
mechanics.” Their letters to Octave Chanute, a respected civil engineer and 
aviation enthusiast of the day, reveal the Wright Brothers to have been 
learned men well versed in basic concepts such as work, energy, statics, and 
dynamics. A three-view drawing of their first airplane is presented in Figure 
1.2. 

On September 18, 1901, Wilbur Wright was invited to deliver a lecture 
before the Western Society of Engineers at a meeting in Chicago, Illinois. 
Among the conclusions reached by him in that paper were: 
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1. “That the ratio of drift to lift in well-shaped surfaces is less at angles of 
incidence of five degrees to 12 degrees than at an angle of three degrees” 
(“Drift” is what we now call “drag.”) 

2. “That in arched surfaces the center of pressure at 90 degrees is near the 
center of the surface, but moves slowly forward as the angle becomes less, 
till a critical angle varying with the shape and depth of the curve is reached, 
after which it moves rapidly toward the rear till the angle of no lift is 
found.” 

3. “That a pair of superposed, or tandem surfaces, has less lift in proportion 
to drift than either surface separately, even after making allowance for 
weight and head resistance of the connections.” 

These statements and other remarks (see Ref. 1.1) show that the Wright 
Brothers had a good understanding of wing and airfoil behavior well beyond 
that of other experimenters of the time. 

Following their first successful flights at Kitty Hawk. North Carolina, in 
1903, the Wright Brothers returned to their home in Dayton, Ohio. Two years 
later they were making flights there, almost routinely, in excess of 30 km and 
30 min while others were still trying to get off the ground. 

Most of the success of the Wright Brothers must be attributed to their 
own research, which utilized their wind tunnel and numerous experiments 
with controlled kites and gliders. However, their work was built, to some 
degree, on the gliding experiments of Otto Lilienthal and Octave Chanute. 
Beginning in 1891, Lilienthal, working near Berlin, Germany, made ap¬ 
proximately 2000 gliding flights over a 5-yr period. Based on measurements 
obtained from these experiments, he published tables of lift and drag 
measurements on which the Wright Brothers based their early designs. 
Unfortunately, Lilienthal had no means of providing direct aerodynamic 
control to his gliders and relied instead on kinesthetic control, whereby he 
shifted his weight fore and aft and side to side. On August 9, 1896, as the 
result of a gust, Otto Lilienthal lost control and crashed from an altitude of 
approximately 15 m. He died the next day. During 1896 and 1897, Octave 
Chanute, inspired by Lilienthal*s work, designed and built several gliders that 
were flown by others near Miller, Indiana. Chanute recognized LilienthaPs 
control problems and was attempting to achieve an “automatic” stability in 
his designs. Chanute*s principal contribution was the addition of both vertical 
and horizontal stabilizing tail surfaces. In addition, he went to the “box,” or 
biplane, configuration for added strength. Unfortunately, he also relied on 
kinesthetic control. 

When the Wright Brothers began their gliding experiments in the fall of 
1900, they realized that adequate control about all three axes was one of the 
major prerequisites to successful flight. To provide pitch control (i.e., nose up 
or down), they resorted to an all-movable horizontal tail mounted in front of 


the wing. Yaw control (i.e., turning to the left or right) was accomplished by 
means of an all-movable vertical tail mounted behind the wing. Their method 
of roll control (i.e., lowering one side of the wing and raising the other) was 
not as obvious from photographs as the controls about the other two axes. 
Here, the Wright Brothers devised a means of warping their “box” wing so 
that the angle of incidence was increased on one side and decreased on the 
other. The vertical tail, or rudder, was connected to the wing-warping wires so 
as to produce what pilots refer to today as a coordinated turn. 

The Wright Brothers were well ahead of all other aviation enthusiasts of 
their era. In fact, it was not until 3 yr after their first flight that a similar 
capability was demonstrated, this by Charles and Gabriel Voisin in Paris, 
France (Ref. 1,2). On March 30, 1907, Charles Voisin made a controlled flight 
of approximately 100 m in an airplane similar in appearance to the Wright 
flyer. A second machine built by the Voisin Brothers for Henri Farman, a 
bicycle and automobile racer, was flown by Farman later that year on flights 
that exceeded 2000 m. By the end of that year at least five others succeeded in 
following the Wright Brothers’ lead, and aviation was on its way. 

Today we are able to explain the results of the early experimenters in a 
very rational way by applying well-established aerodynamic principles that 
have evolved over the years from both analysis and experimentation. These 
developments have their beginnings with Sir Isaac Newton, who has been 
called the first real fluid mechanician (Ref. 1.3). In 1687 Newton, who is 
probably best known for his work in solid mechanics, reasoned that the 
resistance of a body moving through a fluid is proportional to the fluid 
density, the velocity squared, and the area of the body. 

Newton also postulated the shear force in a viscous fluid to be propor¬ 
tional to the velocity gradient. Today, any fluid obeying this relationship is 
referred to as a Newtonian fluid. 

In 1738, Daniel Bernoulli, a Swiss mathematician, published his treatise, 
“Hydrodynamics,” which was followed in 1743 by a similar work produced 
by his father, John Bernoulli. The Bernoulli made important contributions to 
understanding the behavior of fluids. In particular, John introduced the 
concept of internal pressure, and he was probably the first to apply momen¬ 
tum principles to infinitesimal fluid elements. 

Leonhard Euler, another Swiss mathematician, first put the science of 
hydrodynamics on a firm mathematical base. Around 1755, Euler properly 
formulated the equations of motion based on Newtonian mechanics and the 
works of John and Daniel Bernoulli. It was he who first derived along a 
streamline the relationship that we refer to today as “Bernoulli’s equation.” 

The aerodynamic theories of the 1800s and early 1900s developed from 
the early works of these mathematicians. In 1894 the English engineer, 
Frederick William Lanchester, developed a theory to predict the aerodynamic 
behavior of wings. Unfortunately, this work was not made generally known 


Table 1.1 Largest Aircraft Examples Starting with the Wright Brothers 









Table 1.1 (continued) 
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until 1907 in a book published by Lanchester. By then the Wright Brothers 
had been flying for 3 yr. Much of the knowledge that they had laboriously 
deduced from experiment could have been reasoned from Lanchester’s 
theory. In 1894, Lanchester completed an analysis of airplane stability that 
could also have been of value to the Wrights. Again, this work was not 
published until 1908. 

Lanchester’s wing theory was somewhat intuitive in its development. In 
1918 Ludwig Prandtl, a German professor of mechanics, presented a mathe¬ 
matical formulation of three-dimensional wing theory; today both men are 
credited with this accomplishment. Prandtl also made another important 
contribution to the science with his formalized boundary layer concept. 

Around 1917 Nikolai Ergorovich Joukowski (the spelling has been angli¬ 
cized), a Russian professor of rational mechanics and aerodynamics in 
Moscow, published a series of lectures on hydrodynamics in which the 
behavior of a family of airfoils was investigated analytically. 

The work of these early hydro- and aerodynamicists contributed little, if 
any, to the progress and ultimate success of those struggling to fly. However, 
it was the analytical base laid by Euler and those who followed him on which 
the rapid progress in aviation was built. 

After 1908, the list of aviators, engineers, and scientists contributing to 
the development of aviation grew rapidly. Quantum improvements were 
accomplished with the use of flaps, retractable gear, the cantilevered wing, 
all-metal construction, and the turbojet engine. This impressive growth is 
documented in Table 1.1. Note that in less than 10 yr from the Wright 
Brothers' first flight, the useful load increased from 667 N (1501b) to more 
than 13,300 N (3000 lb). In the next 10 yr, the useful load increased by a factor 
of 10 and today is more than 1.78 x 10* N (400,000 lb) for the Lockheed C-5A. 

Our state of knowledge is now such that one can predict with some 
certainty the performance of an airplane before it is ever flown. Where 
analytical or numerical techniques are insufficient, sophisticated experimental 
facilities are utilized to investigate areas such as high-lift devices, complicated 
three-dimensional flows in turbomachinery, and aerothermodynamics. 


A BRIEF INTRODUCTION TO THE TECHNOLOGY 
OF AERONAUTICS 

Consider the airplane in steady, climbing flight shown in Figure 1.3. The 
term steady means that the airplane is not accelerating; hence, all forces and 
moments on the aircraft must be in balance. To be more precise, one states 
that the vector sum of all forces and moments on the airplane must be zero. 
To depict the angles more clearly, all forces are shown acting through the 
center of gravity (eg). Although the resultant of all the forces must pass 
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Figure 1.3 Forces and moments on an airplane in a steady climb. 

through the center of gravity, it is not generally true that any one of the 
forces, with the exception of W, must satisfy this condition. 

In this figure V represents the velocity vector of the airplane’s center of 
gravity. This vector is shown inclined upward from the horizontal through the 
angle of climb, 0 C . The angle between the horizontal and the thrust line is 
denoted as 6. If this line is taken to be the reference line of the airplane, then 
one states that the airplane is pitched up through this angle. The angle 
between the reference line and the velocity vector, 6 - d c , is referred to as the 
angle of attack of the airplane. Later we will use other angles of attack 
referenced to the wing geometry; thus, one must be careful in interpreting lift 
and drag data presented as a function of the angle of attack. 

The thrust, T, is the propelling force that balances mainly the aerody¬ 
namic drag on the airplane. T can be produced by a propeller, a turbojet, or a 
rocket engine. 

The lift, L, is defined as the component of all aerodynamic forces 
generated by the aircraft in the direction normal to the velocity vector, V. In 
level flight this means principally the upward vertical force produced by the 
wing. Generally, however, it includes the tail and fuselage forces. For exam¬ 
ple, in landing many aircraft require a downward force on the horizontal tail 
in order to trim out the nose-down pitching moment produced by the wing 
flaps. This trimming force can be significant, requiring a wing lift noticeably in 
excess of the airplane’s weight. 

Similar to the lift, the drag, D, is defined as the component of all 
aerodynamic forces generated by the airplane in the direction opposite to the 
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velocity vector, V. This force is composed of two principal parts, the parasite 
drag and the induced drag. The induced drag is generated as a result of 
producing lift; the parasite drag is the drag of the fuselage, landing gear, 
struts, and other surfaces exposed to the air. There is a fine point concerning 
the drag of the wing to be mentioned here that will be elaborated on later. 
Part of the wing drag contributes to the parasite drag and is sometimes 
referred to as profile drag. The profile drag is closely equal to the drag of the 
wing at zero lift; however, it does increase with increasing lift. This increase 
is therefore usually included as part of the induced drag. In a strict sense this 
is incorrect, as will become clearer later on. 

W is the gross weight of the airplane and, by definition, acts at the center 
of gravity of the airplane and is directed vertically downward. It is composed 
of the empty weight of the airplane and its useful load. This latter weight 
includes the payload (passengers and cargo) and the fuel weight. 

The pitching moment, M, is defined as positive in the nose-up direction 
(clockwise in Figure 1.3) and results from the distribution of aerodynamic 
forces on the wing, tail, fuselage, engine nacelles, and other surfaces exposed 
to the flow. Obviously, if the airplane is in trim, the sum of these moments 
about the center of gravity must be zero. .. 

We know today that the aerodynamic forces on an airplane are the same 
whether we move the airplane through still air or fix the airplane and move 
the air past it. In other words, it is the relative motion between the air and 
airplane and not the absolute motion of either that determines the aerody¬ 
namic forces. This statement was not always so obvious. When he learned of 
the Wright Brothers’ wind tunnel tests. Octave Chanute wrote to them on 
October 12, 1901 (Ref. 1.1) and referred to “natural wind.” Chanute con-y 
jectured in his letter: 

“It seems to me that there may be a difference in the result whether the air 
is impinged upon by a moving body or whether the wind impinges upon the 
same body at rest . In the latter case each molecule , being driven from behind, 
tends to transfer more of its energy to the body than in the former case when 
the body meets each molecule successively before it has time to react on its 
neighbors.” 

Fortunately, Wilbur and Orville Wright chose to believe their own wind 
tunnel results. 

Returning to Figure 1.3, we may equate the vector sum of all forces to 
zero, since the airplane is in equilibrium. Hence, in the direction of flight, 

T cos ($ — 6 C ) — D — W sin 0 C = 0 (1.1) 

Normal to this direction, 

W cos 0 C — L — T sin (0 — 0 C ) = 0 


( 1 . 2 ) 
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These equations can be solved for the angle of climb to give 


6 C = tan 1 


T cos (e -0 C )-D 
L + T sin (8 — B c ) 


(1.3) 


In this form, 8 C appears on both sides of the equation. However, let us assume 
a priori that 8 C and (6 — 6 C ) are small angles. Also, except for very high 
performance and V/STOL (vertical or short takeoff and landing) airplanes, the 
thrust for most airplanes is only a fraction of the weight. Thus, Equation 1.3 
becomes 


T-D 


W 


(1.4) 


For airplanes propelled by turbojets or rockets, Equation 1.4 is in the 
form that one would normally use for calculating the angle of climb. 
However, in the case of airplanes with shaft engines, this equation is modified 
so that we can deal with power instead of thrust. 

First, consider a thrusting propeller that moves a distance S in time t at a 
constant velocity, V. The work that the propeller performs during this time is, 
obviously, 

work — TS 

Power is the rate at which work is performed; hence, 

T 5 

power ~T — 

But S/t is equal to the velocity of advance of the propeller. Hence the power 
available from the propeller is given by 

P avail = TV (1.5) 

Similarly, the power required by a body traveling through the air with a 
velocity of V and having a drag of D will be 

P req‘d = DV 

Thus, returning to Equation 1.4, by multiplying through by WV, we get 

W(V0 c ) = P avajJ -P rcq . d (1.6) 

The quantity V6 C is the vertical rate of climb, V c . The difference between the 
power that is required and that available is referred to as the excess power, 
Pxs. Thus Equation 1.6 shows that the vertical rate of climb can be obtained 
by equating the excess power to the power required to lift the airplane’s 
weight at the rate V c . In operating an airplane this means the following. A 
pilot is flying at a given speed in steady, level flight with the engine throttle 
only partially open. If the pilot advances the throttle while maintaining a 




constant airspeed, the power from the engine will then be in excess of that 
required for level flight, and the airplane will climb. 

Suppose, instead of keeping it constant, the pilot, while opening the 
throttle, allows the airspeed to increase in such a manner as to maintain a 
constant altitude. When a wide open throttle (WOT) condition is reached the 
maximum power available is equal to the power required. This is the con¬ 
dition for maximum airspeed, “straight and level.” 

From this brief introduction into airplane performance, it is obvious that 
we must be able to estimate the aerodynamic forces on the airplane before we 
can predict its performance. Also, a knowledge of the characteristics of its 
power plant-propulsor combination is essential. 

In addition to performance, the area of “flying qualities” is very im¬ 
portant to the acceptance of an airplane by the customer. Flying qualities 
refers primarily to stability and control, but it also encompasses airplane 
response to atmospheric disturbances. 

Let us briefly consider the pitching moment M, shown in Figure 1.3. This 
moment, which must be zero for steady, trimmed flight, results mainly from 
the lift on the wing and tail. In addition, contributions arise from the fuselage, 
nacelles, propulsor, and distribution of pressure over the wing. Suppose now 
that the airplane is trimmed in steady, level flight when it is suddenly 
disturbed (possibly by a gust or an input from the pilot) such that it pitches up 
by some amount. Before it can respond, the airplane’s path is still essentially 
horizontal, so that the angle between the velocity vector and the plane’s axis 
is momentarily increased. It will be shown later that, at a given airspeed, the 
moment, M, is dependent on this angle, defined previously as the angle of 
attack. Since the moment was initially zero before the airplane was disturbed, 
it follows that, in general, it will have some value other than zero due to the 
increase in angle of attack. Suppose this increment in M is positive. In this 
case the tendency would then be for the angle of attack to increase even 
further. This is an unstable situation where the airplane, when disturbed, 
tends to move even further from its steady-state condition. Thus, for the 
airplane to exhibit a more favorable, stable response, we desire that the 
increment in M caused by an angle of attack change be negative. 

This is about as far as we can go without considering in detail the 
generation of aerodynamic forces and moments on an airplane and its 
components. The preceding discussion has shown the importance of being 
able to predict these quantities from both performance and flying qualities 
viewpoints. The following chapters will present detailed analytical and 
experimental material sufficient to determine the performance and stability 
and control characteristics of an airplane. 

As you study the material to follow, keep in mind that it took the early 
aviation pioneers a lifetime to accumulate only a fraction of the knowledge 
that is yours to gain with a few months of study. 
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The primary system of units to be used 
Internationale) system. Smce lh ' b T dish s i te m is presented in Appendix 

United States, a compmson «o the Engtahsystem P number of 

the SI system, see Reference 1.4. 


problems 

,, Calculate the rate of climb 

o( 0.25 and a >“L«rsper Id. C^Tnt practice is to express 
fje of S clfmbinf«t pm — Wha, would be your answer tn these 

u of the systems (ball and track) pictured below are in equilibrium? 

Which are stable? 



ssss.' ’ITS £%Z —-■ - *— “ 

flight. If this angle is denoted as fc show that 


e c = tan' 


TcosJfr^D 
L+ T sin 0-T 


1.5 A student pushes against the shie of a budding with a force of 6N for a 
period of 4 hr. How much work ™> *>“• j, a , an altitud e of 1500 m 

1.6 An aircraft has a lift-to-drag ra ' 16km (994 miles) ahead. Will 

(4921 ft) when the engine fads. An airport is to Km t* 

the pdot be able to reach it? 
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TWO 

FLUID MECHANICS 


This chapter will stress the principles in fluid mechanics that are especi¬ 
ally important to the study of aerodynamics. For the reader whose pre¬ 
paration does not include fluid mechanics, the material in this chapter should 
be sufficient to understand the developments in succeeding chapters. For a 
more complete treatment, see any of the many available texts on fluid 
mechanics (e.g.. Refs. 2.1 and 2.2). 

Unlike solid mechanics, one normally deals with a continuous medium in 
the study of fluid mechanics. An airplane in flight does not experience a 
sudden change in the properties of the air surrounding it. The stream of water 
from a firehose exerts a steady force on the side of a burning building, unlike 
the impulse on a swinging bat as it connects with the discrete mass of the 
baseball. 

In solid mechanics, one is concerned with the behavior of a given, finite 
system of solid masses under the influence of force and moment vectors 
acting on the system. In fluid mechanics one generally deals not with a finite 
system, but with the flow of a continuous fluid mass under the influence of 
distributed pressures and shear stresses. 

The term fluid should not be confused with the term liquid, since the 
former includes not only the latter, but gases as well. Generally a fluid is 
defined as any substance that will readily deform under the influence of 
shearing forces. Thus a fluid is the antonym of a solid. Since both liquids and 
gases satisfy this definition, they are both known as fluids. A liquid is 
distinguished from a gas by the fact that the former is nearly incompressible. 
Unlike a gas, the volume of a given mass of liquid remains nearly constant, 
independent of the pressure imposed on the mass. 


FLUID STATICS AND THE ATMOSPHERE 

Before treating the more difficult case of a fluid in motion, let us consider 
a fluid at rest in static equilibrium. The mass per unit volume of a fluid is 
defined as the mass density, usually denoted by p. The mass density is a 
16 


constant for liquids, but it is a function of temperature, T, and pressure, p, for 
gases. Indeed, for a gas, p, p, and T are related by the equation of state 

p = pRT (2.1) 

R is referred to as the gas constant and has a value of 287.3 m 2 /°K-sec 2 for air 
at normal temperatures. In Equation 2.1, T is the thermodynamic or absolute 
temperature in degrees Kelvin. T and the Celsius temperature, t, are related 
by 

T = t + 273.15 (2.2) 

A container filled with a liquid is pictured in Figure 2.1a. A free-body 
diagram of a small slug of the fluid is shown in Figure 2.lb. This slug has a 
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unit cross-sectional area and a differential length of dh. Acting downward 
over the upper surface is the static pressure, p, while acting upward over the 
lower face is this same pressure plus the rate of increase of p with depth 
multiplied by the change in depth, dh. The static pressure also acts inward 
around the sides of the element, but this contributes nothing to the balance of 
forces in the vertical direction. In addition to the pressure forces, the weight 
of the fluid element, pgdh, acts vertically downward; g is the gravitational 
constant. 

Summing forces on the element in the vertical direction leads to 



Integrating Equation 2.3 from h = 0 at the surface to any depth, h, results 
in the static pressure as a function of the depth. 

P = Pa+pgh (2.4) 

where p a is the atmospheric pressure at the free surface. 

A manometer is a device frequently used to measure pressures. It is 
based on Equation 2.4. Consider the experimental setup pictured in Figure 
2.2. Here, a device known as a pitot-static tube is immersed in and aligned 
with a gas flow. The impact of the gas being brought to rest at the nose of the 
tube produces a pressure higher than that along the sides of the tube. This 
pressure, known as the total pressure, is transmitted through a tube to one 
side of a U-shaped glass tube partially filled with a liquid. Some distance back 
from the nose of the pitot-static tube the pressure is sampled through a small 
opening that is flush with the sides of the tube. This opening, if it is far 
enough back from the nose, does not disturb the flow so that the pressure 
sampled by it is the same as the static pressure of the undisturbed flow. This 
static pressure is transmitted to the right side of the glass U-tube manometer. 
The total pressure, being higher than the static pressure, causes the liquid in 
the left side of the U-tube to drop while the level on the right side rises. 

If we denote p as the static pressure and p + Ap as the total pressure, the 
pressure at the bottom of the U-tube can be calculated by Equation 2.4 using 
either the right or left side of the tube. Equating the results from the two sides 
gives 

p + Ap + pgh 0 -p + pg(Ah + h 0 ) 
or 

Ap = pg Ah (2.5) 

Hence, the difference of the liquid levels in the two sides of the manometer 
is a direct measure of the pressure difference applied across the manometer. 
In this case we could then determine the difference between the total pressure 
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Gas flow 



Figure 2.2 Pitot-static tube connected to a U-tube liquid manometer. 


and static pressure in the gas flow from which, as we will see later, the 
velocity of the gas can be calculated. 

Now consider the variation of static pressure through the atmosphere. 
Again the forces acting on a differential mass of gas will be treated in a manner 
similar to the development of Equation 2.3 for a liquid. However, h will be 
taken to be the altitude above the ground and, since the gravitational attrac¬ 
tion is now opposite to the direction of increasing h, the sign of Equation 2.3 
changes. For the atmosphere. 


dp 

dh 


= -PS 


( 2 . 6 ) 


The mass density, p, is not a constant in this case, so that Equation 2.6 
cannot be integrated immediately. In order to perform the integration the 
equation of state, Equation 2.1, is substituted for p, which leads to 

dp _ g dh 
p RT 


(2.7) 
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From experimental observation, the variation of temperature with al¬ 
titude is known or, at least, a standard variation has been agreed on. Up to an 
altitude of 11 km, the temperature is taken to decrease linearly with altitude at 
a rate, known as the lapse rate, of 6.51 °C/km. Thus, Equation 2.7 becomes 
dp _ dT 1 
p T R(dTldh) 
or 

S = 0 52561 ( 2 . 8 ) 

where 6 is the ratio of the static pressure at altitude to the pressure at sea 
level and 0 is the corresponding absolute temperature ratio. 

Using the equation of state, the corresponding density ratio, a, is 
obtained immediately from Equation 2.8. 

5 

a = e 

or 

a = 0 42561 (2.9) 

Using the standard lapse rate and a sea level temperature of 288.15°K, 6 as a 
function of altitude is given by 

6- 1— 0.02256 h (2.10) 

where h is the altitude in kilometers. 

The lower region of the atmosphere up to an altitude for which Equations 
2.8 to 2.10 hold is referred to as the troposphere. This is the region in which 
most of today’s flying is done. Above 11 km and up to an altitude of 
approximately 23 km, the temperature is nearly constant. This region forms 
the lower part of the stratosphere. Through the remainder of the stratosphere, 
the temperature increases, reaching approximately 270 °K at an altitude of 
around 50 km. 

Figure 2.3 presents graphs (taken from Ref. 2.3) of the various properties 
of the standard atmosphere as a function of altitude. Each property is 
presented as a ratio to its standard sea level value denoted by the subscript 
0. In addition to p, p, and T, the acoustic velocity and kinematic viscosity 
are presented. These two properties will be defined later. 

One normally thinks of altitude as the vertical distance of an airplane above 
the earth’s surface. However, the operation of an airplane depends on the 
properties of the air through which it is flying, not on the geometric height. 
Thus the altitude is frequently specified in terms of the standard atmosphere. 
Specifically, one refers to the pressure altitude or the density altitude as the 
height in the standard atmosphere corresponding t*.. the pressure or density, 
respectively, of the atmosphere in which the airplane is operating. An air¬ 
plane’s altimeter is simply an absolute pressure gage calibrated according to 
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Figure 2.3 The standard atmosphere. 

the standard atmosphere. It has a manual adjustment to allow for variations in 
sea level barometric pressure. When set to standard sea level pressure 
(760 mm Hg, 29.92 in. Hg), assuming the instrument and static pressure source 
to be free of errors, the altimeter will read the pressure altitude. When set to 
the local sea level barometric pressure (which the pilot can obtain over the 
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radio while in flight), the altimeter will read closely the true altitude above sea 
level. A pilot must refer to a chart prescribing the ground elevation above sea 
level in order to determine the height above the ground. 


FLUID DYNAMICS 

We will now treat a fluid that is moving so that, in addition to gravita¬ 
tional forces, inertial and shearing forces must be considered. 

A typical flow around a streamlined shape is pictured in Figure 2.4. Note 
that this figure is labled “two-dimensional flow”; this means simply that the 
flow field is a function only of two coordinates (x and y, in the case of Figure 
2.4) and does not depend on the third coordinate. For example, the flow of 
wind around a tall, cylindrical smokestack is essentially two-dimensional 
except near the top. Here the wind goes over as well as around the stack, and 
the flow is three-dimensional. As another example. Figure 2.4 might represent 
the flow around a long, streamlined strut such as the one that supports the 
wing of a high-wing airplane. The three-dimensional counterpart of this shape 
might be the blimp. 

Several features of flow around a body in general are noted in Figure 2.4. 
First, observe that the flow is illustrated by means of streamlines. A stream¬ 
line is an imaginary line characterizing the flow such that, at every point along 
the line, the velocity vector is tangent to the line. Thus, in two-dimensional 



Figure 2.4 Two-dimensional flow around a streamlined shape- 


flow, if y(x) defines the position of a streamline, y(x) is related to the x and y 
components of the velocity, u(x) and u(x), by 


dy v(x ) 
dx u(x) 


(2-11) 


Note that the body surface itself is a streamline. 

In three-dimensional flow a surface swept by streamlines is known as a 
stream surface. If such a surface is closed, it is known as a stream tube. 

The mass flow accelerates around the body as the result of a continuous 
distribution of pressure exerted on the fluid by the body. An equal and 
opposite reaction must occur on the body. This static pressure distribution, 
acting everywhere normal to the body’s surface, is pictured on the lower half 
of the body in Figure 2.4. The small arrows represent the local static pressure, 
p, relative to the static pressure, p 0 , in the fluid far removed from the body. 
Near the nose p is greater than p 0 ; further aft the pressure becomes negative 
relative to p 0 . If this static pressure distribution, acting normal to the surface, 
is known, forces on the body can be determined by integrating this pressure 
over its surface. 

In addition to the local static pressure, shearing stresses resulting from 
the fluid’s viscosity also give rise to body forces. As fluid passes over a solid 
surface, the fluid particles immediately in contact with the surface are brought 
to rest. Moving away from the surface, successive layers of fluid are slowed 
by the shearing stresses produced by the inner layers. (The term “layers” is 
used only as a convenience in describing the fluid behavior. The fluid shears 
in a continuous manner and not in discrete layers.) The result is a thin layer of 
slower moving fluid, known as the boundary layer, adjacent to the surface. 
Near the front of the body this layer is very thin, and the flow within it is 
smooth without any random or turbulent fluctuations. Here the fluid particles 
might be described as moving along in the layer on parallel planes, or laminae; 
hence the flow is referred to as laminar. 

At some distance back from the nose of the body, disturbances to the flow 
(e.g., from surface roughnesses) are no longer damped out. These disturbances 
suddenly amplify, and the laminar boundary layer undergoes transition to a 
turbulent boundary layer. This layer is considerably thicker than the laminar one 
and is characterized by a mean velocity profile on which small, randomly 
fluctuating velocity components are superimposed. These flow regions are 
shown in Figure 2.4. The boundary layers are pictured considerably thicker than 
they actually are for purposes of illustration. For example, on the wing of an 
airplane flying at 100 m/s at low altitude, the turbulent boundary 1.0 m back from 
the leading edge would be only approximately 1.6 cm thick. If the layer were still 
laminar at this point, its thickness would be approximately 0.2 cm. 

Returning to Figure 2.4, the turbulent boundary layer continues to 
thicken toward the rear of the body. Over this portion of the surface the fluid 
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is moving into a region of increasing static pressure that is tending to oppose 
the flow. The slower moving fluid in the boundary layer may be unable to 
overcome this adverse pressure gradient, so that at some point the flow 
actually separates from the body surface. Downstream of this separation 
point, reverse flow will be found along the surface with the static pressure 
nearly constant and equal to that at the point of separation. 

At some distance downstream of the body the separated flow closes, and 
a wake is formed. Here, a velocity deficiency representing a momentum loss 
by the fluid is found near the center of the wake. This decrement of 
momentum (more precisely, momentum flux) is a direct measure of the body 
drag (i.e., the force on the body in the direction of the free-stream velocity). 

The general flow pattern described thus far can vary, depending on the 
size and shape of the body, the magnitude of the free-stream velocity, and the 
properties of the fluid. Variations in these parameters can eliminate transition 
or separation or both. 

One might reasonably assume that the forces on a body moving through a 
fluid depend in some way on the mass density of the fluid, p, the size of the 
body, /, and the body’s velocity, V. If we assume that any one force, F, is 
proportional to the product of these parameters each raised to an unknown 
power, then 

F * p °v b l c 

In order for the basic units of mass, length, and time to be consistent, it 
follows that 



Considering M, L, and T in order leads to three equations for the unknown 
exponents a, b, and c from which it is found that a = 1, b = 2, and c = 2. 
Hence, 

F <x P V 2 / 2 (2.12) 

For a particular force the constant of proportionality in Equation 2.12 is 
referred to as a coefficient and is modified by the name of the force, for 
example, the lift coefficient. Thus the lift and drag forces, L and D, can be 
expressed as 

L = $pV 2 SC L (2.13 a) 

D = 2P V 2 S C d (2.13f») 

Note that the square of the characteristic length, I 2 , has been replaced by 
a reference area, S. Also, a factor of 1/2 has been introduced. This can be 
done, since the lift and drag coefficients, C L and Cd, are arbitrary at this 
point. The quantity pV 2 /2 is referred to as the dynamic pressure, the 
significance of which will be made clear shortly. 
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For many applications, the coefficients Cl and Cd remain constant for a 
given geometric shape over a wide range of operating conditions or body size. 
For example, a two-dimensional airfoil at a 1° angle of attack will have a lift 
coefficient of approximately 0.1 for velocities from a few meters per second 
up to 100 m/s or more. In addition, C L will be almost independent of the size 
of the airfoil. However, a more rigorous application of dimensional analysis 
[see Buckingham’s 7T theorem (Ref. 2.1)] will result in the constant of 
proportionality in Equation 2.12 possibly being dependent on a number of 
dimensionless parameters. Two of the most important of these are known as 
the Reynolds number, R, and the Mach number, M, defined by, 


R 


Vtp 

F 


(2.14a) 


where / is a characteristic length, V is the free-stream velocity, p. is the 
coefficient of viscosity, and a is the velocity of sound. The velocity of sound 
is the speed at which a small pressure disturbance is propagated through the 
fluid; at this point, it requires no further explanation. The coefficient of 
viscosity, however, is not as well known and will be elaborated on by 
reference to Figure 2.5. Here, the velocity profile is pictured in the boundary 
layer of a laminar, viscous flow over a surface. The viscous shearing produces 
a shearing stress of t w on the wall. This force per unit area is related to the 
gradient of the velocity u( y) at the wall by 



Actually, Equation 2.15 is applicable to calculating the shear stresses 
between fluid elements and is not restricted simply to the wall. Generally, the 
viscous shearing stress in the fluid in any plane parallel to the flow and away 
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from the wall is given by the product of p. and the velocity gradient normal to 
the direction of flow. 

The kinematic viscosity, v , is defined as the ratio of fx to p. 

V= IL 

P 

v is defined as a matter of convenience, since it is the ratio of p to p that 
governs the Reynolds number. The kinematic viscosity for the standard 
atmosphere is included in Figure 2.3 as an inverse fraction of the standard sea 
level value. 

A physical significance can be given to the Reynolds number by multi¬ 
plying numerator and denominator by V and dividing by /. 



In the following material (see Equation 2.28) the normal pressure will be 
shown to be proportional to pV 2 whereas, from Equation 2.15, fiV/l is 
proportional to the shearing stress. Hence for a given flow the Reynolds 
number is proportional to the ratio of normal pressures (inertia forces) to 
viscous shearing stresses. Thus, relatively speaking, a flow is less viscous than 
another flow if its Reynolds number is higher than that of the second flow. 

The Mach number determines to what extent fluid compressibility can be 
neglected (i.e., the variation of mass density with pressure). Current jet 
transports, for example, can cruise at Mach numbers up to approximately 0.8 
before significant compressibility effects are encountered. 

At lower Mach numbers, two flows are geometrically and dynamically 
similar if the Reynolds numbers are the same for both flows. Hence, for 
example, for a given shape, C D for a body 10 m long at lOOm/s will be the 
same as C D for a 100-m long body at 10 m/s. As another example, suppose 
transition occurs 2 m back from the leading edge of a flat plate aligned with a 
flow having a velocity of 50 m/s. Then, at 25 m/s transition would occur at a 
distance of 4 m from the leading edge. Obviously the effects of R and M on 
dimensionless aerodynamic coefficients must be considered when interpreting 
test results obtained with the use of small models. 

For many cases of interest to aerodynamics the pressure field around a 
shape can be calculated assuming the air to be inviscid and incompressible. 
Small corrections can then be made to the resulting solutions to account for 
these “real fluid” effects. Corrections for viscosity or compressibility will be 
considered as needed in the following chapters. 

Conservation of Mass 

Fluid passing through an area at a velocity of V has a mass flow rate 
equal to pAV. This is easily seen by reference to Figure 2.6. Here flow is 



T 


Figure 2.6 Mass flow through a surface. 

pictured along a streamtube of cross-sectional area A. The fluid velocity is 
equal to V. At time t = 0, picture a small slug of fluid of length, 1, about to 
cross a reference plane. At time //V, this entire slug will have passed through 
the reference plane. The volume of the slug is Al, so that a mass of pAI was 
transported across the reference plane during the time IIV. Hence the mass 
rate of flow, m, is given by 

_ pAI 
m (l/V) 

= pAV (2.16) 

Along a streamtube (which may be a conduit with solid walls) the quantity pAV 
must be a constant if mass is not to accumulate in the system. For incompressible 
flow, p is a constant, so that the conservation of mass leads to the continuity 
principle 

A V = constant 

AV is the volume flow rate and is sometimes referred to as the flux. 
Similarly, pAV is the mass flux. The mass flux through a surface multiplied by 
the velocity vector at the surface is defined as the momentum flux. Generally, 
if the velocity vector is not normal to the surface, the mass flux will be 

pAV ■ n 

with the momentum flux written as 

(pAV • n)V 

here n is the unit vector normal to the surface and in the direction in which 
the flux is defined to be positive. For example, if the surface encloses a 
volume and the net mass flux out of the volume is to be calculated, n would 
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be directed outward from the volume, and the following integral would be 
evaluated over the entire surface. 

J JpVndS 

Consider the conservation of mass applied to a differential control 
surface. For simplicity, a two-dimensional flow will be treated. A rectangular 
contour is shown in Figure 2.7. The flow passing through this element has 
velocity components of u and v in the center of the element in the x and y 
directions, respectively. The corresponding components on the right face of 
the element are found by expanding them in a Taylor series in x and y and 
dropping second-order and higher terms in Ax. Hence the mass flux out 
through the right face will be 

[ pu+ ^rT] iy 

Writing similar expressions for the other three faces leads to the net mass flux 
out being 

[^) + a^)] Aj . Ay 

The net mass flux out of the differential element must equal the rate at which 
the mass of the fluid contained within the element is decreasing, given by 

Ay) 


y 



Figure 2.7 A rectangular differential control surface. 
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Since Ax and Ay are arbitrary, it follows that, in general, 
dp d(pu ) | a(pu) _ Q 
dt dx dy 

In three dimensions the preceding equation can be written in vector notation 


where V is the vector operator, del, defined by 


.8 d J _ 1 d 
'Tx +l ^ +v Tz 


Any physically possible flow must satisfy Equation 2.17 at every point in the 
flow. 

For an incompressible flow, the mass density is a constant, so Equation 
2.17 reduces to 

VV = 0 (2.18) 

The above is known as the divergence of the velocity vector, div V. 


The Momentum Theorem 

The momentum theorem in fluid mechanics is the counterpart of New¬ 
ton’s second law of motion in solid mechanics, which states that a force 
imposed on a system produces a rate of change in the momentum of the 
system. The theorem can be easily derived by treating the fluid as a collection 
of fluid particles and applying the second law. The details of the derivation 
can be found in several texts (e.g.. Ref. 2.1) and will not be repeated here. 

Defining a control surface as an imaginary closed surface through which 
a flow is passing, the momentum theorem states: 

“The sum of external forces {or moments) acting on a control surface 
and internal forces {or moments) acting on the fluid within the control surface 
produces a change in the flux of momentum {or angular momentum) through 
the surface and an instantaneous rate of change of momentum (or angular 
momentum) of the fluid particles within the control surface. 

Mathematically, for linear motion of an inviscid fluid, the theorem can be 
expressed in vector notation by 

— j j pn dS + B = Is / pv(v * n > dS+ ii / J v / pv dr (2,19) 

In Equation 2.19, n is the unit normal directed outward from the surface. 
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S, enclosing the volume, V. V is the velocity vector, which generally depends 
on position and time. B represents the vector sum of all body forces within 
the control surface acting on the fluid, p is the mass density of the fluid 
defined as the mass per unit volume. 

For the angular momentum, 

Q -J JpOTX.XV.rtdS+i/JJp(VXr)<2r (2.20) 

Here, Q is the vector sum of all moments, both internal and external, acting 
on the control surface or the fluid within the surface, r is the radius vector to 
a fluid particle. 

As an example of the use of the momentum theorem, consider the force 
on the burning building produced by the firehose mentioned at the beginning 
of this chapter. Figure 2.8 illustrates a possible flow pattern, admittedly 
simplified. Suppose the nozzle has a diameter of 10 cm and water is issuing 
from the nozzle with a velocity of 60 m/s. The mass density of water is 
approximately 1000 kg/m 3 . The control surface is shown dotted. 
Equation 2.19 will now be written for this system in the x direction. 
Since the flow is steady, the partial derivative with respect to time of the 
volume integral given by the last term on the right side of the equation 
vanishes. Also, B is zero, since the control surface does not enclose any 



Figure 2.8 A jet of water impacting on a wall. 
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bodies. Thus Equation 2.19 becomes 

-jjpndS=jjp\(\-n)dS 

Measuring p relative to the atmospheric static pressure, p is zero 
everywhere along the control surface except at the wall. Here n is directed to 
the right so that the surface integral on the left becomes the total force 
exerted on the fluid by the pressure on the wall. If F represents the magnitude 
of the total force on the wall, 

-iF — J j pV(V ■ n) dS 

For the fluid entering the control surface on the left, 

V — 60! 
n = -i 

For the fluid leaving the control surface, the unit normal to this cylindrical 
surface has no component in the x direction. Hence, 

-iF = - J j (1000»0l(-60) dS 

--SSxlO'i J s JdS 

The surface integral reduces to the nozzle area of 7.85 x 10 -? m 2 . Thus, 
without actually determining the pressure distribution on the wall, the total 
force on the wall is found from the momentum theorem to equal 28.3 kN. 


Euler's Equation of Motion 

The principle of conservation of mass, applied to an elemental control 
surface, led to Equation 2.17, which must be satisfied everywhere in the flow. 
Similarly, the momentum theorem applied to the same element leads to 
another set of equations that must hold everywhere. 

Referring again to Figure 2.7, if p is the static pressure at the center of 
the element then, on the center of the right face, the static pressure will be 


P 


. dp Ax 
dx 2 


This pressure and a similar pressure on the left face produce a net force in the 
x direction equal to 
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Since there are no body forces present and the fluid is assumed inviscid, 
the above force must equaJ the net momentum flux out plus the instantaneous 
change of fluid momentum contained within the element. 

The momentum flux out of the right face in the x direction will be 
T . ■ d(pn) Ax~1/ ,5 m Ax\. 

Out of the upper face the corresponding momentum flux will be 




Similar expressions can be written for the momentum flux in through the left 
and bottom faces. 

The instantaneous change of the fluid momentum contained within the 
element in the x direction is simply 


Thus, equating the net forces in the x direction to the change in 
momentum and momentum flux and using Equation 2.17 leads to 

dt dx dy p dx K ’ 

Generalizing this to three dimensions results in a set of equations known as 
Euler’s equations of motion. 


Notice that if u is written as u(x, y, z, t), the left side of Equation 2.22 is 
the total derivative of u. The operator, d( )jdt, is the local acceleration and 
exists only if the flow is unsteady. 

In vector notation Euler’s equation can be written 


If the vector product of the operator V is taken with each term 
Equation 2.23, Equation 2.24 results. 
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> is the curl of the velocity vector, V x V, and is known as the vorticity. 


One can conclude from Equation 2.24 that, for an inviscid fluid, the 
vorticity is constant along a streamline. Since, far removed from a body, the 
flow is usually taken to be uniform, the vorticity at that location is zero; 
hence, it is zero everywhere. 

Bernoulli’s Equation 

Bernoulli’s equation is well known in fluid mechanics and relates the 
pressure to the velocity along a streamline in an inviscid, incompressible flow. 
It was first formulated by Euler in the middle 1700s. The derivation of this 
equation follows from Euler’s equations using the fact that along a streamline 
the velocity vector is tangential to the streamline. 


First, multiply Equation 2.22a through by dx and then substitute Equation 
2.26 for v dx and w dx. Also, the first term of the equation will be set equal to 
zero; that is, at this time only steady flow will be considered. 


Similarly, multiply Equation 2.22b by dy. Equation 2.22c by dz, and 
substitute Equation 2.26 for u dy, w dy and u dz, v dz, respectively. Adding the 
three equations results in perfect differentials for p and V 2 , V being the 
magnitude of the resultant velocity along the streamline. This last term results 
from the fact that 

— - * 

U dx 2 dx 
and 


Thus, along a streamline, Euler’s equations become 
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If p is not a function of p (i.e., the flow is incompressible). Equation 2.27 can 
be integrated immediately to give 

P + ipV 2 — constant (2.28) 

If the flow is uniform at infinity, Equation 2.28 becomes 

P + 2 pV 2 = constant = p *+%pVj (2.29) 

Here V is the magnitude of the local velocity and p is the local static 
pressure. V« and p* are the corresponding free-stream values. Equation 2.29 
is known as Bernoulli’s equation. 

The counterpart to Equation 2.29 for compressible flow is obtained by 
assuming pressure and density changes to follow an isentropic process. For 
such a process, 

p/p v = constant (2.30) 

y is the ratio of the specific heat at constant pressure to the specific heat at 
constant volume and is equal approximately to 1.4 for air. Substituting 
Equation 2.30 into Equation 2.27 and integrating leads to an equation some¬ 
times referred to as the compressible Bernoulli’s equation. 

— = constant (2.31) 

2 y - 1 p 

This equation can be written in terms of the acoustic velocity. First it is 
necessary to derive the acoustic velocity, which can be done by the use of the 
momentum theorem and continuity. Figure 2.9 assumes the possibility of a 
stationary disturbance in a steady flow across which the pressure, density, 
and velocity change by small increments. In the absence of body forces and 
viscosity, the momentum theorem gives 


— dp = (p + dp)(u + du) 2 — pu 2 

But, from continuity, 

(p + dp)(u + du) = pu 



Figure 2.9 A stationary small disturbance in a steady compressible flow. 
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If the small disturbance is stationary in the steady flow having a velocity of u , 
then obviously u is the velocity of the disturbance relative to the fluid. By 
definition, it follows that u, given by Equation 2.32, is the acoustic velocity. 
By the use of Equation 2.30, the acoustic velocity is obtained as 


An alternate form, using the equation of state (Equation 2.1), is 
a = ( yRT) ]12 

Thus Equation 2.31 can be written 

V 2 a 2 

-z- H- t = constant 

2 y- 1 


(2.34) 

(2.35) 


The acoustic velocity is also included in Figure 2.3 for the standard atmos¬ 
phere. 


Determination of Free-Stream Velocity 

At low speeds (compared to the acoustic velocity) a gas flow is essen¬ 
tially incompressible. In this case, and for that of a liquid. Equation 2.29 
applies. If the fluid is brought to rest so that the local velocity is zero then, 
from Equation 2.29, the local pressure, referred to in this case as the 
stagnation or total pressure, is equal to the sum of the free-stream static 
pressure, p«, and pVj/2. This latter term is called the dynamic pressure and is 
frequently denoted by the symbol q. Thus, 

V. = p(Po-P»)j (2.36) 

where p 0 is the total pressure, also referred to as the stagnation or reservoir 
pressure. The pitot-static tube shown in Figure 2.2 measures (po-p«) and is 
probably the most common means used to determine airspeed. However, 
notice that Equation 2.36 contains the mass density that must be determined 
before the airspeed can be calculated. This is most readily achieved by 
measuring, in addition to the difference between the stagnation pressure and 
the static pressure, the static pressure itself and the temperature. The density 
is then determined from the equation of state (Equation 2.1). 

At higher speeds (and we will now examine what is high) Equation 2.29 
no longer holds, so that V* must be determined from Equation 2.31. 

v 1 , 7 p y Po 

2 y - 1 p y - 1 po 
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At this point the subscript <» has been dropped, so that p, p, and V without a 
subscript may be local or free-stream values. 

Remembering that yplp is the square of the acoustic velocity, the 
preceding equation becomes 

M 2 { 1 _ 1 po p 

2 y - 1 y - 1 p po 
Using Equation 2.3, this can be written as 

(2.37) 

The dynamic pressure, q, is defined as 

q = 1 ipV 2 (2.38) 

which can be written in terms of the Mach number as 
q - \ypM- 

Combining this with Equation 2.37 gives 


P«~P _ 2 
q yM‘ 


(2.39) 


The square root of Equation 2.39 is presented graphically in Figure 2.10. 
The departure of this function from unity is a measure of the error to be 
incurred in calculating the airspeed from the incompressible Bernoulli equa¬ 
tion. Below a Mach number of 0.5 the error is seen to be less than 3%. 


Determination of True Airspeed 

During training, a pilot soon learns that the airspeed that appears on the 
airspeed indicator is not the true airspeed. Instead, in order to determine the 
true airspeed, the pilot must also read the altimeter and outside air tem¬ 
perature. The pilot then resorts to a small hand calculator or, in some 
instances, adjusts the dial on the airspeed indicator accordingly to allow for 
the atmospheric properties. 

The airspeed indicator is nothing more than an accurate differential pres¬ 
sure gage calibrated according to Equation 2.31. This equation can be put in 
the form 

v ^[(^-r-i] 

The airspeed indicator measures the difference, p« — p (sometimes called the 
compressible dynamic pressure), but then is calibrated to obtain V by 
assuming standard sea level values for the acoustic velocity and the free- 


0 0.2 0.4 0.6 0.8 1.0 

Mach number, M 

Figure 2.10 Relationship between reservoir pressure and dynamic pressure as a 
function of Mach number. 



and can be obtained from Equation 2.40. 0 and 8 are the temperature and 
pressure ratios, respectively. 

If compressibility can be neglected, the airspeed will be given by Equa- 
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tion 2.36. In this case, using the standard sea level value of p to calibrate the 
airspeed indicator dial leads to the equivalent airspeed. 

v J2(P,r_P)] m 
' L Psl J 



where cr is the density ratio. 

Finally, the indicated airspeed, V, is defined simply as the airspeed that 
the pilot reads from the airspeed indicator. If the indicator has no mechanical 
error (instrument error ) and if the static source is located so that it samples 
the true free-stream static pressure (otherwise a position error incurs), 
Vi = Vcai. Furthermore, if the Mach number is not too high, Vi = V ca j = V f . 

As an example in determining true airspeed, suppose a pilot reads an 
indicated airspeed of 180 m/s for an OAT of 239 °K and an altimeter reading 
of 6000 m. This altitude, according to Equations 2.8 and 2.10, corresponds to a 
pressure ratio 5 of 0.466. The measured temperature ratio, 6, is equal to 0.826. 
Hence cr = 0.564. According to Equation 3.42, the true airspeed will be 
239.7 m/s. From Equation 2.41, the true airspeed is calculated to be 231.6 m/s. 
Thus, using the incompressible relationship to estimate the true airspeed from 
the indicated airspeed results in a speed a few percent higher than that 
obtained from the calibrated airspeed relationship. 

To be precise one calculates the true airspeed from the calibrated 
airspeed and then determines the equivalent airspeed from its definition. 
Equation 2.42. In the previous example, this results in a V e of 173.9 m/s, a 
value 3.4% less than the calibrated airspeed. 


POTENTIAL FLOW 

For a steady, inviscid, incompressible flow, Euler’s equations of fluid 
motion reduce to two relatively simple relationships that govern the velocity 
vector. 

div V = V • V = 0 (2.43a) 

curl V = V x V = 0 (2.43b) 

The first equation satisfies conservation of mass; the second one assures 
that the dynamics of the flow is treated correctly. 

In addition to satisfying Equation 2.43 one must assure that any mathe¬ 
matical description of the flow field around a given body shape satisfies the 
boundary condition that there be no velocity normal to the body at all points 
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on its surface. If 
hold. 


is the unit vector normal to the surface, the following must 
V ■ n = 0 (2.44) 


Velocity Potential and Stream Function 


To assist in the solution of Equation 2.43, two functions are introduced. 
The first of these is known as the velocity potential, <£, and is defined such 
that 


u 


_ d<£ 
dx 


or, generally, 



V= V</> 


(2.45) 


Equation 2.43 satisfies identically Equation 2.43b. However, in order to 
satisfy Equation 2.43a, it follows that </> must be a harmonic function; that is, 
V 2 </> = 0 (2.46) 


The operator, V 2 , known as the Laplacian, is defined as 


a 2 a 2 _a^_ 
ax 1 + ay 2 + a? 


A flow for which Equation 2.43 is satisfied, and hence <f> can be defined, 
is known as a potential flow. The resulting fluid motion is described as being 
irrotational. This follows since, in the limit at a point, the curl of the velocity 
vector, which is zero, is equal to twice the rotational or angular velocity. 

The stream function, ifr, is related to the velocity components by 



(2.47) 


can only be defined for two-dimensional, or axisymmetric, flow. To obtain a 
particular component, the partial derivative of i (r is taken in the direction 
normal to the velocity and to the left as one looks in the direction of the 
velocity. 

A line element is pictured in Figure 2.11 with flow passing through it. This 
element is a segment of an arbitrary line connecting two points A and B. The 
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differential flux through this element will be 

dQ = V • n ds 
But 

n = (j dy - j dx)lds 
V = i m + j t> 

Substituting n and V into dQ and using Equation 2.47 results i 


That is, the change in the stream function between two points is equal to the 
flux between the points. It follows that i/> is a constant along a streamline. 
This can be shown by noting that along a streamline 


dtp = dQ = udy-vdx 
Combined, the two relationships give 

dip = 0 
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or 

i/f = constant (along a streamline) 

The stream function, as a measure of the flux, satisfies identically 
Equation 2.43a. For an irrotational flow, however, in order to meet Equation 
2.43b it follows that ip must also be harmonic. 

? V = 0 (2.49) 

In a manner similar to the derivation of Equation 2.48, the change in <p 
between two points can also be easily obtained. If 

<t> = 4>{x, y) 

then 

d * = to ix+ i$ d > 

= udx + vdy 

or, using vector notation, 

<£(B) - tf>(/\) = J V-dR (2.50) 

where R is the radius vector to the curve along which the integration is being 
performed, as shown in Figure 2.11. dR is then the differential vector along 
the curve, directed positively, with a magnitude of ds. 

As an example in the use of <p and ip, consider the uniform flow pictured 
in Figure 2.12. For this flow, 

u = U = constant 
v =0 

ip will be taken to be zero along the x-axis. This choice is arbitrary, since 
the values of both tp and \p can be changed by a constant amount without 
affecting the velocity field obtained from their derivatives. Equation 2.48 will 
be zero if the integral is performed along a line for which y is a constant. 



Figure 2.12 Uniform flow in the x direction. 
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Thus ip changes only in the y direction. Integrating Equation 2.48 in this 
direction gives 

iHy) - >K0) =J (iU) i dy 
or 

ip=Uy 

If the uniform flow contains, in addition to U, a constant velocity 
component V in the y direction, ip becomes 

tp = Uy-Vx (2.51) 

The minus sign in Equation 2.51 is in accordance with the positive 
direction of n, as shown in figure 2.11. n is directed to the right as one looks in 
the direction of point B from point A. In this case the integration for the 
second term is in the positive x direction, so n equals -j. 

In a more formal manner, 4> will be derived for this same flow using 
Equation 2.50. 

V=iU+jV 
R = ix + j y 

Hence 

dR= i dx + j dy 

so that, taking <p = 0 at x, y — 0, 

<f>(x, y) = f Udx + Vdy 
J o 

= Ux + Vy (2.52) 

Elementary Flow Functions 

If and <f>2 are functions satisfying Equation 2.50 then, because this 
equation is linear, their sum will also satisfy Equation 2.50. In general, both 
the velocity potential and stream function can be constructed by summing 
less complicated functions. 

<P (x, y) = <pi(x, y) (2.53a) 

y) = 2 y) (2.53b) 

Equation 2.53 represents the real benefit to be gained in describing a flow 
in terms of <p and ip. This statement will become obvious as the developments 
proceed. 

The simple flows from which more complicated patterns can be 
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developed are referred to as elementary flow functions. There are three of 
them: uniform rectilinear flow, vortex, and source. The first of these has 
already been covered with <f> and ip given by Equations 2.52 and 2.51, 
respectively. 


Vortex 


A vortex is pictured in Figure 2.13. This flow in two dimensions is purely 
circular around a point with no radial velocity component. Denoting the 
tangential velocity component by v e , the problem is to find v e as a function of 
r that will satisfy the set of Equations 2.42a and 2.42h. is to be independent 
of 6. 

In polar coordinates, 


dv e | V 8 1 dv r 
' dr + r r d$ 


V. v='^+^+^ 
r 36 dr r 


where r and 6 are the polar coordinates, with v, being the radial component of 
velocity and v e the tangential component. 

Since v r is zero in Figure 2.13 and v e is independent of 0, Equation 2.54b 
is satisfied identically and, from Equation 2.54a 

<fc, +Ii = 0 

dr r 

or, after integrating, 

rv g = constant 
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Thus, for potential flow, the tangential velocity around a vortex must 
vary inversely with the radial distance from the center of the vortex. 

The strength of a vortex, denoted by y, is measured by integrating the 
tangential velocity completely around the vortex. The value of this integral is 
independent of the path providing it encloses the singular point at the center 
of the vortex. 

y = j> V • dR (2.55) 

This closed-line integral of the velocity is known as the circulation. 
Evaluating Equation 2.55 on a constant radius leads to the relationship 
between the tangential velocity around a vortex, the radius, and the vortex 
strength. 


Equation 2.55 is a well-known relationship and can be easily remembered 
from the definition of -y. 

(ft and for a vortex follow immediately from Equations 2.50, 2.48, and 

2.56. 


<HB)-<HA) = f"^irde) 


If 0 is measured relative to zero and 0(B) is taken to be any value of 9, 


The stream function for a vortex is found from 




Letting i }/(A) be zero and r(A ) be an arbitrary radius, a, leads to 


T 2 p a K 

The minus sign results from the choice of positive coordinate directions. 


The source is the counterpart of a vortex. Here the flow, pictured in 
Figure 2.14, is again symmetrical about the center, but it is entirely radial with 
no tangential velocity component. 
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Equations 2.56 and 2.59, which define the velocities around vortices and 
sources, can be extended to three dimensions. If Q is the strength of a 
three-dimensional source, this flux will equal the product of the radial velocity 
and the surface area through which the velocity is passing. Thus, one can 
write a r immediately as 


Biot~Savart Law 

The three-dimensional velocity field associated with a vortex line is 
considerably more complicated and is given by the Biot-Savart law. The 
derivation of this law is beyond the scope of this text. Figure 2.15a illustrates 




^ X 

Figure 2.15b The Biot-Savart law for a straight-line vortex. 


a portion of a vortex line about which at any point the circulation, F\ is 
constant. If v, is the velocity vector induced at any point, P, in the field by the 
vortex line, the Biot-Savart law states 


dvi 


J_rxdR 
4it frp 


(2.63) 


This is the most general form of the Biot-Savart law. dR is the derivative 
of the radius vector from the origin to the vortex line and is thus the directed 
differential distance along the line, r is the radius vector from the point P to 
the line element dR. The positive direction of the circulatory strength, T, is 
defined according to the right-hand rule. The x, y , z orthogonal coordinate 
system is also right-handed. 

A special form of the Biot-Savart law for a straight-line vortex segment 
found in many texts can be obtained by integrating Equation 2.63. Referring 
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to Figure 2A5b, for convenience the line vortex is placed on the x-axis and 
lies between 0 and x. The z-axis will project out of the paper according to the 
right-hand rule. The circulation Y is taken to be positive in the x direction 
which means it will be clockwise when viewed in that direction. For this 
figure, 

R = ix 
OP + r = R 

OP = i x p + j y p 

Thus 



a, p, and h are defined in Figure 2A5b. Notice that the velocity has only a z 
component. As the line becomes infinite in length, the angles a and p 
approach zero, and Equation 2.64 reduces to the expression for the velocity 
around a two-dimensional point vortex given by Equation 2.56. 


The Calculation of Flows for Well-Defined Body Shapes 

The flow functions described thus far are basic functions. By combining 
these functions a multitude of more complicated flows can be described. 
When combining these functions, the velocities will add vectorially. This is 
obvious from Equation 2.53 a, since 

grad <f> = grad <£, + grad <f> 2 + • • • 
or 

v = v, + v 2 +- • • 

As an example of the use of these functions, consider the classic 
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Figure 2.16 Source-sink combination in a uniform flow. 

two-dimensional case illustrated in Figure 2.16. Here a source and a sink (a 
negative source) of equal strength are placed a unit distance from either side 
of the origin in a uniform flow. Three velocities from the source, sink, and 
uniform flow are shown added vectorially at P giving the resultant velocity, 
V 

It is easily verified that the entire resulting flow is symmetrical about both 
the x- and y-axes. To the left of the source along the x-axis a distance of x„ 
from the origin a stagnation point exists where the resultant velocity is zero. 
The location of this point can be found from Equation 2.59 for the velocity 
from a source (or sink). 

0 — u _ —5-i-—- 

U U 2tt(x 0 -1) 2tt(x 0 +1) 

or 

X ,i =1 + JL (2.65) 

It will be seen that this point lies on a dividing streamline that is closed 
and that separates the flow leaving the source and entering the sink from the 
uniform flow. The resultant streamline pattern can be constructed by cal¬ 
culating the velocity at many points in the field and fairing streamlines tangent 
to these vectors. However, a more direct way is to form the stream function 
and then solve for y(x) for constant values of </». Adding the iff functions for 
the uniform flow, source, and sink, one obtains 

iff = ^(6 l -0 2 )+Uy ( 2 - 66 ) 

where and 0 2 are shown in Figure 2.16. Because of the multivaluedness of 
the tangent function, one must be careful in evaluating this expression. 
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At the stagnation point, $i = 0 2 and y = 0, so that </» = 0- Since this point 
lies on the dividing streamline, a value of if/ = 0 will define this streamline. 
Hence 


Since the flow is symmetrical, we only need to calculate the streamline 
shapes in one quadrant, say the one for which x and y are both positive. In 
this quadrant. 


Hence x and y along the dividing streamline are related by 


r ‘r^ 7 "'( , - 2 7 3 ’) 


This can be solved explicitly for x. 

X ' = , ' y!_ tan{ 1 rll ^VUylq)]} (2 ' 69) 

Notice that as y approaches zero in Equation 2.68, x approaches the 
value given by Equation 2.65. 

Since any streamline can be replaced by a solid boundary, the dividing 
streamline represents a closed body. If it is assumed that the body’s fineness 
ratio (i.e., thickness to length) is small, then at its maximum thickness point 
(x = 0), Equation 2.68 can be written approximately as 


(l-2^f*) 


yo is the semithickness of the body corresponding to x 0 , the semilength. 
Hence the fineness ratio of the body, tfl, is related to ql U by 

- =-- m (2.71) 

/ 2[1 + (q/7rf/)] 32 1 

This classical body shape is referred to as the Rankine oval. The streamline 
pattern can be determined as a function of the streamline position far from 
the body. If y„ is the location of a particular streamline away from the body 
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then, for this particular streamline, 

if/ = y»U 

Equating this to Equation 2.66, one obtains a relationship between x and y 
along the streamline as a function of y * and qlU that can be solved explicitly 
for x. 

x = {' " y2 ~ tan u[l - 2 (UlqKy - yj]} <2 ' 72) 

This relation was used to calculate the streamline patterns of Figure 2.17. 
Only the flow external to the dividing streamline is shown. 

For comparison, this figure presents the streamline patterns (in one 
quadrant only) for 20% and 50% thick ovals. For each thickness, qlU was 
chosen on the basis of Equation 2.69. Because this equation assumes til to be 
small, the thickness ratio of the one shape is slightly greater than 50%. 

Having defined the shape of a body and the velocity field around it, the 
next obvious point of interest is the pressure distribution over the body’s 
surface. In order to remove the dependence of the predicted pressure dis¬ 
tribution on the free-stream pressure and velocity, the pressure distribution is 
normally presented in coefficient form, C p , where C p is defined according to 



Figure 2.17 Calculated streamlines for 20 and 50% thick Rankine ovals. 
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From Equation 2.29, C p is found from the ratio of local velocity to free- 
stream velocity. 

C '=‘-(vS (2-74) 

Returning now to the Rankine oval, note first that C p = 1 at the stagnation 
point where V is zero. Moving away from the nose along the if/ = 0 stream¬ 
line, the velocity increases to a maximum at some location. Depending on the 
fineness ratio, this point of maximum velocity may or may not be at the 
maximum thickness of the body; that is, x = 0. 

Although one could work with </r, knowing the source (and sink) strength 
the easiest approach is to calculate the u and v components directly by adding 
the components attributed to each elementary flow function. In this case it 
will be found that 

— = 1 i q f x+l *-1 1 
U 0 2ttU 0 L(x + 1) 2 + y 2 (x — 1) 2 + y 2 J 

v _ g r y 

U 0 27 tU 0 l(x + l) 2 + y : 

The pressure coefficient is then calculated from 

c ° = l -[(v-) 2+ (fr$] < 2 - 75 > 

The pressure distribution along the surfaces of 20 and 50% thick Rankine 
ovals have been calculated using the preceding equations, and the results are 
presented in Figure 2.18. It is not too surprising to find that, for the 20% oval, 
the minimum pressure occurs near the nose, where the curvature is the 
greatest. For the 50% thick oval the minimum pressure occurs approximately 
halfway from the center of the body to the nose, but it is nearly flat over the 
middle 70% of the body’s length. 


The Circular Cylinder 

The flow field around a circular cylinder and resulting pressure dis¬ 
tribution can be determined as a limiting case of the Rankine oval. In Figure 
2.16 the source and sink were placed at x = - 1.0 and x = 1.0, respectively. 
We will now move them instead toward the origin, but increase their strengths 
in inverse proportion to the distance between them. In the limit as the 
distance between the source and sink goes to zero, a so-called source-sink 
doublet is obtained. 

Letting 2e equal the distance between the source and sink and m the 
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The tangential velocity along the surface of the cylinder is found from 
Equation 2.47 by differentiating ijf with respect to r and evaluating the result 
at r = R. In this way v e is found to be 

v„ = 2U sin 6 (2.77) 

The pressure coefficient distribution, from Equation 2.45, is thus predic¬ 
ted to be 

C p = 1 - 4 sin 2 0 (2.78) 

In Chapter Four it will be seen that Equation 2.78 agrees fairly well with 
experimental results over the front half of the cylinder, but departs from 
actual measurements over the rear portion as the result of viscosity. 

A point vortex of strength y can be placed at the origin without altering 
the streamline representing the surface of the cylinder. If this is done, 
Equation 2.77 becomes 

v g = 2U sin 0 + (2.79) 

Relative to p a the pressure on the surface of the cylinder will be 

p -p n = |pf/ 2 -|p ^2U sin 0 + 2 ^] (2.80) 

Referring to Figure 2.19, the net vertical force, or lift, on the cylinder 
resulting from the pressure distribution will be 



figure 2.19 Circular cylinder with circulation. 
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or, from Equation 2.78, this reduces to 

L = pUy (2.81) 

This is referred to as the Kutta-Joukowski law. Although derived here 
specifically for a circular cylinder, it can be applied to other shapes where y 
represents generally the circulation around the shape. This will be amplified 
further in Chapter Three. 

The net horizontal force, or drag, on the cylinder is found from 
D = j pR cos 6 d0 

Using Equation 2.80, the drag is found to be zero, a result that is true in 
general for a closed body in steady potential flow. This result is known as 
D’Alembert’s paradox, after Jean le Rond D’Alembert, a French mathemati¬ 
cian who first reached this conclusion around 1743. 

The Numerical Calculation of Potential Flow Around 
Arbitrary Body Shapes 

The preceding has demonstrated how particular body shapes can be 
generated by the superposition of elementary flow functions. This procedure 
can be generalized and the inverse problem can be solved where the body 
shape is prescribed and the elementary flow functions that will generate the 
body shape are found. 

The concept of a point source or a point vortex can be extended to a 
continuous distribution of these functions. Consider first the two-dimensional 
source distribution shown in Figure 2.20. Here q is the source strength per 
unit length. 

Consider the closed contour shown dashed in Figure 2.20, having a length 
of Ax and a vanishing small height. The total flux through this surface must 
equal q Ax. Close to the surface the u velocity components from the elemen¬ 
tal sources will cancel so that only a v component remains. Thus 
2v Ax = q Ax 
or 



In Reference 2.4 this relationship is used to determine the flow about 
arbitrary shapes. Thus, unlike the Rankine oval, the body shape is specified, 
and the problem is to find the distribution of singularities to satisfy the 
condition that the velocity everywhere normal to the body surface to be zero. 
This particular problem is referred to as the Neumann problem. Essentially 
the numerical solution of the problem proceeds by segmenting the body 
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Figure 2.20 Distributed sources in two-dimensional flow. 


surface and distributing a unit source strength, q„ over the ith element. The 
normal velocity induced at the middle of the ith element by q, is obtained 
immediately from Equation 2.82. The contribution to the velocity at the ith 
element from another element is calculated by assuming the total source 
strength at the second element to be a point source located at the middle of 
that element. Taking n elements and letting i = l, 2,3 ,,n leads to a set of 
n linear simultaneous algebraic equations for the unknowns, q„ q 2 , Qi,..., q„. 

Consider in more detail this approach for two-dimensional flow. Figure 
2.21 shows two elements along the surface of a body. The ith element is the 
control element over which the unit source strength q f is distributed. At the 
Jth element a point source is located having a strength equal to q, AS f , ASj 
being the length of the yth element. The free-stream velocity U 0 is shown 
relative to the body x-axis at an angle of attack of a. 

At the center of the ith element the normal velocity components from 
each source and the free stream must vanish. Hence 

U 0 sin (0i - a ) = y + 2 qf u j * i (2.83) 

Cn is an influence coefficient, which accounts for the geometry of the 
body shape in determining the normal velocity induced at the ith element by 
the source at the /th element. 

If £ and rj; correspond to the midpoints of the ith element then, for Figure 

2 . 21 : 


(ft ~ <f>u) 

l-nfij 


(2.8 4a) 
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y 



Figure 2.21 Numerical solution of the Neumann problem by a distribution of 
sources. 



(2.84 b) 


(2.8 Ac) 

e , = tan -.|±iri 

(2.8 4d) 

Ci = j (Xi + x, +1 ) 

(2.8 4e) 

1, 

’)i = 2'K + Jiu) 

(2.84/) 


Having thus determined the source strengths, q if the resultant velocity at 
any location can be determined by adding vectorially the free-stream velocity 
to the contributions from all of the sources. Finally, the pressure distribution 
can be determined from Equation 2.74. 

This numerical procedure applied to a circular cylinder with a unit radius 
is illustrated in Figure 2.22. Here, only eight elements are shown. For this 
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Figure 2.22 Approximation to a circular cylinder with straight-line segments. 



Figure 2.23 Convergence of the numerical result to the exact solution for 
circular cylinder. 
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where 



The numerical calculation of the pressure distribution around a circular 
cylinder is compared in Figure 2.23 with the exact solution given by Equation 
2.78. As the number of segments increases, the approximate solution is seen 
to approach the exact solution rapidly. In Chapter Three this numerical 
method will be extended to include distributed vortices in addition to sources. 
In this way the lift of an arbitrary airfoil can be predicted. 


SUMMARY 

This chapter has introduced some fundamental concepts in fluid 
mechanics that will be expanded on and applied to explaining the aerody¬ 
namic behavior of airplane components in succeeding chapters. Potential flow 
methods will be used extensively with corrections given for Reynolds and 
Mach numbers. 


PROBLEMS 

2.1 Prove that the resultant static force on the face of a dam acts at the 
centroid of the dam’s area. 

2.2 Show that the incompressible Bernoulli’s equation (Equation 2.28) 
becomes p + pgh +ipV 2 = constant for a liquid, the weight of which is 
significant in comparison to the static pressure forces. (h is the depth of 
the streamline relative to an arbitrary horizontal reference plane.) 

2.3 A pilot is making an instrument approach into the University Park 
Airport, State College, Pennsylvania, for which the field elevation is 
listed at 378 m (1241ft) above sea level. The sea level barometric 
pressure is 763.3 mm Hg (30.05 in. Hg), but the pilot incorrectly sets the 
altimeter to 758.2 mm Hg (29.85 in. Hg). Will the pilot be flying too high 
or too low and by how much? [Note. Standard sea level pressure is equal 
to 760 mm Hg (29.92 in. Hg)]. 

2.4 Set to standard sea level pressure, an altimeter reads 2500 m (8200 ft). 
The outside air temperature (OAT) reads -15 °C (5°F). What is the 
pressure altitude? What is the density altitude? 

2.5 By integrating the pressure over a body’s surface, prove that the buoyant 
force on the body when immersed in a liquid is equal to the product of 
the volume of the displaced liquid, the liquid’s mass density, and the 
acceleration due to gravity. 


2.6 The hypothetical wake downstream of a two-dimensional shape is pic¬ 
tured below. This wake is far enough away from the body so that the 
static pressure through the wake is essentially constant and equal to the 
free-stream static pressure. Calculate the drag coefficient of the shape 
based on its projected frontal area. 



2.7 An incompressible flow has velocity components given by u = o>y and 
v = (ox, where w is a constant. Is such a flow physically possible? Can a 
velocity potential be defined? How is a> related to the vorticity? Sketch 
the streamlines. 

2.8 Derive Bernoulli’s equation directly by applying the momentum theorem 
to a differential control surface formed by the walls of a small stream- 
tube and two closely spaced parallel planes perpendicular to the velocity. 

2.9 A jet of air exits from a tank having an absolute pressure of 152,000 Pa 
(22 psi). The tank is at standard sea level (SSL) temperature. Calculate 
the jet velocity if it expands isentropically to SSL pressure. 

2.10 A light aircraft indicates an airspeed of 266 km/hr (165.2 mph) at a 
pressure altitude of 2400 m (7874 ft). If the outside air temperature is 
-10°C, what is the true airspeed? 

2.11 Prove that the velocity induced at the center of a ring vortex (like a smoke 
ring) of strength T and radius R is normal to the plane of the ring and has a 
magnitude of I72R. 

2.12 Write a computer program to solve the Biot-Savart equations numerically. 
This can be done by dividing a line vortex into finite, small straight-line 
elements. At a desired location the velocities induced by all of the elements 
can then be added vectorially to give the total resultant velocity. Check 
your program by using it to solve Problem 2.11. 
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THE GENERATION OF LIFT 


Lift is the component of the resultant aerodynamic forces on an airplane 
normal to the airplane’s velocity vector. Mostly, the lift is directed vertically 
upward and sustains the weight of the airplane. There are exceptions 
however. A jet fighter with a thrust-to-weight ratio close to unity in a steep 
climb may be generating very little lift with its weight being opposed mainly 
by the engine thrust. 

The component that is the major lift producer on an airplane and on 
which this chapter will concentrate is the wing. Depending on the airplane’s 
geometry, other components can contribute to or significantly affect the lift, 
including the fuselage, engine nacelles, and horizontal tail. These latter 
components will be considered, but to a lesser extent than the wing. 

WING GEOMETRY 

The top view, or planform, of a wing is shown in Figure 3.1. The length, 
b, from one wing tip to the other is defined as the wingspan. The chord, c, at 
some spanwise station, y, is the distance from the wing’s leading edge to its 
trailing edge measured parallel to the plane of symmetry in which the 
centerline chord, c 0 , lies. The chord generally varies with y so that, for 
purposes of characterizing wing geometry, a mean chord, c, is defined as the 
value that, when multiplied by the span, results in the planform area, S. 



The aspect ratio of a wing, A, is a measure of how long the span is with 
respect to the mean chord. Thus 

A=^ (3.2) 

For a rectangular planform where the chord is constant, this reduces to 

A=^- (3.3) 

c 
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As shown in Figure 3.1, a wing planform may be tapered and swept back. 
The taper ratio. A, is defined as the ratio of the tip chord, c„ to the midspan 
chord, Co. 



The sweep angle, A, is frequently measured relative to the quarter-chord 
line of the wing, that is, a line defined by the locus of points a quarter of the 
distance from the leading edge to the trailing edge. A, on occasion, is also 
measured relative to the leading edge. 

Usually the center portion of a wing is enclosed by the fuselage. In such 
an instance the wing’s aspect ratio and taper ratio are determined by ignoring 
the fuselage and extrapolating the planform shape into the centerline. The 
midspan chord in this instance is thus somewhat fictitious. The wing root is 
defined as the wing section at the juncture of the wing and fuselage. 
Occasionally, in the literature, one will find wing geometry characterized in 
terms of the wing root chord instead of the midspan chord. 

Approximately the aft 25 to 30% of a wing’s trailing edge is movable. On 
the outer one-third or so of the span the trailing edge on one side of the wing 
deflects opposite to that on the other. These oppositely moving surfaces are 
called ailerons; ailerons provide a rolling moment about the airplane’s long¬ 
itudinal axis. For example, when the aileron on the left wing moves down 
and the one on the right moves up, a moment is produced that tends to lift the 


left wing and lower the right one; this is a maneuver necessary in making a 
coordinated turn to the right. 

The inner movable portions of the wing’s trailing edge on both sides of 
the wing are kn own as the flaps. For takeoff and landing the flaps are lowered 
the same on both sides. There is no differential movement of the flaps on the 
left and right sides of the wing. The purpose of the flaps is to allow the wing 
to develop a higher lift coefficient than it would otherwise. Thus, for a given 
weight, the airplane can fly slower with the flaps down than with them up. 
Flaps, including leading edge flaps and the many different types of trailing 
edge flaps, will be discussed in more detail later. 

For some applications both ailerons are lowered to serve as an extension 
to the flaps. In such a case they are referred to as drooped ailerons, or 
flaperons. When flaperons are employed, additional roll control is usually 
provided by spoilers. These are panels that project into the flow near the 
trailing edge to cause separation with an attendant loss of lift. 

In order to understand and predict the aerodynamic behavior of a wing, it 
is expedient to consider first the behavior of two-dimensional airfoils. An 
airfoil can be thought of as a constant chord wing of infinite aspect ratio. 


AIRFOILS 

A considerable amount of experimental and analytical effort has been 
devoted to the development of airfoils. Much of this work was done by the 
National Advisory Committee for Aeronautics (NACA), the predecessor of 
the National Aeronautics and Space Administration (NASA). Reference 3.1 is 
an excellent summary of this effort prior to 1948. More recently NASA and 
others have shown a renewed interest in airfoil development, particularly for 
application to helicopter rotor blades, general aviation aircraft, and aircraft 
operating at transonic speeds of Mach 0.7 or higher. 

The development of an unflapped airfoil shape is illustrated in Figure 3.2. 
First, in Figure 3.2a, the chord line, c, is drawn. Next in Figure 3.2b, the 
camber line is plotted up from the chord a small distance z, which is a 
function of the distance from the leading edge. Next, as shown in Figure 3.2c, 
the semithickness is added to either side of the camber line. Also, the nose 
circle is centered on a tangent to the camber line at the leading edge and 
passes through the leading edge. Finally, an outer contour is faired around the 
skeleton to form the airfoil shape. Observe that the chord line is the line 
joining the ends of the mean camber line. 

The early NACA families of airfoils were described in this way, with the 
camber and thickness distributions given as algebraic functions of the 
chordwise position. However, for certain combinations of maximum thick- 
ness-to-chord ratios, maximum camber-to-chord ratios, and chordwise posi- 
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tion of maximum camber, tabulated ordinates for the upper and lower 
surfaces are available (Ref. 3.1). 

Before discussing the various families of airfoils in detail, we will 
generally consider the aerodynamic characteristics for airfoils, all of which 
can be influenced by airfoil geometry. 

To begin, an airfoil derives its lift from the pressure being higher on the 
lower surface of the airfoil than on the upper surface. If a subscript l denotes 
lower surface and “u” denotes upper surface, then the total lift (per unit 
span) on the airfoil will be 

L~l (pi-pjdx (3.5) 

The moment about the leading edge, defined positive nose up, will be 

M le = — j x(pi~p u )dx (3.6) 

In accord with Equation 2.12, the lift and moment can be expressed in 
terms of dimensionless coefficients. 


(3.7) 
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mLE (l/2)pVV 

Note that lowercase subscripts are used to denote coefficients for a 
two-dimensional airfoil, whereas uppercase subscripts are used for the three- 
dimensional wing. 

Writing 

Pi — Pu _ Pi ~Po _ Pu~Pe 

(l/2)pV 2 (1 m P V- (1 12)pv> 

and redefining x as the distance in chord lengths from the leading edge. 
Equations 3.5 and 3.6 become 


C, = f(C„- C P J dx 
us = ~j\(C n -C r ,)dx 


where the upper and lower pressure coefficients are defined according to 
Equation 2.73. 

The moment calculated from Equation 3.10 can be visualized as being 
produced by the resultant lift acting at a particular distance back from the 
leading edge. As a fraction of the chord, the distance x cp to this point, known 
as the center of pressure, can be calculated from 

-x cp C f - C mLE (3.11) 

Knowing x cp , the moment coefficient about any other point, x, along the 
airfoil can be written, referring to Figure 3.3, as 

Cn, =-(x cp -x)C ( (3.12) 


c, 
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It will be shown later that a point exists on an airfoil called the 
aerodynamic center about which the moment coefficient is constant and does 
not depend on C,. Denoting the location of the aerodynamic center by x ac . 
Equation 3.12 can be solved for the location of the center of pressure. 



Do not confuse the aerodynamic center with the center of pressure. 
Again, the aerodynamic center is the location about which the moment is 
constant, the center of pressure is the point at which the resultant lift acts. 

The progressive development of an airfoil shape is illustrated by 
.reference to Figure 3.4a and 3.4b. Historically airfoils developed ap¬ 
proximately in this manner. Consider first the simple shape of a thin, flat 
plate. 

Beginning with Figure 3.4a if the angle of attack of a thin, flat plate is 
suddenly increased from zero, the flow will appear for a moment as shown. 
Because of near-symmetry, there is practically no lift produced on the plate. 
However, because of viscosity, the flow at the trailing edge cannot continue 
to turn the sharp edge to flow upstream. Instead, it quickly adjusts to the 
pattern shown in Figure 3.4b. Here the flow leaves nearly tangent to the 
trailing edge. This condition is known as the Kutta condition after the 
German scientist, W. M. Kutta, who in 1902 first imposed the trailing edge 
condition in order to predict the lift of an airfoil theoretically. In Figure 3.4b 
observe that there is one streamline that divides the flow that passes over the 
plate from that below. Along this “dividing streamline,” the flow comes to 
rest at the stagnation point, where it joins perpendicular to the lower surface 
of the plate near the leading edge. As the flow progresses forward along this 
line, it is unable to adhere to the surface around the sharp leading edge and 
separates from the plate. However, it is turned backward by the main flow 
and reattaches to the upper surface a short distance from the leading edge. 
The resulting nonsymmetrical flow pattern causes the fluid particles to ac¬ 
celerate over the upper surface and decelerate over the lower surface. Hence, 
from Bernoulli’s equation, there is a decrease in air pressure above the plate 
and an increase below it. This pressure difference acting on the airfoil 
produces a lift. 

If the angle of attack of the plate is too great, the separated flow at the 
leading edge will not reattach to the upper surface, as shown in Figure 3.4c. 
When this occurs, the large separated region of unordered flow on the upper 
surface produces an increase in pressure on that surface and hence a loss in 
lift. This behavior of the airfoil is known as stall. Thus the limit in Q, that is, 
C /mM , is the result of flow separation on the upper surface of the airfoil. 

To improve this condition, one can curve the leading edge portion of the 
flat plate, as shown in Figure 3.4d, to be more nearly aligned with the flow in 



Figure 3.4 Progressive development of airfoil shapes, (a) Flat plate at sudden 
angle of attack—no lift, (b) Flat plate at angle of attack in steady flow and generating 
lift, (c) Flat plate experiencing leading edge separation and loss of lift (stall). ( d ) Flat 
plate with curved leading edge to prevent leading edge separation, (e) Airfoil with 
thickness and camber to delay stall, (f) Airfoil with trailing edge separation. 
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that region. Such a shape is similar to that used by the Wright Brothers. This 
solution to the separation problem, as one might expect, is sensitive to angle 
of attack and only holds near a particular design angle. However, by adding 
thickness to the thin, cambered plate and providing a rounded leading edge, 
the performance of the airfoil is improved over a range of angles, with the 
leading edge separation being avoided altogether. Thus, in a qualitative sense, 
we have defined a typical airfoil shape. Camber and thickness are not needed 
to produce lift (lift can be produced with a flat plate) but, instead, to increase 
the maximum lift that a given wing area can deliver. 

Even a cambered airfoil of finite thickness has its limitations, as shown in 
Figure 3.4/. As the angle of attack is increased, the flow can separate initially 
near the trailing edge, with the separation point progressively moving forward 
as the angle of attack continues to increase. 

The degree to which the flow separates from the leading or trailing edge 
depends on the Reynolds number and the airfoil geometry. Thicker airfoils 
with more rounded leading edges tend to delay leading edge separation. This 
separation also improves with increasing values of the Reynolds number. 

Leading edge separation results in flow separation over the entire airfoil 
and a sudden loss in lift. On the other hand, trailing edge separation is 
progressive with angle of attack and results in a more gradual stalling. The 
situation is illustrated in Figures 3.5 and 3.6 (taken from Ref. 3.1). In Figure 3.5 
note the sharp drop in C, at an a of 12° for R = 3 x 10 6 , whereas for 
R = 9x 10 6 , the lift curve is more rounded, with a gradual decrease in Q 
beyond an a of 14°. In Figure 3.6, for a thicker airfoil with the same camber, 
the lift increases up to an angle of approximately 16° for all R values tested. 
At this higher angle, even for R = 9 x 10*, it appears that leading edge 
separation occurs because of the sharp drop in Q for a values greater than 
16°. From a flying qualities standpoint, an airfoil with a well-rounded lift curve 
is desirable in order to avoid a sudden loss in lift as a pilot slows down the 
airplane. However, other factors such as drag and Mach number effects must 
also be considered in selecting an airfoil. Hence, as is true with most design 
decisions, the aerodynamicist chooses an airfoil that represents the best 
compromise to conflicting requirements, including nonaerodynamic con¬ 
siderations such as structural efficiency. 

Figures 3.5 and 3.6 illustrate other characteristics of airfoil behavior that 
will be considered in more detail later. Observe that the lift curve, Q versus 
a, is nearly linear over a range of angles of attack. Notice also that the slope, 
dCJda, of the lift curve over the linear portion is unchanged by deflecting the 
split flap. The effect of lowering the flap or, generally, of increasing camber is 
to increase Q by a constant increment for each a in the linear range. Thus the 
angle of attack for zero lift, a 0 (, is negative for a cambered airfoil. In the case 
of the 1408 airfoil pictured in Figure 3.5, a 0 / equals —12.5°, with the split flap 
deflected 60°. 


Figure 3.5 Characteristics of the NACA 1408 airfoil. 
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If a is increased beyond the stall G will again begin to increase before 
dropping off to zero at an a of approximately 90°. The second peak in Q is 
generally not as high as that which occurs just before the airfoil stalls. S. P. 
Langley, in his early experiments, noted these two peaks in the Q versus a 
curve but chose to fair a smooth curve through them. Later, the Wright 
Brothers observed the same characteristics and were troubled by Langley’s 
smooth curve. After searching Langley’s original data and finding that he, too, 
had a “bump” in the data, Wilbur Wright wrote to Octave Chanute on 
December 1, 1901. 

“// he ( Langley ) had followed his observations, his line would probably have 
been nearer the truth. I have myself sometimes found it difficult to let the lines run 
where they will, instead of running them where I think they ought to go. My 
conclusion is that it is safest to follow the observations exactly and let others do 
their own correcting if they wish ” (Ref. 1.1). 

To paraphrase the immortal Wilbur Wright, “Do not ‘fudge’ your data—it may 
be right.” 


AIRFOIL FAMILIES 
NACA Four-Digit Series 

Around 1932, NACA tested a series of airfoil shapes known as the 
four-digit sections. The camber and thickness distributions for these sections 
are given by equations to be found in Reference 3.1. These distributions were 
not selected on any theoretical basis, but were formulated to approximate 
efficient wing sections in use at that time, such as the well-known Clark-Y 
section. 

The four-digit airfoil geometry is defined, as the name implies, by four 
digits; the first gives the maximum camber in percent of chord, the second the 
location of the maximum camber in tenths of chord, and the last two the 
maximum thickness in percent of chord. For example, the 2412 airfoil is a 
12% thick airfoil having a 2% camber located 0.4c from the leading edge. The 
2412 airfoil is pictured in Figure 3.7 along with other airfoils yet to be 
described. 

NACA Five -Digit Series 

The NACA five-digit series developed around 1935 uses the same thick¬ 
ness distribution as the four-digit series. The mean camber line is defined 
differently, however, in order to move the position of maximum camber 
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Whitcomb-type supercritical airfoil 


NASA GA(W)-1 

Figure 3.7 Comparison of various airfoil shapes. 

forward in an effort to increase C (inu . Indeed, for comparable thicknesses and 
cambers, the Q mtx values for the five-digit series are 0.1 to 0.2 higher than 
those for the four-digit airfoils. The numbering system for the five-digit series 
is not as straightforward as for the four-digit series. The first digit multiplied 
by 3/2 gives the design lift coefficient of the airfoil. The next two digits are 
twice the position of maximum camber in percent of chord. The last two 
digits give the percent thickness. For example, the 23012 airfoil is a 12% thick 
airfoil having a design C; of 0.3 and a maximum camber located 15% of c 
back from the leading edge. This airfoil is also pictured in Figure 3.7. 
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NACA l-Series (Series 16) 

The NACA 1-series of wing sections developed around 1939 was the first 
series based on theoretical considerations. The most commonly used 1-series 
airfoils have the minimum pressure located at the 0.6c point and are referred 
to as series-16 airfoils. The camber line for these airfoils is designed to 
produce a uniform chordwise pressure difference across it. In the thin airfoil 
theory to follow, this corresponds to a constant chordwise distribution of 
vorticity. 

Operated at its design Q, the series-16 airfoil produces its lift while 
avoiding low-pressure peaks corresponding to regions of high local velocities. 
Thus the airfoil has been applied extensively to both marine and aircraft 
propellers. In the former application, low-pressure regions are undesirable 
from the standpoint of cavitation (the formation of vaporous cavities in a 
flowing liquid). In the latter, the use of series-16 airfoils delays the onset of 
deleterious effects resulting from shock waves being formed locally in regions 
of high velocities. 

Series-1 airfoils are also identified by five digits as, for example, the 
NACA 16-212 section. The first digit designates the series; the second digit 
designates the location of the minimum pressure in tenths of chord. Following 
the dash, the first number gives the design Q in tenths. As for the other 
airfoils, the last two digits designate the maximum thickness in percent of 
chord. The 16-212 airfoil is shown in Figure 3.7. 

NACA 6Series 

The 6-series airfoils were designed to achieve desirable drag, com¬ 
pressibility, and C| mw( performance. These requirements are somewhat 
conflicting, and it appears that the motivation for these airfoils was primarily 
the achievement of low drag. The chordwise pressure distribution resulting 
from the combination of thickness and camber is conducive to maintaining 
extensive laminar flow over the leading portion of the airfoil over a limited 
range of Q values. Outside of this range, C d and C, mM values are not too much 
different from other airfoils. 

The mean lines used with the 6-series airfoils have a uniform loading 
back to a distance of x/c = 2. Aft of this location the load decreases linearly. 
The a = 1 mean line corresponds to the uniform loading for the series-16 
airfoils. 

There are many perturbations on the numbering system for the 6-series 
airfoils. The later series is identified, for example, as 

NACA 65i-212 a = 0.6 
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Figure 3.8 Aerodynamic characteristics of the NACA 65, -212, 
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Here 6 denotes the series; the numeral 5 is the location of the minimum 
pressure in tenths of chord for the basic thickness and distribution; and the 
subscript 1 indicates that low drag is maintained at Q values of 0.1 above and 
below the design Q of the 0.2, denoted by the 2 following the dash. Again, the 
last two digits specify the percentage thickness. If the fraction, a, is not 
specified, it is understood to equal unity. The 65 r 212 airfoil is shown in Figure 
3.7. 

Lift and drag curves for the 65,-212 airfoil are presented in Figure 3.8. 
Notice the unusual shape of C d versus C,, where the drag is significantly lower 
between C, values of approximately 0 to 0.3. In this region, for very smooth 
surfaces and for Reynolds numbers less than 9 x 10 s , extensive laminar flow is 
maintained over the surface of the foil with an attendent decrease in the skin 
friction drag. This region, for obvious reasons, is known as the “drag bucket.” 
In practice this laminar flow, and resulting low drag, is difficult to achieve 
because of contamination by bugs or by structurally transmitted vibration that 
perturbs the laminar boundary layer, causing transition. Chapter Four will 
discuss the drag of these airfoils in more detail. 


MODERN AIRFOIL DEVELOPMENTS 

Systematic series of airfoils have given way, at least in part, to speci¬ 
alized airfoils designed to satisfy particular requirements. These airfoils are 
synthesized with the use of sophisticated computer programs such as the one 
described in Reference 3.5, which will be discussed in more detail later. One 
such special purpose airfoil is the so-called supercritical airfoil reported on in 
References 3.6 and 3.7. This airfoil has a well-rounded leading edge and is 
relatively flat on top with a drooped trailing edge. For a constant thickness of 
12%, wind tunnel studies indicate a possible increase of approximately 15% in 
the drag-divergence Mach number for a supercritical airfoil as compared to a 
more conventional 6-series airfoil. In addition, the well-rounded leading edge 
provides an improvement in at low speeds over the 6-series, which has 
sharper leading edges. 

A qualitative explanation for the superior performance of the super¬ 
critical airfoil is found by reference to Figure 3.9. At a free-stream Mach 
number as low as 0.7 or so depending on the shape and C,, a conventional 
airfoil will accelerate the flow to velocities that are locally supersonic over the 
forward or middle portion of its upper surface. The flow then decelerates 
rapidly through a relatively strong shock wave to subsonic conditions. This 
compression wave, with its steep positive pressure gradient, causes the 
boundary layer to thicken and, depending on the strength of the shock, to 
separate. This, in turn, causes a significant increase in the drag. The mi nim um 
value of the free-stream Mach number for which the local flow becomes 



supersonic is referred to as the critical Mach number. As this value is 
exceeded by a few hundredths, the shock wave strengthens sufficiently to 
cause the drag to rise suddenly. This free-stream Mach number is known as 
the drag-divergence Mach number. 

The supercritical airfoil also accelerates the flow to locally supersonic 
conditions at free-stream Mach numbers comparable to the 1- or 6-series 
airfoils. However, the supercritical airfoil is shaped, so that around its design 
lift coefficient, the flow decelerates to subsonic conditions through a dis¬ 
tribution of weak compression waves instead of one strong one. In this way 
the drag-divergence Mach number is increased substantially. 
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Although the possibility of such airfoils was known for some time, their 
successful development in modern times is attributed to R. T. Whitcomb. A 
Whitcomb-type supercritical airfoil is pictured in Figure 3.7. 

Tested at low speeds, the supercritical airfoils were found to have good 
O nu values as well as low Q values at moderate lift coefficients. As a result, 
another family of airfoils evolved from the supercritical airfoils, but for 
low-speed applications. These are the “general aviation” airfoils, designated 
GA(W) for general aviation (Whitcomb). The GA(W)-1 airfoil is the last of the 



Figure 3.10a Effect of Reynolds number on section characteristics of the 
GA(W)-1 airfoil Model smooth, M = 0.15. 
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(b) 

Bgoro 3.11 Comparison of section characteristics of NASA GA(W)-1 airfoil and 
NACA 652-415 and 653-418 airfoils. M = 0.20; fl«6xl0 6 . (a) Variation of C, and C m 
with a. (6) Variation of C a with C,. 
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airfoils pictured in Figure 3.7. Test results for this airfoil are reported in 
Reference 3.8, where its C, mu values are shown to be about 30% higher than 
those for the older NACA 65-series airfoils. In addition, above Q values of 
around 0.6, its drag is lower than the older laminar flow series with standard 
roughness. These data are presented in Figure 3.10 for the GA(W)-1 airfoil. 
Comparisons of CJ mix and C d for this airfoil with similar coefficients for other 
airfoils are presented in Figures 3.11 and 3.12. 

Observe that the performance of the GA(W)-1 airfoil is very Reynolds 
number-dependent, particularly C; mai , which increases rapidly with Reynolds 
number from 2 to 6 million. At the time of this writing, the GA(W) airfoil is 
beginning to be employed on production aircraft. The same is true of the 
supercritical airfoil. Indeed, the supercritical airfoil is being used on both the 
Boeing YC-14 and McDonnell-Douglas YC-15 prototypes currently being tested 
for the advanced medium STOL transport (AMST) competition. At the time of 
this writing, NASA is adopting a new nomenclature for the GA(W) airfoils. They 
will be designated by LS (low speed) or MS (medium speed) followed by four 
digits. For example, the GA(W)-1 airfoil becomes LS(1)-0417. The (1) designates 
a family. The 04 refers to a design lift coefficient of 0.4, and 17 is the maximum 
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Fjgure 3.12 Comparison of maximum lift coefficient of the GA(W)-1 airfoil with 
other NACA airfoils M = 0.15. 
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thickness in percent of chord. For more information on modern airfoils, consult 
Reference 3.36. 

PREDICTION OF AIRFOIL BEHAVIOR 

In Chapter Two it was noted that the concepts of a point vortex and a 
point source could be extended to a continuous distribution of the elementary 
flow functions. In that chapter a distribution of sources in a uniform flow was 
found to produce a nonlifting body of finite thickness. In the case of the 
circular cylinder, the addition of a vortex also produced lift. 

Comparable to the continuous distribution of sources pictured in Figure 
2.20, consider a similar distribution of vortices as illustrated in Figure 3.13. 
Such a distribution is referred to as a vortex sheet. If y is the strength per unit 
length of the sheet, y Ax will be the total strength enclosed by the dashed 
contour shown in the figure. The contour is taken to lie just above and below 
the sheet. Ax is sufficiently small so that the velocity tangent to the sheet, v, 
can be assumed to be constant. Because of the symmetry to the flow provided 
by any one segment of the sheet, the tangential velocity just below the sheet 
is equal in magnitude but opposite in direction to that just above the sheet. 
From Equation 2.55, relating circulation to the strength of a vortex, it follows 
that 

y Ax = 2v Ax 
or 



Note the similarity of this relationship to that expressed by Equation 2.82. 
However, in the case of Equation 3.14, the velocity is tangent to the vortex 
sheet whereas, for Equation 2.82, the velocity is normal to the line on which 
the sources lie. 

Consider now the thin airfoil pictured in Figure 3.14. If the airfoil is 
producing a lift, the pressure on the lower surface is greater than that on the 
upper. Hence, from Bernoulli’s equation, the velocity on the upper surface is 
greater than the velocity on the lower surface. Letting this difference in 



L —- t Ax = total strength of 


Figure 3.13 Distributed vortices in a two-dimensional flow (vortex sheet). 
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Figure 3.14 The velocity difference across a lifting thin airfoil. 


velocity across the airfoil equal 2v, the upper and lower velocities can be 
written as 

V u = V+v 

and 

V,= V-v 

Thus the flow field around the airfoil is the same as that which would be 
produced by placing, in a uniform flow of velocity V, a vortex sheet of unit 
strength 2v along the airfoil. 

The contribution to the lift of a differential length of the airfoil will be 
dl = (p, ~p u ) dx 

Or, using Bernoulli’s equation, this becomes, 
dl = pV(2v) dx 

Since 2v is the unit vortex strength, the Kutta-Joukowski law (Equation 2.81) 
is found to hold for the airfoil element. 

dl = p Vy dx 

or, integrating Equation 3.18 over the entire chord, 

l = pVT (3.15) 

where T is the total circulation around the airfoil given by 

T =f'ydx (3.16) 

In order to predict the lift and moment on the airfoil, one must find the 
chordwise distribution of y(x) that will produce a resultant flow everywhere 
tangent to the mean camber line (thin airfoil approximation). In addition, the 
Kutta condition is applied at the trailing edge to assure that the flow leaves 
the trailing edge tangent to the mean camber line at that point. This is a 
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necessary condition; otherwise, the resulting flow will appear similar to Figure 
3.4 a with the lift being equal to zero. 

An analytical solution to the thin airfoil will be obtained later but, first, let 
us consider a numerical approach to predicting the lift and moment of an 
airfoil. 

As a gross approximation to the distributed vorticity along the airfoil, the 
distribution will be replaced by only one vortex of unknown strength, T. 
However, T will be placed at a particular point on the airfoil, at the quarter- 
chord point. The boundary condition and the Kutta condition will be satisfied 
at only one point, the three-quarter-chord point. This approximation, known 
as Weissinger’s approximation, is illustrated in Figure 3.15 for a flat-plate 
airfoil. 

The velocity induced at 3c/4 by f placed at c/4 will be 



Assuming a to be a small angle, it follows that 

= Va (3.17) 

or 

r = TTcVa (3.18) 

From the Kutta-Joukowski relationship, L = pVT, so that 

L = pircV 2 a (3.19) 

Expressing lift in terms of the lift coefficient and using Equation 3.19 leads to 
C, = 2 t ra (3.20) 

where a is the angle of attack in radians. 

The expression agrees identically with the theoretical solution of this 
problem that follows. Notice that the result predicts the slope of the lift curve, 
dCJda, to be 2Wrad. Experimentally this figure is usually found to be 



Figure 3.15 Weissinger's approximation to a thin airfoil. 


somewhat less. Figures 3.5 and 3.6, for example, show a value of around 
0.105/deg or 6.02/rad. 

The approximation of Figure 3.15 can be improved on by dividing the 
airfoil chord into a number of equal segments and placing a vortex of 
unknown strength at the quarter-chord point of each segment. The unknown 
strengths are determined by assuring that the normal velocity vanishes at the 
three-quarter-chord point of each segment. With the last control point down¬ 
stream of the last vortex singularity, the Kutta condition is assured. 

To illustrate this numerical solution of the thin airfoil, consider Figure 
3.16. Here, a circular arc airfoil having a unit chord length with a maximum 
camber ratio of z is operating at an angle of attack a. 

If it is assumed that 

2 «1 (3.21) 

the radius of curvature, R, of the airfoil will be related to z approximately by 



The slope of the camber line relative to the chord line (the angle <f> in Figure 
3.16) at any distance x can be determined from the geometry of the figure. 

(f> - (4- 8x)z (3.22) 



Figure 3.16 A circular arc airfoil approximated by two vortices. 
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The component of V normal to the mean camber line and directed 
upward is thus 

»„= V(a-4>) 

It follows that at control points 1 and 2 located at x values of 3c/8 and 7c/8, 
respectively, the two vortices similating the airfoil must induce velocities 
downward given by 

Vi = V(a — z) at 1 (3.23a) 

u, = V(a+3z) at 2 (3.23fc) 

The problem is linearized by Equation 3.21 so that the vortices, yi and y 2 , 
are taken to lie on the chord line. Thus, according to Equation 2.56, the total 
velocities induced at the two control points by the two vortices will be 

v, =~<2y,-2yi) all - (3.24a) 

= ^;(| Ti 4- Z-Yi) at 2 (3.246) 

Equating Equations 3.23 and 3.24 results in 

= + z) (3.23a) 

y2 = ~{«+Sz) (3.256) 

Applying the Kutta-Joukowski law to each vortex results not only in a 
predicted total lift, but also in a moment. In coefficient form the lift and 
moment (about the leading edge) become 

C, = 27r(a + 2z) (3.26a) 

C«= -»(f+£) (3.266) 

The moment coefficient about the leading edge can be transferred to the 
quarter-chord point by using 



Thus, 

c m = < 3 - 27 > 

1 '* 4 

This simple, two-point model results in several important observations 
that are in agreement with more exact solutions. First, note that Equation 3.26a 
shows the lift coefficient to be a linear combination of a and z. Thus, 
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cambering an airfoil will not change the slope of the lift curve. Second, it is 
predicted that the moment about the quarter chord will be independent of a. 
Hence, this point is predicted to be the aerodynamic center. 

As one divides the airfoil into a greater and greater number of elements, 
the resulting y distributions will approach the theoretical pressure distribution 
predicted on the basis of continuous y distributions. The strength, y, of a 
vortex placed at the c/4 point of an element of length Ax will be related to the 
pressure jump, A p, across the element by 

Ap=^2 (3.28) 

Figure 3.17 presents a comparison, for the flat-plate airfoil, between the 
pressure distribution obtained using the foregoing numerical procedure with 
that based on a continuous distribution of y along the chord. It is seen that 
the numerical results rapidly converge to the continuous solution as the 
number of elements increases. In preparing this figure it should be noted that 
Ap, given by Equation 3.28, has been expressed in coefficient form and plotted 
at the location of each point vortex. Figure 3.18 presents a similar comparison 
for the circular arc airfoil. In this case a is taken to be zero, avoiding the 
infinitely negative C p at the leading edge. 



Figure 3.17 Comparison of numerical calculation of chordwise lift distribution 
with analytical prediction for a flat-plate airfoil at 10° angle of attack. 
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Figure 3.18 Comparison of numerical calculation of chordwise lift distribution 
with analytical prediction for a 4% cambered, circular arc airfoil. 


The numerical model predicts the lift in exact agreement with more 
precise analytical models. However, the moment coefficient, given by Equa¬ 
tion 3.27, is only three-quarters of that obtained by analytical means. Figure 
3.19 shows that the exact value is approached rapidly, however, as the 
number of segments increases. 

As indicated by Figure 3.19, the exact value of the moment coefficient 
about the aerodynamic center (c/4) for the circular arc airfoil is given by 

C m>c = -7T2 (3.29) 

Using Equation 3.13, the location of the center of pressure can be found 
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Figure 3.19 Numerical calculation of moment coefficient compared with analy¬ 
tical prediction for 4% cambered, circular arc airfoil. 


Observe that as Q decreases, the center of pressure moves aft, ap¬ 
proaching infinity as Q goes to zero. This movement of the center of pressure 
is opposite to what was believed to be true by the early pioneers in aviation. 
The Wright Brothers were probably the first to recognize the true nature of 
the center-of-pressure movement as a result of their meticulous wind tunnel 
tests. 

Analytical solutions to the thin airfoil can be found in several texts (e.g.. 
Ref. 3.2 and 3.3). Here, the airfoil is replaced by a continuous distribution of 
vortices instead of discrete point vortices, as used with the numerical solu¬ 
tion. 

Referring to Figure 3.20, without any loss of generality, the airfoil is 
taken to have a unit chord lying along the x-axis with the origin at the leading 
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edge. The shape of the camber line is given by z(x ), and it is assumed that 


With this assumption the problem is linearized and made tractable by replac¬ 
ing the airfoil with a vortex sheet of unit strength y(x) lying along the chord 
line instead of along the camber line. 

At the point x 0 , the downward velocity induced by an elemental vortex of 
strength y(x) dx located at x, according to Equation 2.56, will be given by 


or, integrating over the chord. 


In order to satisfy the boundary condition that the flow be tangent 
everywhere to the mean camber line, it follows that, to a small angle 
approximation. 


w(xq) / dz\ 

V “ \dx},„ 


Thus, given a and z(x), the following integral equation must be solved 
for y(x). 

Sv/:^r-«-(£)„ 

In addition, y(x) must vanish at the trailing edge in order to satisfy the 
Kutta condition. Otherwise, the induced velocity will be infinite just down¬ 
stream of this point. 

Equation 3.33 is solved by first transforming to polar coordinates. 

Letting 


Equation 3.33 becomes 


y(0) de 
cos 6 — cos 6 


On the basis of the more sophisticated method of conformal mapping 
(e.g., see Ref. 3.4), it is known that y(x) is generally singular at the leading 
edge approaching infinity as l/x. Thus we will assume a priori that Equation 
3.35 can be satisfied by a y(0) distribution of the form 


(3.36) 


PREDICTION OF AIRFOIL BEHAVIOR 91 


Using the relationships 


^ ^cos (n - 1)0 — cos (n + 1)0 j = sin nO sin 0 

f 7 * cos nOdd sin n0 o 
Jo cos 6 - cos 0 O sin 0 O 


Equation 3.35 becomes 


- 2 A„ cos nO = a — -p 


Multiplying both sides of the preceding equation by cos m0 (m =0, 1, 2,..., 
n ,...) and integrating from 0 to tt leads to 


a 1 f w dz M 

Ao = «-*l Tx dd 

A„=— f ^-cos n6d$ 
tt Jo dx 


Thus, knowing the shape of the mean camber line, the coefficients A 0 , A u 
A 2 ,... can be determined either in closed form or by graphical or numerical 
means (see Ref. 3.1). Having these coefficients, Q and C m can then be easily 
determined from the Kutta-Joukowski relationship. 

The lift and moment about the leading edge are given by 


- = pVy(x) dx 
f.= J pVy(x)xdx 


From these and using Equation 3.36, 


C mLE =-f(Ac + A,-f) (3.40) 

It follows that C m about the quarter-chord point is independent of a, so that 
this point is the aerodynamic center, with the moment coefficient being given 
by 

Since a is contained only in the A 0 coefficient, it can be concluded 
immediately without considering the actual form of z(x) that Q is given by a 
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linear combination of a and a function of z. Thus, camber changes can be 
expected to affect the angle of zero lift but not the slope of the lift curve. 

Reference to airfoil data, such as that presented in Figures 3.5 and 3.6, 
will show that the predictions of thin airfoil theory are essentially correct. 
There is a range of angles of attack over which the lift coefficient varies 
linearly with a. The slope of this lift curve is usually not as high as the theory 
predicts, being approximately 4 to 8% less than the theoretical value. For 
many purposes an assumed value of 0.1 G/deg is sufficiently accurate and is a 
useful number to remember. Experimental data also show the aerodynamic 
center to be close to the quarter-chord point. The effects of camber on G> 
dQlda, and Gn ac are also predicted well. 

Recently large numerical programs have been developed to predict the 
performance of airfoils that incorporate Reynolds numbdr and Mach number 
effects. These are typified by Reference 3.5, which will be described briefly. 
This program begins by calculating the potential flow around the airfoil. In 
order to allow for both finite thickness and circulation, the airfoil contour is 
approximated by a closed polygon, as shown in Figure 3.21. A continuous 
distribution of vortices is then placed on each side of the polygon, with the 
vortex strength per unit length, y, varying linearly from one corner to the next 
and continuous across the comer. Figure 3.22 illustrates this model for two 
sides connecting corners 3, 4, and 5. Control points are chosen midway 
between the comers. The values of the vortex unit strengths at the corners are 
then found that will induce velocities at each control point tangent to the 
polygon side at that point. Note, however, that if there are n corners and 
hence n + 1 unknown y values at the corners, the n control points provide 
one less equation than unknowns. This situation is remedied by applying the 
Kutta condition at the trailing edge. This requires that y n+ i = —y it assuring 
that the velocities induced at the trailing edge are finite. 

Having determined the vortex strengths, the velocity field and, hence, the 



Figure 3.21 Approximation of airfoil contour by closed polygon. 




Figure 3.22 Vortex distributions representing airfoil contour. 

pressure distribution around the airfoil can be calculated. This result is then 
used to calculate the boundary layer development over the airfoil, including 
the growth of the laminar layer, transition, the growth of the turbulent layer, 
and possible boundary layer separation. The airfoil shape is then enlarged 
slightly to allow for the boundary layer thickness and the potential flow 
solutions are repeated. The details of this iterative procedure are beyond the 
scope of this text. 


MAXIMUM LIFT 


Airfoil theory based on potential flow methods predicts the lift of an 
airfoil in its linear range but does not provide any information concerning 
maximum lift capability. As discussed previously, G ro „, is determined by flow 
separation, which is a “real fluid” effect. Separation is difficult to predict 
analytically, so the following material on G raaa is mainly empirical. 

Typically, conventional airfoils without any special high-lift devices will 
deliver G nu values of approximately 1.3 to 1.7, depending on Reynolds 
number, camber, and thickness distribution. The appreciable dependence of 
C !aax on R shown in Figure 3.12 for the GA(W)-1 airfoil is typical of other 
airfoils. Figure 3.23 presents G m „ as a function of J? and thickness ratio for 
NACA four-digit airfoils having a maximum camber of 2%, located 40% of 
the chord back from the leading edge. At intermediate thickness ratios of 
around 0.12, the variation of with R parallels that of the 17% thick 
GA(W)-1 airfoil. Note, at least for this camber function, that a thickness ratio 
of 12% is about optimum. This figure is taken from Reference 3.14. This same 
reference presents the following empirical formula for G raax f° r NACA four¬ 
digit airfoils at an R of 8 x 10 6 . 


= 1-67 + 7.8 pz 


« , (0. 123 + 0.022 - 0.5z - tf 

2.6 pii 


(3.42) 


1 
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Figure 3.23 Variation of C, wtx with thickness ratio of NACA 24xx airfoils for 
various Reynolds numbers. (B. W. McCormick, Aerodynamics of V/STOL Flight , 
Academic Press, Inc, 1967. Reprinted by permission of Academic Press, Inc.) 


t, z, and p are thickness, maximum camber, and position of maximum camber, 
respectively, expressed as a fraction of the chord. For example, for a 2415 
airfoil, 

t =0.15 
2 = 0.02 
p =0.40 

so that according to Equation 3.42, 

CL, = 1.70 

For a plain wing (unffapped), there is little effect of aspect ratio or taper 
ratio on C w Even the presence of a fuselage does not seem to have much 
effect. As the angle of attack of a wing increases, C L ^ is reached and any 
further increase in a will result in a loss of lift. Beyond the wing is said 
to be stalled. Although taper ratio does not significantly affect the overall wing 
Ctw ( an d wing twist) significantly affect what portion of the wing stalls 
first. As the taper ratio is decreased, the spanwise position of initial stall moves 
progressively outboard. This tendency is undesirable and can be compensated 
for by “washing out” (negative twist) the tips. One usually wants a wing to 
stall inboard initially for two reasons. First, with inboard stall, the turbulence 
shed from the stalled region can shake the tail, providing a built-in stall 
warning device. Second, the outboard region being unstalled will still provide 


aileron roll control even though the wing has begun to stall. The lift charac¬ 
teristics of three-dimensional wings will be treated in more detail later. 

Flaps 

An examination of all of the airfoil data presented in Reference 3.1 
discloses that the greatest value of Q mx one can expect at a high Reynolds 
number from an ordinary airfoil is around 1.8. This maximum value is 
achieved by the NACA 23012 airfoil. Another 12% thick airfoil, the 2412, 
delivers the second highest value, 1.7. 

In order to achieve higher CL, values for takeoff and landing without 
unduly penalizing an airplane’s cruising performance, one resorts to the use of 
mechanical devices to alter temporarily the geometry of the airfoil. These 
devices, known as flaps, exist in many different configurations, the most 
common of which are illustrated in Figure 3.24. In addition to the purely 
mechanical flaps, this figure depicts flaps that can be formed by sheets of air 
exiting at the trailing edge. These “jet flaps” can produce CL, values in 
excess of those from mechanical flaps, provided sufficient energy and 
momentum are contained in the jet. Frequently one uses the tqrms “powered” 
and “unpowered” to distinguish between jet and mechanical flaps. 

The effect of a mechanical flap can be seen by referring once again to 
Figure 3.6. Deflecting the flap, in this case a split flap, is seen to shift the lift 
curve upward without changing the slope. This is as one might expect from 
Equation 3.26, since deflecting the flap is comparable to adding camber to the 
airfoil. 

Some flap configurations appear to be significantly better than others 
simply because, when deflected, they extend the original chord on which the 
lift coefficient is based. One can determine if the flap is extensible, such as the 
Fowler or Zap flaps in Figure 3.24, by noting whether or not the slope of the 
lift curve is increased with the flap deflected. Consider a flap that, when 
deflected, extends the original chord by the fraction x. The physical lift curve 
would have a slope given by 


^=|pV 2 (l +x)cC,„ (3.43) 

since (l+x)c is the actual chord. C la does not depend significantly on 
thickness or camber; hence, the lift curve slope of the flapped airfoil based on 
the unflapped chord, c, would be 

Q o (flapped) = (1 + x)Q a (unflapped) 

Now the maximum lift, expressed in terms of the extended chord and Ci mil , 
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Figure 3.24 Flap configurations. 


(based on that chord) would be, 

Thus Cf^ based on the original chord becomes 

CLx = U + x)Ci waXe 
or 

G^-CUAl+x) 

_ r Cl (unflapped) 

_ Q a (flapped) 


(3.44) 



Figure 3.25 Performance of plain flaps, (a) Variation of maximum section lift 
coefficient with flap deflection for several airfoil sections equipped with plain 
flaps. ( b ) Variation of optimum increment of maximum section lift coefficient 
with flap chord ratio for several airfoil sections equipped with plain flaps, (c) 
Effect of gap seal on maximum lift coefficient of a rectangular Clark-Y wing 
equipped with a full-span 0.20c plain flap. A = 6, Ft = 0.6 x 10 4 . 







Figure 3.26 Variation of maximum section lift coefficient with flap deflection for 
three NACA 230-series airfoils equipped with split flaps of various sizes. R = 
3.5 x 10®. (a) NACA 23012 airfoil section, (fa) NACA 23021 airfoil section, (c) NACA 
23030 airfoil section. 


Figures 3.25 to 3.30 and Tables 3.1 and 3.2 present section data on plain, 
split, and slotted flaps as reported in Reference 3.15. With these data one 
should be able to estimate reasonably accurately the characteristics of an 
airfoil section equipped with flaps. 

A study of this data suggests the following: 

Plain Flaps 

1. The optimum flap chord ratio is approximately 0.25. 

2. The optimum flap angle is approximately 60°. 

3. Leakage through gap at flap nose can decrease C (niM by approximately 0.4. 

4. The maximum achievable increment in is approximately 0.9. 

Split Flaps 

1. The optimum flap chord ratio is approximately 0.3 for 12% thick airfoils, 
increasing to 0.4 or higher for thicker airfoils. 

2. The optimum flap angle is approximately 70°. 

3. The maximum achievable increment in Q miX is approximately 0.9. 

4. increases nearly linearly with log i? for 0.7 x 10 6 < R < 6 x 10 6 . 

5. The optimum thickness ratio is approximately 18%. 
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Figure 3.27 Variation of maximum section lift coefficient with Reynolds number 
for several NACA airfoil sections with and without 0.20c split flaps deflected 60°. 
(a) Smooth airfoil, (b) Airfoil with leading edge roughness. 


Slotted Flaps 

1. The optimum flap chord ratio is approximately 0.3. 

2. The optimum flap angle is approximately 40° for single slots and 70° for 
double-slotted flaps. 

3. The optimum thickness ratio is approximately 16%. 

4. is sensitive to flap (and vane) position. 

5. The maximum achievable increment in is approximately 1.5 for single 
slots and 1.9 for double slotted flaps. 



Figure 3.28 ' Contours of flap and vane positions for maximum section lift coefficient for several airfoil sections equipped 
with double-slotted flaps, (a) NACA 23012 airfoil section; 8, =60° ( b ) NACA 23021 airfoil section; 5, = 60°. (c) NACA 
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Figure 3.29 Maximum section lift coefficients for several NACA airfoil sections 
with double-slotted and split flaps. 


Referring to Equation 3.44 and Figure 3.30, it is obvious that some of the 
superior performance of the double-slotted flap results from the extension of 
the chord. From the figure, (flapped) is equal to 0.12 CJ deg as compared to 
the expected unflapped value of approximately 0.1. Hence, based on the 
actual chord, the increment in for the double-slotted flap is only 1.6. 
However, this is still almost twice that afforded by plain or split flaps and 
points to the beneficial effect of the slot in delaying separation. 

Figure 3.31 (taken from Ref. 3.15) presents pitching moment data for 
flapped airfoil sections. The lift and moment are taken to act at the aerody¬ 
namic center of the airfoil, located approximately 25% of the chord back from 
the leading edge. The moment is positive if it tends to increase the angle of 
attack. 

From Figure 3.31, the lowering of a flap results in an incremental pitching 
moment. In order to trim the airplane a download must be produced on the 
horizontal tail. The wing must now support this download in addition to the 
aircraft’s weight. Hence the effective increment in lift due to the flap is less 
than that which the wing-flap combination produces alone. This correction 
can typically reduce ACi^ by 0.1 to 0.3. 



Figure 3.30 Section lift characteristics of an NACA 63, 4-421 (approximately) 
airfoil equipped with a double-slotted flap and several slot-entry-skirt extensions, 
(a) No skirl extension; R-2.4X10-. ( b) Partial skirt extension; R =2.4x10’. (c) 
Partial skirt extension; R = 6.0x10'. (d) Full skirt extension: R = 2.4x10“. 




































Table 3.1 Maximum Lift Coefficients of Airfoil Sections Equipped with Single-Slotted Flaps 


Slot-entry configuration 



Airfoil Section cy/c 


Slot-Entry 

configuration 


Flap 

Nose Sf 

Shape C,_ (deg) 


Optimum 
y/ Position 


R 


Clark Y 

Clark Y 

Clark Y 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23012 

23021 

23021 

23021 

23021 

23021 

23021 


0.20 1.00 b 

0.30 1.00 b 

0.40 1.00 b 

0.10 0.93 a 

0.15 1.00 b 

0.25 1.00 b 

0.256 0.800 a 

0.257 0.83 a 

0.257 0.83 a 

0.267 1.00 b 

0.30 0.90 c 

0.30 0.90 c 

0.30 0.90 c 

0.30 0.90 b 

0.30 1.00 b 

0.40 0.715 b 

0.40 0.715 a 

0.15 1.00 b 

0.15 1.00 b 

0.25 1.00 b 

0.257 0.827 b 

0.257 0.827 a 

0.257 0.827 b 

0.257 0.827 a 


A 2.44 30 

A 2.83 40 

A 3.10 40 

A 2.25 50 

A 2.68 30 

A 3.22 40 

B 2.76 50 

A 2.81 50 

A 2.83 40 

A 2.90 30 

A 2.92 50 

A 2.92 40 

A 2.93 30 

A 2.88 40 

A 3.29 40 

A 2.87 50 

A 2.90 50 

A 2.59 60 

A 2.66 60 

A 3.17 40 

B 2.69 60 

B 2.74 60 

A 2.71 60 

A 2.82 50 


0 -0.025 

0 -0.025 

0 -0.025 

0.004 0.005 

0 0.015 

0 0.015 

0.005 0.018 

0.005 0.016 

0.013 0.024 

0 0.025 

0.002 0.010 

0.002 0.020 

0.002 0.030 

0.002 0.020 

0 0.015 

0.015 0.015 

0.015 0.015 

0 0.015 

0.050 0.030 

0.025 0.015 

0 0.015 

0 0.015 

0.005 0.020 

0 0.025 


Yes 0.61 

Yes 0.61 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

No 3.5 

No 3.5 

No 3.5 

No 3.5 

No 3.5 

No 3.5 

Yes 3.5 

Yes 3.5 

No 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 


23021 0.40 

23021 0.40 

23030 0.257 

23030 0.257 

23030 0.40 

23030 0.40 

63,4-420 0.25 

3,4-421 (approximately) 0.243 
65-210 0.25 

65-210 0.25 

65-210 0.25 

5(ti2)Alll (approxi- 0.35 

mately) 

5(-213 (approximately) 0.336 
5(215)-114 0.259 

5.2- 221 (approximately) 0.263 

6(215)-116, a =0.6 0.25 

6.2- 116, a — 0.6 0.25 

6.2- 216, a = 0.6 0.30 

6,2-216, a =0.6 0.25 

6.2- 216, a = 0.6 0.25 

6.2- 118 0.25 


0.715 

0.715 

0.860 

0.860 

0.775 

0.775 

0.88 

0.835 

0.84 

0.90 

0.975 

0.839 


0.889 

0.915 

0.832 

0.824 

0.827 

0.90 

0.824 

0.834 

0.90 


A 2.79 50 

A 2.88 50 

B 2.59 60 

B 2.68 60 

B 2.82 50 

B 2.90 50 

B 3.00 35 

A 3.21 40 

A 2.47 45 

A 2.48 41.3 

A 2.45 35 

A 2.69 35 

A 2.63 40 

A 2.80 40 

B 2.83 30 

B 2.70 55 

A 2.69 45 

A 2.92 37 

A 2.89 40 

A 2.88 45 

2.68 32.5 


0.015 0.025 

0.015 0.045 

0.025 0.040 

-0.005 0.040 

0.025 0.060 

0.025 0.060 

0.018 0.045 

0 0.027 

0.009 0.010 

0.014 0.009 

0.004 0.020 

-0.020 0.032 

0.019 0.046 

0.019 0.038 

0.025 0.046 

0 0.028 

0.017 0.038 

0 0.016 

0.023 0.040 

0.011 0.031 


Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

Yes 3.5 

No 6.0 

Yes 9.0 

Yes 6.0 

Yes 6.0 

Yes 6.0 

Yes 9.0 

No 9.0 

No 9.0 

Yes 9.95 

No 6.0 

No 6.0 

No 6.0 

Yes 5.1 

Yes 5.1 

No 6.0 



Typical single—slotted flap configuration, 
me ns ions are given in fractions of airfoil chord.) 


(All dir 





Table 3.2 Maximum Lift Coefficients of Airfoil Sections Equipped with Double-Slotted Flaps 


Airfoil Section 

c,lc 

CjC 

C„ 

C,_, 

s, 

8, 

*r 

yi 

x. 

y» 

Optimum 

Position 

R 

23012 

0.10 

0.189 

0.83 

2.99 

70 

40 

0.009 

0.009 

0.014 

0.024 

Yes 

3.5 X I0‘ 

23012 

0.257 

0.227 

0.715 

3.47 

70 

30 

0.014 

0.012 

0.015 

0.035 

No 

3.5 

23021 

0.257 

0.227 

0.715 

3.56 

60 

30 

0.019 

0.024 

0.025 

0.065 

No 

3.5 

23030 

0.257 

0.260 

0.715 

3.71 

80 

40 

0.049 

0.050 

0.045 

0.040 

No 

3.5 

23012 

0.257 

0.117 

0.826 

3.30 

60 

25 

-0.016 

0.010 

-0.004 

0.017 

Yes 

3.5 

23021 

0.257 

0.147 

0.827 

3.32 

70 

30 

0.017 

0.027 

0.007 

0.024 

Yes 

3.5 

63-210 

0.25 

0.075 

0.84 

2.91 

50 

25 

0.022 

0.024 

0.024 

0.018 

Yes 

6.0 

63,4-421 (approximately) 

0.195 

0.083 

0.87 

3.30 

55 

14 

0.038 

0.012 

0.009 

0.016 

No 

6.0 

64-208 

0.25 

0.075 

0.84 

2.51 

45 

30 

0.015 

0.015 

0.015 

0.019 

Yes 

6.0 

64-208 

0.25 

0.056 

0.84 

2.40 

50 

25 

0.018 

0.014 

0.015 

0.024 

Yes 

6.0 

64-210 

0.25 

0.075 

0.84 

2.82 

55 

30 

0.023 

0.006 

0.012 

0.018 

Yes 

6.0 

64,-212 

0.25 

0.075 

0.84 

3.03 

50 

30 

0.021 

0.020 

0.010 

0.019 

Yes 

6.0 

64,A212 

0.229 

0.083 

0.833 

2.83 

55 

26 

0.044 

0.005 

0.004 

0.014 

Yes 

6.0 

65-210 

0.25 

0.075 

0.84 

2.72 

50 

25 

0.025 

0.011 

0.009 

0.024 

Yes 

6.0 

65(216)-215, a =0.8 

0.248 

0.0% 

0.82 

3.38 

70 

12 

0.024 

0.010 

0.025 

0.032 

No 

6.3 

65 } -118 

0.244 

0.10 

0.864 

3.35 

65 

23 

0.038 

0.007 

0.009 

0.025 

Yes 

6.0 

65j-418 

0.236 

0.106 

0.851 

3.50 

65 

21 

0.027 

0.007 

0.012 

0.028 

Yes 

6.0 

654-421 

0.236 

0.109 

0.85 

3.08 

51 

20 

0.029 

0.017 

0.012 

0.024 

Yes 

2.2 

66-210 

0.25 

0.075 

0.84 

2.64 

55 

25 

0.029 

0.023 

0.012 

0.022 

Yes 

6.0 

66-210 

0.25 

0.100 

0.84 

2.72 

60 

25 

0.027 

0.039 

0.024 

0.021 

Yes 

6.0 

66^-214 (approximately) 

0.227 

0.085 

0.854 

3.00 

55 

20 

0.044 

0.009 

0.004 

0.025 

Yes 

9.0 

1410 

0.25 

0.075 

0.84 

3.06 

50 

25 

0.026 

0.016 

0.012 

0.019 

Yes 

6.0 


1 




Definition of Flap Geometry given in Table 3.2 
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(b) 

Figure 3.31 Influence of flap configuration on pitching moment coefficient, (a) 
Plain flaps, (ft) Split flaps, (c) Slotted flaps. 



Cf 

Flap chord ratio,— 


Figure 3.31 (Continued) 


In a high-wing airplane, lowering the flaps can cause the nose to pitch up. 
This is due to the moment produced about the center of gravity from the 
increase in wing drag because of the flaps. Based on the wing area, the 
increment in wing drag coefficient, ACj> due to the flaps is given ap¬ 
proximately by, 

A C D = 1.7(c//c) U8 (S//S) sin 2 5y (plain and split) (3.45) 

= 0.9(c f icy 3 *(SflS) sin 2 S/ (slotted) (3.46) 

If the wing is located a height of h above the center of gravity, a 
balancing upload is required on the tail. The effect of trim on Cr l l t ^ for a 
complete airplane will be discussed in more detail later. 

Flap Effectiveness in the Linear Range 

Frequently one needs to estimate the increment in C t below stall, AC;, 
produced by a flap deflection. Not only is this needed in connection with the 
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wing lift, but A Q is required also in analyzing the effectiveness of movable 
control surfaces, which frequently resemble plain flaps. 

If an airfoil section has a lift curve slope of C, a and lowering its flap 
produces an increment of A C h the angle of zero lift, is decreased by 

4 “" = Tp (3 - 47> 

The rate of decrease of « 0i per unit increase in the flap angle S f is referred to 
as the flap effectiveness factor, t. Thus, for a flapped airfoil, the lift 
coefficient can be written as 


c, = c la (a + rS f ) (3.48) 

where a is the angle of attack of the airfoil’s zero lift line with the flap 
undeflected. 

Theoretically r is a constant for a given flap geometry but, unfortunately, 
flap behavior with S f is rather nonlinear and hence r must be empirically 
corrected by a factor tj to account for the effects of viscosity. Including tj. 
Equation 3.48 becomes, 

C , = C/ a (a + tij 5/) (3.49) 

The functions r and rj can be obtained from Figures 3.32 and 3.33. Figure 
3.33 is empirical and is based on data from References 3.15, 3.17, 3.19, and 
3.20. Although there is some scatter in the data, as faired, the comparisons 
between the various types of flaps are consistent. The double-slotted flap 
delays separation on the upper surface, so that the decrease in flap effective¬ 
ness occurs at higher flap angles than for the other flap types. The same can 
be said of the slotted flap relative to the plain and split flaps. The plain flap is 
fairly good out to about 20° and then apparently the flow separates from the 
upper surface and the effectiveness drops rapidly, approaching the curve for 
split flaps at the higher flap angles. In a sense the flow is always separated on 
the upper surface of a split flap. Thus, even for smalt flap angles, the effective 
angular movement of the mean camber line at the trailing edge of an airfoil 
with a split flap would only be about half of the flap displacement. 

In the case of the double-slotted flap it should be emphasized that this 
curve in Figure 3.33 is for an optimum flap geometry. The trailing segment of 
the flap is referred to as the main flap and the leading segment is called the 
vane. In applying Equation 3.49 and Figures 3.32 and 3.33 to the double- 
slotted flap, the total flap chord should be used together with the flap angle of 
the main flap. Usually, the deflection angle of the vane is less than that for the 
main flap for maximum lift performance. 

Figure 3.32 is based on the thin airfoil theory represented by Equation 
3.39. As an exercise, derive the expression for r given on the figure, r can 
also be obtained using the numerical methods that led to Figures 3.17 and 



Figure 3.32 Flap effectiveness factor. 


3.18. As another exercise, apply Weissinger’s approximation to the flapped 
airfoil using only two point vortices to represent the airfoil. Placing one 
vortex on the quarter chord of the flap and the other on the quarter chord of 
the remainder of the airfoil leads to 


. 3(3-2 c f )c f 

4(1 - c,)c, + 3 


(3.50) 


where c { is the fraction of chord that is flapped. Equation 3.50 is ap¬ 
proximately 10% lower than Figure 3.32 for crfc values of around 0.25. 

The angle of attack at which the flapped airfoil stalls is generally less than 
that for the plain airfoil. Hence, the increment in because of the flap is 
not as great as the increment in Q at an angle below the stall. Denoting these 
increments by AC^ and AC ( , respectively, it is obvious that the ratio 
ACj^/ACj must depend on crfc. If crfc, for example, is equal to 1.0, in a sense the 
entire airfoil is the flap and AC^ x must be zero. Systematic data on AC^/ACj 
are sparse. Figure 3.34 has been drawn based on a limited number of data 
points and should be used with discretion. 
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Figure 3.33 Correction factor to flap effectiveness factor r. Note that curves 
apply for thickness ratios of approximately 12% and flap chord fractions of 40% or 
less. 



Figure 3.34 C, max increment ratio as a function of flap chord ratio. 


As an example, in using Figures 3.32, 3.33, and 3.34, consider the 
prediction of C, max for a 23012 airfoil equipped with a 30% chord split flap 
deflected 60° and operating at a Reynold’s number of 3.5 x 10 6 . From Figure 
3.32, r = 0.66 for c f lc = 0.3 and from Figure 3.33, *7=0.35 for a split flap 
deflected 60°. Hence from Equation 3.49, A C, is equal to 

AC, = C,tt)8 

= (0.105)(0.66)(0.35)(60) (3.51) 

= 1.46 

In Equation 3.51, C, a of 0.015 is obtained from Reference 3.1. Using Figure 
3.34, the ratio of AC Jmi to AC, is obtained as 0.66. Hence, 

AC U> - 0.96 

According to Figure 3.27, C, mtt for a plain 23012 airfoil equals 1.65 at R = 3.5 
x 10*. Thus, for the flapped airfoil, C, ma> is predicted to be 1.65 + 0.96, or 2.61. 
This result compares closely with Figure 3.26a. If the procedure is repeated 
for other flap angles, close agreement is also obtained with the figure. 
However for a flap chord ratio of 0.1, the predicted values of C,, nas based on 
Figures 3.32 to 3.34 are higher than those shown in Figure 3.26a. 

Leading Edge Devices 

In order to avoid leading edge separation, particularly at low Reynolds 
numbers or for airfoils with relatively sharp leading edges, special high-lift 
devices can also be incorporated into the leading edge to supplement the 
benefits of trailing edge flaps. These are illustrated in Figure 3.35. The fixed 
slot and extensible slat have been in use for some time, whereas the Kruger-type 
nose flap was first employed on the turbojet transport. 

As the name implies, the fixed slot is just that—a narrow channel through 
which the air can flow from the lower surface to the upper surface. This 
channeling of the flow allows the airfoil to operate at higher angles of attack 
before the upper surface of the leading edge separates than otherwise would 
be the case. Increments in C, mt of approximately 0.1 or 0.2 are achieved by 
the fixed slot. It is a moot question as to why this delay in the separation 
occurs. As in the case of slots with trailing edge flaps, the explanation has 
been offered in the past that the flow through the slot feeds energy into the 
slower moving boundary layer, thereby decreasing its tendency to separate. 
More recently, however, in a Wright Brothers’ Lecture (Ref. 3.16) Smith, in 
examining numerical results on multielement airfoils, concluded that im¬ 
proved stall performance from slots is most likely the result of more favor¬ 
able pressure gradients being produced on one airfoil element by the other. 

The extensible slat is similar in its performance to the slot, but it is 
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Kruger-type leading edge flap 

Figure 3.35 Various methods for delaying leading edge separation. 


considerably more efficient because it can be positioned to optimize its 
contribution to C/ roij . The mechanically extended slat is finding increased 
application, particularly with the use of thinner airfoil sections for high-speed 
applications. Figure 3.36 presents some data on slats taken from Reference 
3.17. Here a NACA 64A010 airfoil was tested using a slat in combination with 
split and double-slotted trailing edge flaps. The slat is seen to improve C /otti 
significantly, producing increments in of approximately 0.9, 0.8, and 0.6 
for the no-flap, split-flap, and double-slotted flap configurations, respectively. 
Unlike the trailing edge flap, the primary effect of the slat is seen to be an 
extension of the lift curve without the slat; that is, opening the slat does not 
change C, by a large increment at a fixed angle of attack. The same is true of 
leading edge flaps and is not unexpected in view of Figure 3.32. 

The performance of a leading edge flap is presented in Figure 3.37 for the 
same airfoil as for Figure 3.36. Comparing the two figures, it is obvious that 
the two leading edge devices are nearly comparable in performance. 

Figure 3.38 shows a section of a sophisticated Kruger-type flap. As this 
flap swings down and forward, it assumes the curved shape that is shown. 
With this optimum shaping, its performance probably exceeds to some extent 
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Figure 3.38 Flexible fiberglass leading edge flap used on the Boeing 747 and 
YC-14 airplanes. 

the data presented in Figures 3.39 and 3.40. Figure 3.39 (taken from Ref. 3.18) 
shows that without a trailing edge flap the Kruger flap gives an increment in 
C u , to the 64-012 airfoil of only 0.4. However, the plain airfoil has a higher 

to begin with than that of Figures 3.36 and 3.37. Hence, the total C, max for 
the Kruger-flapped airfoil without a trailing edge flap is about the same as for 
the other two leading edge devices. However, with the split flap, the Kruger 
flap produces a combined C finai equal to 3.0, which is 0.3 to 0.4 higher than the 
corresponding data of Figures 3.36 and 3.37. 

The data of Figure 3.40 (taken from Ref. 3.21) are based on Kruger’s 
original work. 

The Optimum Airfoil for High Lift 

Stratford, in References 3.23 and 3.24, examined both theoretically and 
experimentally the possibility of diffusing a turbulent boundary layer in such a 
way as to produce zero wall shear. Known as “imminent separation pressure 
recovery,” Stratford found that it is indeed possible, with the proper pressure 
gradient, to maintain a velocity profile along a diffuser such that du(y)ldy is 
equal to zero at the wall. u(y) is the velocity in the boundary layer parallel to 
the wall and is a function of the distance, y, from the wall. With the velocity 
gradient at the wall equal to zero, the boundary layer is just on the verge of 
separating, since a negative value of this gradient will result in reverse flow, 
as illustrated in Figure 3.41. 
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Section angle of attack q 0 , deg 

Figure 3.39 Section lift characteristics for the NACA 64i-012 airfoil section 
equipped with a 0.10c upper-surface leading edge flap alone, and in combination 
with a 0.20c trailing edge split flap. R = 6.0 x 10®. 

In the abstract to Reference 3.24, Stratford states: 

“No fundamental difficulty wax encountered in establishing the flow and it 
had, moreover, a good margin of stability. The dynamic head in the zero skin 
friction boundary layer was found to be linear at the wall (i.e., u « y 1G ), as 
predicted theoretically in the previous paper. (Author's note, Stratford is 
referring to Ref. 3.23.) 


Figure 3.40 Characteristics of Kruger flaps, (a) Illustration of Kruger's type of 
nose flap and simple, hinged, leading edge flap, (b) Effect of flap angle on 
maximum lift coefficient, (c) Effect of flap chord on maximum lift coefficient. 
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The flow appears to achieve any specified pressure rise in the shortest 
possible distance and with probably the least possible dissipation of energy for a 
given initial boundary layer. Thus, an airfoil which could utilize it immediately 
after transition from laminar flow would be expected to have a very low drag." 

The Stratford imminent separation pressure recovery was adopted for 
airfoils by Liebeck and Ormsbee (Ref. 3.25) and was extended later by 
Liebeck (Ref. 3.22). Using variational calculus, optimum chordwise pressure 
distributions for the upper and lower surfaces are prescribed that are modified 
slightly by additional constraints not present in the optimization process. 
Specifically, the optimum C p distributions are modified in order to (1) close 
the airfoil contour at the trailing edge, (2) round the leading edge to allow 
operation over an angle-of-attack range, and (3) satisfy the Kutta condition at 
the trailing edge. 

The resulting modified form of the optimum pressure distribution is 
compared with the optimum distribution in Figure 3.42. Beginning at the 
stagnation point, the flow is accelerated up to a so-called rooftop region over 
which the velocity, and hence the pressure, is constant. Following the rooftop 
region, the Stratford pressure recovery distribution is employed to reduce the 
velocity over the upper surface to its value at the trailing edge. 

One such airfoil design is presented in Figure 3.43 (taken from Ref. 3.22). 
Included on the figure is the pressure distribution to which the airfoil was 
designed. Test data on this airfoil obtained at a Reynolds number of 3 x 10‘ 
are presented in Figure 3.44a and 3.441?. Although this configuration is 
referred to by the reference as a-“turbulent rooftop” case, transition does not 
occur until the start of the Stratford pressure recovery. In this case the 
performance of the airfoil is seen to be good from the standpoint of Ci mx and 
Cj. The drag coefficient remains below the value of 0.01 over a wide range of 
Q values from 0.6 to 1.6. 

Artificially producing transition near the leading edge severely com¬ 
promises C Jmii and Cd, as shown in Figure 3.44b. Still, by comparison with the 
standard NACA airfoils, the Liebeck airfoil appears to offer superior per¬ 
formance at low speeds and, in the future, may find application to general 
aviation aircraft. One possible drawback in this regard is the sharp drop in its 
lift curve at stall. 

Powered-Lift Systems 

Figure 3.45 (taken from Ref. 3.26) presents the growth of C Ltoi over the 
years since the Wright Brothers’ success. The two points labeled K and L are 
somewhat misleading, since these two aircraft were experimental in nature 
and used distributed suction over the wing to delay separation. From this 
figure and the preceding information on flapped and plain airfoils, of 




Figure 3.42 Pressure distributions for Liebeck airfoils, (a) Optimum airfoil pres¬ 
sure distribution according to variational analysis, (b) Modified form of optimum 
pressure distribution for airfoil upper and lower surfaces (not to scale). (R. H. 
Liebeck and A. I. Ormsbee, "Optimization of Airfoils for Maximum Lift", AIAA 
Journal of Aircraft, 1970. Reprinted from the Journal of Aircraft by permission of 
the American Institute of Aeronautics and Astronautics.) 

slightly over 3 is probably the best that can be achieved without the addition 
of power. Although two-dimensional airfoils with double-slotted flaps can do 
better than this, as will be seen later, their full potential cannot be achieved 
when applied to an airplane. Generally, the flaps cannot be applied over the 
entire span of the wing. In addition to this loss in an added penalty 
results from the fuselage and tail download required for trim. 

Cl^ values considerably above those achievable with flap systems dis¬ 
cussed so far are possible by the expenditure of power. Most of the powered- 
flap systems presently under consideration bear a resemblance, or can be 
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Figure 3.43 Liebeck airfoil with its pressure distribution. (R. H. Liebeck, ‘Class 
of Airfoils Designed for High Lift”, AIAA Journal of Aircraft, 1973. Reprinted from 
the Journal of Aircraft by permission of the American Institute of Aeronautics and 
Astronautics.) 

related in their performance, to the jet flap. Thus, in order to understand 
better the performance of systems such as upper surface blowing (USB), 
externally blown flaps (EBF), augmentor wing, and circulation control, we 
will begin with the jet flap shown in Figure 3.46. A thin sheet of air exits the 
trailing edge at a downward angle of 5 relative to the airfoil zero lift line. This 
line is shown at an angle of attack a. If r c is the total circulation around the 
airfoil then, assuming a and 8 to be small, the total lift on the airfoil will be 

L = p vr c + mjV t (a + 8) (3.52) 

where ntj is the mass flux in the jet and Vj is the jet velocity. 

As the jet leaves the airfoil, it gets turned in the direction of the 
free-stream velocity. In order to redirect the flux of jet momentum, it follows 
that a pressure difference must exist across the jet. This pressure difference, 
Ap, can be related to m } Vj and the radius of curvature R by the use of Figure 
3.46b. Applying the momentum theorem in the direction of curvature, 

ApRd = rrijVjd 
or 

(3.53) 



Figure 3.44a Experimental lift curve and drag polar for airfoil of Figure 3.43. (R. 
H. Liebeck, "Class of Airfoils Designed for High Lift”, AIAA Journal of Aircraft, 
1973. Reprinted from the Journal of Aircraft by permission of the American 
Institute of Aeronautics and Astronautics.) 



Figure 3.44b Effect of transition location on the lift and drag of a Liebeck airfoil. 
(R. H. Liebeck, “Class of Airfoils Designed for High Lift”, AIAA Journal of Aircraft, 
1973. Reprinted from the Journal of Aircraft by permission of the American 
Institute of Aeronautics and Astronautics.) 
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Figure 3.45 History of maximum lift coefficients for mechanical lift systems. (F. 
A. Cleveland, “Size Effects in Conventional Aircraft’’, AIAA Journal of Aircraft, 
1970 Reprinted from the Journal of Aircraft by permission of the American 
Institute of Aeronautics and Astronautics.) 
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Since the jet exerts a force on the fluid, it can be replaced by an 
equivalent continuous vortex sheet that exerts the same force. Letting y } be 
the strength per unit length of the sheet. 


or 


A pR0 — pVyRd 


(3.54) 


Measuring the position of the jet, y, positively downward, the radius of 
curvature and y for a nearly horizontal jet are related by 


R 


d 2 y 

dx 2 


(3.55) 


Combining Equations 3.53, 3.54, and 3.55 gives 

_ -mjVj d 2 y 
y ‘~ pV Hx- 


(3.56) 


Equation 3.56 relates the jet vortex strength to the shape of the sheet and 
the jet momentum flux. 

The total circulation of the jet vortex sheet can be obtained by integrating 
Equation 3.56 from x = 0 to ». 


_ mjVj dy f” 
pV dx Jo 


= (a + 5) (3.57) 

Combining Equations 3.52 and 3.57 shows that the Kutta-Joukowski rela¬ 
tionship holds for the jet-flapped airfoil if the circulation is taken around both 
the airfoil and the jet. 

L = p^(r c + r / ) (3.58) 

The boundary value problem is then posed where the airfoil and jet sheet 
are each replaced by an unknown vortex distribution. Distributions must then 
be found that will induce a velocity at each point on the airfoil and combining 
with the free-stream velocity to give a resultant velocity tangent to the airfoil. 
Along the sheet the following must hold. 


The details of the solution are beyond the scope of this text and can be 
found in Reference 3.19. 

Although the solution of Reference 3.19 is not in closed form, the results 
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can be expressed in a relatively simple way. As with a physical flap, the 
increment in C, because of a change in angle of attack and flap angle can be 
expressed as a linear combination of the two angles. 


C, = C\a + C, s S 



(3.59) 


C k 


dC, 
d 8 


The derivatives C, a and C, a are a function of the ratio of the jet 
momentum flux to the product of the free-stream dynamic pressure and a 
reference area. This ratio, known as the momentum coefficient, C„, is defined 
for a two-dimensional airfoil by 


If, in addition, the jet exits ahead of the trailing edge and blows over and 
is deflected by a physical flap having a chord of q (see the blown flap of 
Figure 3.24), then C lg is also a function of crfc. For a pure jet flap (C,/C = 0), 
C;„ and C|, are given by 


C, t = [4ttC;(1 + 0.151C m ,q + 0.139C ( J] 1/2 (3.61) 

~ 2ir(l + O.I51Q 1/2 + 0.219Q) (3.62) 

For qlc values other than zero, C is is given in Figure 3.47. The curve 
labeled qlc =1.0 in this figure corresponds to Equation 3.62 since, for this 
case, C, s = C t . 

Data concerning Ci mia for jet-flapped airfoils is sparse. Generally, the jet 
flap follows the predictions of Figure 3.41 fairly closely, since the jet fixes the 
Kutta condition and provides some control over the boundary layer to 
prevent separation. As a preliminary estimate for C,^, Reference 3.3 
recommends the use of the relationship presented in Figure 3.48. Here the 
difference in the angle of attack for stall, with and without blowing, is 
presented as a function of C M . 

The negative pitching moment of the jet flaps is higher than the moment 
for conventional flaps for two reasons. First, the jet reaction acts at the 
trailing edge; second, the jet-flapped airfoil maintains lift all the way back to 
t e trading edge. As with the lift, C M can be written as a linear combination of 
a and S. 

C M = C Ma a + C Mi 8 


In this equation dCjJdCi can be obtained from Figure 3.49. 



Figure 3.47 Blown flap effectiveness. (B. W. McCormick, Aerodynamics of 
V/STOL Flight, Academic Press, Inc. 1967. Reprinted by permission of Academic 
Press, Inc. 



Figure 3.48 Effect of on a for C w (B. W. McCormick, Aerodynamics of 
V/STOL Flight, Academic Press, Inc. 1967. Reprinted by permission of Academic 
Press, Inc.) 
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Figure 3.49 Effect of on pitching moments for blown flaps. (B. W. McCormick, 
Aerodynamics of V/STOL Flight, Academic Press, Inc. 1967. Reprinted by per¬ 
mission of Academic Press, Inc.) 


To illustrate the use of the foregoing relationship for the jet-flapped 
airfoil, consider the prediction of for the NACA 63,4-421 airfoil (Figure 
3.30) equipped with a 25% blown flap deflected 50°. The jet expands isen- 
tropically from a reservoir pressure of 25 psia and a temperature of 70 °F. The 
airfoil is operating at 50 mph at SSL conditions. It has a chord of 5 ft and the 
jet thickness is 0.2 in. 

We begin by calculating the jet velocity from the compressible Bernoulli 
equation (Equation 2.31). 

Bj >=J3L(a_£) 

y -1 \po p/ 


or, with the use of Equation 2.30, 

•>-feSte[-(£ 

From the equation of state (Equation 2.1), 
Pol Po = RT 0 

Thus, 

„ _ f 2(1.4)(1716)(529.7) __ 


©t 


= 946.9 fps 

The mass density of the expanded jet will be 
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p 0 is calculated from the equation of state. 

25(144) 

1716(529.7) 

- 0.00396 slugs/ft 3 

Thus , 

p, =0.00271 slugs/ft J 

The jet mass flux will be equal to 

m i = PAi v l 

= 0.00271(0.2/12)(946.9) 
= 0.0428 slugs/s 

The free-stream dynamic pressure is 

q = {pV 2 

= 0.0023378(50 xl.467) 2 /2 
= 6.397 psf 


_ (0.0428X946.9) 

6.397(5) 

= 1.27 

From Figures 3.47 and 3.48, 

C, a = 9.09/rad 
C h = 6.5/rad 
- <w(c; = o) = -2° 

From Figure 3.30, C fmai for the unblown airfoil without a flap is equal 
approximately to 1.4. Using Figures 3.33 and 3.34, AC, due to the plain flap 
(C; = 0) is estimated as 


AC, = 2 7T (0.6)(0.035) | 




Thus, for Q = 0 with the flap deflected, C^ M is estimated to be 2.17. At a = 0, 
C, is estimated to be 1.5. Thus, with Ct a = 0.109 (from Figure 3.30), 


= 6.T 
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or, relative to the zero lift line, flaps up, 

a max = 9.1° 

For the operating C M , the angle of attack for the zero lift line at stall is 
estimated to equal 7.1°. 

Thus, 

Ci m „ ~ Gjohm* + Ct, s 


The preceding answer must, of course, be further corrected, using Figure 
3.43 and Equation 3.46, to account for trimming tail loads. Also, it should be 
emphasized that the preceding is, at best, an estimate for preliminary design 
purposes or relative parametric design studies. In the final analysis, model and 
prototype component testing of the blowing system must be performed. 

Credit for the practical application of the jet flap must be given to John 
Attinello. Prior to 1951, all blown systems utilized pressure ratios less than 
critical in order to avoid supersonic flow in the jet. Such systems required 
large and heavy ducting. For his honors thesis at Lafayette College, Attinello 
demonstrated with “homemade” equipment that a supersonic jet would 
adhere to a deflected flap. This was contrary to the thinking of the day that 
not only would a supersonic jet separate from a blown flap, but the losses 
associated with the shock wave system downstream of the nozzle would be 
prohibitive. Later, more sophisticated testing performed by the David Taylor 
Model Basin confirmed Attinello’s predictions (Ref. 3.38) of high lift 
coefficients for supersonic jet flaps. This led to the development of compact, 
lightweight systems using bleed air from the turbojet engine compressor 
section. The Attinello flap system was flight tested on an F9F-4 and produced 
a significant decrease in the stalling speed for an added weight of only 50 lb. 
Following this success, the Attinello flap went into production on the F-109, 
F-4, F8K, AS, and other aircraft, including several foreign models. 


THE LIFTING CHARACTERISTICS OF A FINITE WING 

A two-dimensional airfoil with its zero lift line at an angle of attack of 10° 
will deliver a lift coefficient, C h of approximately 1.0. When incorporated into 
a wing of finite aspect ratio, however, this same airfoil at the same angle of 
attack will produce a wing lift coefficient, Cl, significantly less than 1.0. The 
effect of aspect ratio is to decrease the slope of the lift curve C La as the 
aspect ratio decreases. Figure 3.50 illustrates the principal differences in the 



Figure 3.50 Comparison of NACA 65-210 airfoil lift curve with that of a wing 
using the same airfoil. 

lift behavior of a wing and an airfoil. First, because the wing is twisted so that 
the tip is at a lower angle of attack than the root (washout), the angle for zero 
lift, measured at the root, is higher for the wing by approximately 0.6°. Next, 
the slope of the wing’s lift curve, C L „, is approximately 0.79 of the slope for 
the airfoil. Finally, C^ is only slightly less than C lmix , in the ratio of 
approximately 0.94. These three differences are almost exactly what one 
would expect on the basis of wing theory, which will now be developed. 

The Vortex System for a Wing 

A wing’s lift is the result of a generally higher pressure acting on its lower 
surface compared with the pressure on the upper surface. This pressure 
difference causes a spanwise flow of air outward toward the tips on the lower 
surface, around the tips, and inward toward the center of the wing. Combined 
with the free-stream velocity, this spanwise flow produces a swirling motion 
of the air trailing downstream of the wing, as illustrated in Figure 3.51. This 
motion, first perceived by Lanchester, is referred to as the wing’s trailing 
vortex system. 

Immediately behind the wing the vortex system is shed in the form of a 
vortex sheet, which rolls up rapidly within a few chord lengths to form a pair 
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Figure 3.51 Generation of vortex system by finite aspect ratio wing. 


of oppositely rotating line vortices. Looking in the direction of flight, the 
vortex from the left wing tip rotates in a clockwise direction; the right tip 
vortex rotates in the opposite direction. 

The trailing vortex system, not present with a two-dimensional airfoil, 
induces an additional velocity field at the wing that must be considered in 
calculating the behavior of each section of the wing. 

If the aspect ratio of the wing is large, approximately 5 or higher, the 
principal effect of the trailing vortex system is to reduce the angle of attack of 
each section by a small decrement known as the induced angle of attack, a j. 
In this case Prandtl’s classical lifting line theory (Ref. 3.28) applies fairly well. 
As shown in Figure 3.52, the wing is replaced by a single equivalent vortex 
line, known as the “bound vortex,” since it is in a sense bound to the wing. 
The strength of this vortex, F(y), is related to the lift distribution along the 
wing by the Kutta-Joukowski relationship. 


dL{y) 

dy 


(3.64) 



Figure 3.52 Lifting line model of a wing and trailing vortex system. 

Expressing the lift per unit length of span in terms of the section chord 
length, c(y), and section lift coefficient, C<(y), leads to 

r(y) = My)C,(y)V (3.65) 

With no physical surface outboard of the wing tips to sustain a pressure 
difference, the lift, and hence I\ must vanish at the tips. 

According to the Helmholtz theorem regarding vortex continuity (Ref. 
1.3, p. 120), a vortex line or filament can neither begin nor end in a fluid; 
hence it appears as a closed loop, ends on a boundary, or extends to infinity. 
Thus, it follows that if in going from y to y + dy the bound circulation around 
the wing increases from T to T + d F, a free vortex filament of strength dT, 
lying in the direction of the free-stream velocity, must be feeding into T in 
order to satisfy vortex continuity. This statement may be clarified by 
reference to Figure 3.53. 

The entire vortex system shown in Figure 3.52 can be visualized as being 
closed infinitely far downstream by a “starting” vortex. This vortex, opposite 


r r + dv 



Figure 3.53 Illustration of vortex continuity. 
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in direction to the bound vortex, would be shed from the trailing edge of the 
wing as its angle of attack is increased from zero. 

The trailing vortex system of strength dr induces a downwash, w(y), at 
the lifting line, as mentioned earlier. As pictured in Figure 3.54, this reduces 
the angle of attack by the small angle a-,-. Thus the section lift coefficient will 
be given by 

C, = C 1(i (a-a.) (3.66) 

a being measured relative to the section zero lift line. 

To a small angle approximation, the induced angle of attack, a h is given 
by w/V. The downwash, w, can be determined by integrating the contributions 
of the elemental trailing vortices of strength dr. If the vortex strength dr 
trails from the wing at a location of y, its contribution to the downwash at 
another location y 0 can be found by Equation 2.64 to be 


Thus, a, becomes 


dw(y e ) = 


ar(y) 

4ir()'o-y) 


“M-dvC 


(3.67) 

(3.68) 


Equations 3.65, 3.66, and 3.68 together relate T(y) to c(y) and a(y) so 
that, given the wing geometry and angle of attack, one should theoretically be 
able to solve for T and hence the wing lift. In order to accomplish the 
solution, it is expedient to make the coordinate transformation pictured in 
Figure 3.55. 



*t<y* 


Figure 3.54 A wing section under the influence of the free-stream velocity and 
the downwash. 
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Figure 3.55 Coordinate transformation used in the reduction of the lifting line 
equations. 


Hence, Equation 3.67 becomes 

“'W-d^ Losa-c U < 369 > 

Since more elaborate and comprehensive treatments of wing theory can 
be found in texts devoted specifically to the subject (e.g., see Ref. 3.29), only 
the classical solution for the eliptic T distribution will be covered here. This 
particular case is easily handled and results in the essence of the general 
problem. 

Assume that T is of the form 

r = r 0> /1 - (3.70) 


Using Equation 3.68 this transforms to 

T - T 0 sin $ 

Here, r 0 is obviously the midspan value of the bound circulation. Thus 
Equation 3.69 becomes 

rin r 0 f* cos ede 
ot 2irbVJo cos 0- cos 0o 

The preceding integral was encountered previously in thin airfoil theory 
and has a value of it. Thus, for an elliptic T distribution, a,- and hence the 
downwash is found to be a constant independent of y. 


If the wing is untwisted so that a is also not a function of y then, from 
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Equation 3.66, it follows that the section C ( is constant along the span. Thus, 


= bv 2 C, cdy 
4 J-U2 


Equation 3.65 then gives 

r (y) = ^c(y)C L V 
or. 

Thus it is found that, according to lifting line theory, an untwisted wing with 
an elliptical planform will produce an elliptic T distribution. Such a wing will 
have a constant downwash and section C<. 

Since C, is constant and equal to C L , Equations 3.65, 3.66, and 3.71 can be 
applied at the midspan position in order to determine the slope of the wing lift 
curve. First, from Equations 3.72 and 3.65, 


But, for the planform given by Equation 3.66, c 0 and b are related to the 
aspect ratio by 


Thus, becomes 


Inserted into Equation 3.66, the preceding results in 


C L = C ta 
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or, 

Cl = c 4(i + coi*a \" 

Using the theoretical value of 2 ir Cj/rad derived earlier, the preceding 
becomes 



Equations 3.72 and 3.73 are important results. The induced angle of 
attack is seen to increase with decreasing aspect ratio which, in turn, reduces 
the slope of the lift curve, C^. A wing having a low aspect ratio will require a 
higher angle of attack than a wing with a greater aspect ratio in order to 
produce the same C L . 

It was stated previously that the comparative performance between the 
wing and airfoil shown in Figure 3.50 could be explained theoretically. In this 
case, A = 9.02 so that, on the basis of Equation 3.73, 

C La = 0.819 C,„ 

This result is within about 2% of the experimental results presented in Figure 


As the aspect ratio decreases, the lifting line becomes progressively less 
accurate. For example, for an aspect ratio of 4.0, Equation 3.73 is ap¬ 
proximately 11% higher than that predicted by more exact methods. 

As described in Reference 3.3, a more accurate estimate of C La is 
obtained from 


- A + [2(A + 4)/(A + 2)] 


An alternate to Equation 3.74a is offered by Reference 3.35 and is 
referred to as the Helmbold equation, after the original source noted in the 
reference. The Helmbold equation reads 


‘(CJn) + V(CJir) 2 ^A 2 


Replacing C ta by 2ir in the denominator. 


Equation 3.74a and 3.74b agree within a couple of percent over the range of 
practical aspect ratios and approach each other in the limits of A = 0 or 
A = oo. This holds for high or low aspect ratios and is based on an approximate 
lifting surface theory, which accounts for the chordwise distribution of bound 
circulation as well as the spanwise distribution. 
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The Maximum Lift of a Finite Wing 

The maximum lift coefficient of a finite wing is influenced by several 
factors. Obviously, C^, is strongly dependent on Ci^; that is, the wing’s 
performance depends on it’s airfoil performance. Second, the spanwise extent 
to which the wing is flapped has a significant influence on Cl^. Also in 
estimating Cl^, one must account for the presence of the fuselage, the tail 
download required to trim the aerodynamic pitching moment, and the span- 
wise distribution of loading over the wing. 

The effect of aspect ratio on is slight, as one might expect from the 
preceding considerations on the elliptic wing. The wing lift coefficient and 
section lift coefficients are nearly equal. 

The detailed estimation of a wing’s C LmtK begins with a calculation of its 
spanwise load distribution. There are several methods to be found in the 
literature for doing this. Many of these fall into a class known as vortex 
lattice methods. One of the first of these can be found in Reference 3.30. 

The vortex lattice method is similar to lifting line theory except that 
discrete vortex lines are also distributed in the chordwise direction. As 
illustrated in Figure 3.56, the wing is covered with a mesh of spanwise and 



chordwise bound vortex lines with free vortex lines shed downstream from 
the trailing edge. At every juncture of the vortex lines, Helmholtz’s law of 
vortex continuity must hold. In addition, control points are chosen over the 
wing equal in number to the number of unknown vortex line strengths. The 
unknown vortex strengths are then adjusted to assure that the resultant flow 
at each control point is tangent to the mean surface of the wing. 

The lattice model shown in Figure 3.56, in the manner of Reference 3.30, 
is formed by the superposition of horseshoe-shaped line vortices. One such 
vortex is shown dashed for illustrative purposes. Note that the downstream 
control points are aft of the last bound vortex lines; some points are even 
slightly off of the wing’s surface. Stipulating that the flow at these points 
parallels the camber surface at the trailing edge will satisfy approximately the 
Kutta condition. 

Trends in the behavior of C^ as related to wing geometry can be seen 
by the application of an approximate method for determining spanwise load 
distribution known as Schrenk’s approximation (Ref. 3.31). 

This method deals with two distributions: a basic lift distribution, and an 
additional lift distribution. The basic distribution is that which exists along the 
span when the total lift is equal to zero. Approximately, this lift distribution is 
taken as the average of a constant “zero” distribution and one obtained by 
neglecting any induced velocities. Thus, 

l 2 pV 2 cC, = \ [\pV 2 cCi a (a wv + e) + 0] 
or the basic section lift coefficient, C; fc , will be 


C tb = -y- (<*«* + «) 


(3.78) 


The additional lift distribution results from an angle of attack different 
from and is assumed to be the average of an elliptic distribution and one 
proportional to the planform, both having the same total lift. For the latter, 

cCi « c 
= kc 


But 


L= I ^ qcC, dy 
= qk I cdy 


Thus, k = C L . 

For the elliptic distribution, 
cC L 


qks 




Figure 3.56 Vortex lattice model for a wing. 
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Usually Ci, is defined as the value of Equation 3.79 for a C L of unity; thus, 
Ci = C 4 + C,C. L (3.80) 

The manner in which this equation is used to estimate Cl. is best 
explained by an example. Consider the wing of Figure 3.50. This particular 
wing has a taper ratio of 0.4 and a washout of Y. Using Equations 3.78 and 
3.79, the basic and additional section lift coefficient distributions given in 
Figure 3.57 were calculated. Also graphed on this figure is C, as a function 
of spanwise location. In this instance, C,_, is taken from Figure 3.50 to be a 

re°l!c an | L i" many cases ’ leases toward the tip as the airfoil becomes 
relatively thinner or as the chord lengths become smaller. 

Combining C,, and C in the form of Equation 3.80, Figure 3 57 shows 

ttZaf; l 1 - 22 r r lts ; n a sec,ion c - haif - y «* 5S 

6 n eC ‘.° n , c u. at that location. Any attempt to increase C L 
above this value will therefore cause the wing to stall at this location. Since 
the Ci curve is rather flat in this location, the stalling would be expected to 
spread either side of 2„b equal to 0.5. Thus, to^stimate fteT“ „f“ 
g, one finds the wing C L that results in a section C, somewhere afeng the 
on M com ' S )U r‘ equal “ the Secti ° n C ^“- ,n this insta "cc, the C ^ value 
the C nred re , S dT ra u! y W “ h ‘ he experimental results. Generally, however, 
total iinn e d ^-n ‘*" S method wl11 be somewhat conservative, since the 
to loan Cl may * mcrease somew hat, even though a section of it begins 

of Fiture f 3 r ™ er “ d T* extreme exam P k of the method, consider the wing 

From F ; e , qu,pped Wlth 60% Span ’ 20% chord split Hops deflected 60° 
From Equation 3.43, AC, is estimated to be 


AQ = C^tt] 8 

= 0.108 (0.545)(0.35)(60) 
= 1.236 


ACl ™" empirically, is approximately 0.83 of the preceding equation (Figure 


THE LIFTING CHARACTERISTICS OF A FINITE WING 143 



Figure 3.57 Predicted lift coefficient distributions for wing of Figure 3.50. 


3.34) or 1.026. Thus, Figure 3.57 is revised as shown in Figure 3.58 to increase 
the section C (mM to 2.30 over the inner 60% of the span. 

Aerodynamically the twist of the wing is changed by lowering the flaps. 
Relative to the midspan chord, the zero lift lines of the sections outboard of 
the flaps are reduced in angle of attack by AC,/C, a , or 11.4°. Thus, for this 
flapped wing, e, in degrees, becomes 



For this twist distribution and a taper ratio of 0.4, the angle of attack of the 
midspan zero lift line, a^, for zero lift becomes 4.24°. Thus 

C,„ = 0.054(4.24 + e) 

The additional lift distribution remains the same, since the planform is 
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Figure 3.53 Predicted lift coefficient distributions for wing of Figure 3.50 with 
60% span, 20% chord split flaps deflected 60°. 

unchanged. The predicted C; distributions with the partial span split flaps are 
presented in Figure 3.58. 

Figure 3.58 predicts that the wing will begin to stall just outboard of the 
flaps at a wing C, of 1.63. This result agrees exactly with Reference 3.27 with 
regard to both and the location of the initial stall. This agreement is 

somewhat fortuitous in view of Shrenk’s'approximation, which is obviously 
inexact, since it allows a finite loading at the tip and other discontinuities in 
the cCi distribution. Nevertheless, for preliminary design studies, or in lieu of 
more exact lifting surface methods, Shrenk’s approximation is a useful tool. 

Effect of Fuselage on C ^ 

In working with a wing-fuselage combination, one normally defines the 
wing planform area to include the portion submerged within the fuselage. 
When a lift coefficient is quoted for the combination, it is based on this total 
wing planform area obtained by extrapolating the leading and trailing edges 
into the fuselage centerline. Generally, the fuselage will effect a decrease in 



Figure 3.59 Effect of fuselage on spanwise lift distribution. 


the lift per unit span over the portion of the wing covered by the fuselage. 
This is illustrated in Figure 3.59. The upper distribution is without the 
fuselage. The dashed line on the lower figure is the qualitative drop in cC, due 
to the fuselage. As an approximation, let us assume that the fuselage effects a 
constant drop in cC, over its width proportional to the midspan value of cQ. 
Thus, the lift decrement resulting from the fuselage will be 

A L = —kqC^Sfusi 

Stax is the wing planform area submerged in the fuselage, and C^ is the nearly 
constant section Q near the center of the wing, k is the constant of 
proportionality. Thus the total C L with the fuselage, can be written in 
terms of C L before the fuselage is added as 

(3.81) 

In Reference 3.27, two wings equipped with partial and full-span, split, 
single-slotted, and double-slotted flaps were tested with and without a 
fuselage. The fuselage was circular in cross-section and the wing was mounted 
slightly above the middle of the fuselage. The ratio S t ,..J S was equal to 0.083. 
The results of these tests are plotted in Figure 3.60 and compared with 
Equation 3.79 using kCJCt,= 1.0. Also plotted on Figure 3.60 are test results 
from References 3.32 and 3.33. The ratio S tua JS was nearly the same for these 
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Figure 3.60 Effect of fuselage on C Lmax . 

two references as for Reference 3.27. These data support the form of 
Equation 3.81, at least to the extent that the correction to Q w for the 
fuselage appears to increase linearly with C ^ of the wing alone. The 
correction depends on the cross-sectional shape of the fuselage and seems to 
vanish or even be slightly favorable for a rectangularly shaped section. 
Reference 3.34 also shows the correction to be slight for elliptical shapes 
where the height is greater than the width. 

The decrement in Cl^ x also depends on wing position and appears to be a 
maximum for the midwing configuration. 

Effect of Trim on 

In order to calculate the stalling speed of an airplane in steady flight, one 
must consider that, in addition to the weight, the wing’s lift must support any 
download on the horizontal tail required to trim the airplane around its 
pitching axis. In order to determine this additional trim load, refer to Figure 



Figure 3.61 Longitudinal trim of an airplane. 

3 . 61 . Here, the wing lift L, the tail lift, L T , the pitching moment about the 
wing’s aerodynamic center, M K , and the weight are all shown in a positive 
sense. With the aerodynamic center of the tail located a distance of l T behind 
the center of gravity and the wing’s aerodynamic center a distance of 
(x cg -x ac )c ahead, the tail lift to trim is given by 

. Mac . j C . v 

L T = ~j~ + Ljix ci - x ac ) 

In addition, static equilibrium in the vertical direction requires that 
L+L t = W 

Therefore, it follows that 

W = l[i + |(x«-x„)] + ^ 

In coefficient form this becomes 

Cl -- C Lw [ 1 + - l.j] + C^j- T ( 3 . 82 ) 

Here, C L is taken to mean the trim C L . 



C f . w refers to the untrimmed wing lift coefficient corrected for the 
fuselage. 

It was mentioned earlier that the added drag caused by flaps must 
sometimes be considered in the trim of an airplane. If ACj> denotes this 
increment in the drag coefficient and if the flaps are located a distance of h 
above the center of gravity, Equation 3.82 modified to account for the flap 
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drag becomes 

Cl = G, [ 1 + j- (x ct - *«)] + C M „ j- + AG, | (3.83) 

A Co can be obtained experimentally or estimated on the basis of Equations 
3.45 and 3.46. 

C M „ is normally negative and greater in magnitude than the moments 
resulting from C Lw and ACx> Thus Cl is normally less than C Lw . Since 
Equation 3.83 holds for maximum lift conditions, it follows that the trim CV 
is normally less than the wing Ci™- 

In calculating Cu^ for use in Equation 3.83, the section aerodynamic 
moment determined from Ci^ . including the increment because of the flaps, is 
integrated over the wing excluding the part submerged in the fuselage. 

Estimation oi C L for a Complete Airplane Configuration 

A Piper Cherokee PA-28 is pictured in Figure 3.62. Pertinent dimensions, 
areas, weights, and other data are tabulated on the figure. Extrapolating the 




swept leading edge near the root into the fuselage centerline and accounting 
for the elliptically shaped tips gives a total wing area when the flaps are 
extended of 165.1ft 2 . The area of the wing within the fuselage is 25.3 ft 2 . 
Assuming beforehand, or by iteration, a reasonable value for the “stalling 
speed” of 60 mph leads to a Reynold’s number of approximately 3 x 10* for a 
wing section. For this Reynold’s number. Reference 3.1 shows a value for 
C/nw, 1-45 for the plain 65 2 -415 airfoil with a lift curve slope of 0.106 CJde g. 

Using an 18.5% chord, single-slotted flap deflected 40°, Figures 3.32, 3.33, 
and 3.34 predict a AC, of 1.37 corresponding to an increase of 12.9° in the 
angle of attack of the zero lift line. AC (fras is estimated at 1.33 giving a Q mM of 
2.78 for the flapped wing sections. The derivative dCJdQ is estimated from 
Figure 3.31 to equal -0.20. Since Cm„ c ~ -0.07 for the plain airfoil (Ref. 3.1), 

— —0.34 for the flapped airfoil. 

Accounting for the 2° of washout and the increment of a 0 » caused by flaps 
leads to an value of 6.0° from Equation 3.76. This is the angle of attack of 
the zero lift line at midspan for a zero wing C L . C, a and Ct b can then be 
calculated for the wing alone and are given in Table 3.3. The section 




PA-28 Wing with Flaps Down 40 




aerospace engineering students at The Pennsylvania State University as one 
of the experiments in a course on techniques of flight testing. The other two 
points were calculated from the stalling speeds quoted by the manufacturer in 
the airplane’s flight manual. 
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were obtained at R values of 3 x 10 6 and higher. For remote-piloted vehicles 
(RPV), model airplanes, and the like, Reynolds numbers as low as 5 x 10 4 can 
be encountered. A search of the literature will show little airfoil data available 
in this Reynolds number range. The most reliable low-Reynolds number 
airfoil data appear to be those given in Reference 3.35, where tests of five 
different airfoil shapes are reported for R values as low as 42,000. These tests 
were conducted in a low-turbulence tunnel. 

The five airfoil shapes that were tested in Reference 3.37 are shown in 
Figure 3.64. These are seen to comprise a thin, flat plate, a thin, cambered 
plate, two 12% thick airfoils with 3 and 4% camber, and one 20% thick airfoil 
with 6% camber. The airfoil shapes are similar in appearance to the NACA 
four-digit series. 

The lift curves for these airfoils are presented in Figure 3.65 for four 
different Reynolds numbers. As one might expect, the flat-plate results are 
nearly independent of R since the separation point at the leading edge is well 
defined. To a slightly lesser degree, the same can be said for the cambered 
plate. The form of the lift curves for the three airfoils is seen to change 
substantially, however, over the R range from 4.2 x 10 5 down to 0.42 x 10 5 . 
Particularly at the very lowest Reynolds number, the C, versus a curve is no 
longer linear. The flow apparently separates at all positive angles just down¬ 
stream of the minimum pressure point, near the maximum thickness location. 



The N60R airfoil 


Figure 3.64 Airfoil shapes tested at low Reynolds numbers. 
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Figure 3.70 Drag polar for the 417a airfoil at low Reynolds numbers. 


This explanation is substantiated by Figure 3.66. Here Ci versus a is given for 
the N60 airfoil. As a is first increased up to a value well beyond the stall and 
then decreased, a large hysteresis is seen to exist in the curves for the higher 
Reynolds numbers. Typically, as a is increased, complete separation on the 
upper surface occurs at around 12°. The angle of attack must then be 
decreased to around 5° before the flow will again reattach. At the lowest 
Reynolds number, the lift curve tends to follow the portion of the curves at 
the higher Reynolds numbers after stall has occurred and a is decreasing. 
Thus, above an a of approximately 0°, it would appear that the flow is entirely 
separated from the upper surface for the lower R values of 21,000 and 42,000. 

Aerodynamic drag is considered in more detail in the following chapter. 
Nevertheless, the drag characteristics for these low-Reynolds number tests 
are presented now in Figures 3.67 to 3.71. 


Figure 3.71 Drag polar for the N60 airfoil at low Reynolds numbers. 


PROBLEMS 

3.1 A wing has a taper ratio of 1/3, an area of 20 m 2 , and a tip chord of 1.5 m 
What is its aspect ratio? 

3.2 A thin, cambered airfoil is approximated by two straight-line segments 
as illustrated. Calculate C, and C mic for this airfoil according to Equa 
tions 3.39 and 3.41. 
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3.3 The airfoil of Problem 3.2 can be thought of as a flat-plate airfoil at an 
angle of attack with a 50% chord flap deflected through a given angle. 
What are these two equivalent angles? For this a and zero flap angle, 
what would C, be? Comparing this C, to the value from Problem 3.2, 
calculate the flap effectiveness factor t and compare it with Figure 3.32. 

3.4 Taking a cue from Problems 3.2 and 3.3, derive the equation for t given 
in Figure 3.32. 

3.5 A 23015 airfoil is equipped with a 25% fully extensible, double-slotted 
flap deflected at an optimum angle. It has a 6 ft chord and is operating at 
100 mph at SSL conditions. Estimate from: (a) the summary 
observations listed at the beginning of the section on flaps, (b) the 
numerous tables and graphs of data, and (c) Figures 3.32, 3.33, and 3.34. 

3.6 Estimate C t for a thin flat-plate airfoil at a 5° angle of attack having a 
33% c plain flap deflected 15°. Divide the chord into three equal seg¬ 
ments and model the airfoil with three suitably placed point vortices. 

3.7 The GA(W)-1 airfoil of Figure 3.10a is equipped with a pure jet flap. The 
jet expands isentropically from a reservoir pressure of 170 kPa absolute 
and a temperature of 290 °K. The airfoil is operating at SSL conditions at 
15 m/s. The chord is 3 m long and the jet thickness equals 2.5 mm. 
Calculate C, for an a of 2° and a jet flap angle of 30°. 

3.8 A finite wing is simulated by an approximate lifting line model consisting 
of a bound vortex and two vortices trailing from each side, one from the 
tip and the other one halfway out along the span. Using the midspan and 
three-quarter-span stations as control points, calculate the section lift 
coefficients at these stations for a flat, untwisted rectangular wing with 
an aspect ratio of 6 at an angle of attack of 10°. 

3.9 Use Schrenk s approximation instead of the approximate lifting line 
model to answer Problem 3.8. 

3.10 The wing of Problem 3.1 has a washout of 4° and plain .3c flaps over the 
inboard 60% of the span. Assuming R = 9 x 10 6 and a smooth airfoil with 
the characteristics given by Figure 3.6, calculate C L ^ for a flap angle of 
40°. Do this by comparing section C, and C linii values along the span. 

3.11 Write a computer program to solve the lifting line model illustrated in 
Figure 3.52. This is not as difficult and laborious as it may sound. Place 
symmetrically disposed trailing vortices of strength y, at a distance of 
y;(l,2,3,..., n) from the centerline line. Choose control points of 0, 
(>h + 3'2>/2, (y 2 + y 3 )/2,..., (y„_, + y„)/2. At each control point, the bound 
circulation equals the sum of the vortices shed outboard of the point. 
Also, it is easy to show that C, = 2 TfcV. But C, is given by Equation 
3.66, where a, = w/V. The downwash w can be expressed as a sum of 
contributions from each trailing vortex. Hence these relationships lead to 
a system of n simultaneous equations for the unknown vortex strengths 
Vi, yi, - • -, y n Once these are found, T and C, can be calculated. 
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Check your program by calculating the C t distribution for an elliptic 
wing. You should find that C ( is nearly constant except near the tips, 
where the accuracy of the numerical model deteriorates. An n of 20 
should suffice for this example. 

3.12 A cambered airfoil has an angle of attack for zero lift of -4°. If this 
airfoil is incorporated into an untwisted wing having an elliptical plan- 
form, what will the wing lift coefficient be for an angle of attack of 8°? 
The wing aspect ratio is equal to 5.0. 
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FOUR 

DRAG 


As a child, it was fun to stick your hand out of the car window and feel 
the force of the moving, invisible air. To the aeronautical engineer, however, 
there is nothing very funny about aerodynamic drag. A continuing struggle for 
the practicing aerodynamicist is that of minimizing drag, whether it is for an 
airplane, missile, or ground-based vehicle such as an automobile or train. It 
takes power to move a vehicle through the air. This power is required to 
overcome the aerodynamic force on the vehicle opposite to its velocity 
vector. Any reduction of this force, known as the drag, represents either a 
direct saving in fuel or an increase in performance. 

The estimation of the drag of a complete airplane is a difficult and 
challenging task, even for the simplest configurations. A list of the definitions 
of various types of drag partly reveals why this is so. 

Induced Drag The drag that results from the generation of a trailing vortex 
system downstream of a lifting surface of finite aspect ratio. 

Parasite Drag The total drag of an airplane minus the induced drag. Thus, it is 
the drag not directly associated with the production of lift. The parasite 
drag is composed of many drag components, the definitions of which 
follow. 

Skin Friction Drag The drag on a body resulting from viscous shearing 
stresses over its wetted surface (see Equation 2.15). 

Form Drag (Sometimes Called Pressure Drag) The drag on a body resulting 
from the integrated effect of the static pressure acting normal to its 
surface resolved in the drag direction. 

Interference Drag The increment in drag resulting from bringing two bodies in 
proximity to each other. For example, the total drag of a wing-fuselage 
combination will usually be greater than the sum of the wing drag and 
fuselage drag independent of each other. 

Trim Drag The increment in drag resulting from the aerodynamic forces 
required to trim the airplane about its center of gravity. Usually this takes 
the form of added induced and form drag on the horizontal tail. 


Profile Drag Usually taken to mean the total of the skin friction drag and 
form drag for a two-dimensional airfoil section. 

Cooling Drag The drag resulting from the momentum lost by the air that 
passes through the power plant installation for purposes of cooling the 
engine, oil, and accessories. 

Base Drag The specific contribution to the pressure drag attributed to the 
blunt after-end of a body. 

Wave Drag Limited to supersonic flow, this drag is a pressure drag resulting 
from noncanceling static pressure components to either side of a shock 
wave acting on the surface of the body from which the wave is emana¬ 
ting. 

With the exception of wave drag, the material to follow will consider 
these various types of drag in detail and will present methods of reasonably 
estimating their magnitudes. Wave drag will be discussed in Chapter 6. 


SKIN FRICTION DRAG 

Figure 4.1 depicts a thin, flat plate aligned with the free-stream velocity. 
Frequently the drag of a very streamlined shape such as this is expressed in 
terms of a skin friction drag coefficient, Q, defined by, 



where S H is the wetted surface area that is exposed to the flow. This 
coefficient is presented in Figure 4.1 as a function of Reynolds number for the 
two cases where the flow in the boundary layer is entirely laminar or entirely 
turbulent over the plate. Here the Reynolds number is based on the total 
length of the plate in the direction of the velocity. In a usual application, the 
boundary layer is normally laminar near the leading edge of the plate 
undergoing transition to a turbulent layer at some distance back along the 
surface, as described in Chapter Two. The situation is pictured in Figure 4.1, 
where the velocity profile through the layer is shown. In order to illustrate it, 
the thickness of the layer is shown much greater than it actually is. 

As shown in this figure, a laminar boundary layer begins to develop at the 
leading edge and grows in thickness downstream. At some distance from the 
leading edge, the laminar boundary becomes unstable and is unable to 
suppress disturbances imposed on it by surface roughness or fluctuations in 
the free stream. In a short distance the boundary layer undergoes transition to 
a turbulent boundary layer. Here the layer suddenly increases in thickness 
and is characterized by a mean velocity profile on which a random fluctuating 
velocity component is superimposed. The distance, x, from the leading edge 
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Figure 4.1 Drag of a thin flat plate. 

of the plate to the transition point can be calculated from the transition Reynolds 
number, R x . R x is typically, for a flat plate, of the order of 3* 10 s , R x being 
defined by 

(4.2) 

F 

For very smooth plates in a flow having a low level of ambient turbulence, R x 
can exceed 1 x 10*. 

Since the velocity profile through the boundary layer approaches the 
velocity outside the layer asymptotically, the thickness of the layer is vague. 
To be more definitive, a displacement thickness, 5*, is frequently used to 
measure the thickness of the layer. 5* is illustrated in Figure 4.2 and is 
defined mathematically by 

-fM* (4.3) 

where y is the normal distance from the plate at any location such that, 
without any boundary layer, the total flow past that location would equal the 
flow for the original plate with a boundary layer. To clarify this further, let 8 
be the boundary layer thickness where, for all intents and purposes, u = V. 



Figure 4.2 Displacement thickness. 


Then 


V(S-S*) = J udy (4.4) 

Allowing 8 to become infinite leads to Equation 4.3. 

If we arbitrarily define 8 as the value of y at which u = 0.99 V then, for a 
laminar layer, 


5-5.2 = 5.2x( 


(4.5) 


8* 

For a turbulent layer. 


5 

3 


(4.6) 


(4.7) 


(4.8) 


Observe that relatively speaking, the turbulent boundary layer is more 
uniform, with 6* being only one-eighth of 8 as compared to one-third for the 
laminar layer. 

In order to clarify the use of Figure 4.1 and Equations 4.5 to 4.8, let us 
consider the horizontal tail of the Cherokee pictured in Figure 3.62 at a 
velocity of 60.4 m/s (135 mph) at a 1524 m (5000 ft) standard altitude. We will 
assume that the tail can be approximately treated as a flat plate at zero angle 
of attack. 

From Figure 3.62, the length of the plate is 30 in or 0.762 m. The total 
wetted area, taking both sides and neglecting the fuselage, is 4.65 m 2 (50 ft 2 ). 
From Figure 2.3, at an altitude of 1.52 km, p = 0.1054 kg/m 3 and v — 
1.639 x 10 5 m 2 /s. We will assume that the transition Reynolds number is equal 
to 3 x 10 s . 
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_ 60.4 (0.762) 

1.639 x 10 _ * 

= 2.81 x 10 6 

If the flow over the tail were entirely turbulent then, from Figure 4.1, 

Cf = 0.455 (logio J?i) -2S8 (4.9) 

= 0.00371 

The dynamic pressure q for this case is 
q = pV 2 l2 

0.1054 (60.4) 2 
2 

= 1923 N/m 2 

Hence the total skin friction drag would be 
D = qS„C, 

= 1923 (4.65)(0.00371) 

= 33.17 N 

However, the leading portion of the plate is laminar. The wetted area of 
this portion is equal to 0.497 m 2 . For laminar flow over this portion, 

Q = 1.328 R~ in (4.10) 

= 1.328 (3 x 10 5 )-*' 2 
= 0.00242 

Hence the drag of this portion of the plate is equal to 
D — qCfS w 

= 1923(0.00242)(0.497) 

= 2.31 N 


r 
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If the flow were turbulent over the leading portion of the plate, its Q 
would be 

Q = 0.455 (log.oRr 258 
= 0.455 (logio 3 x 10T 2 - 58 
= 0.00566 

Thus its drag for a turbulent boundary layer would be 
D = qC,S„ 

= (1923)(0.00566)(0.497) 

= 5.35 N 

The above is 5.35-2.31, or 3.04 N higher than the actual drag for laminar 
flow. Hence this difference must be subtracted from the total drag of 33.17 N 
previously calculated assuming the boundary layer to be turbulent over the 
entire plate. Hence the final drag of the total horizontal tail is estimated to be 
D = 33.17 -3.04 
= 30.13 N 
= 6.77 lb. 

The thickness, S, of the laminar boundary layer at the beginning of transition 
can be calculated from Equation 4.5. 

S =5.2 (0.0814)(3 x 10 5 )" 1,2 
= 7.728 x IO" 4 m 
= 0.0304 in. 

The thickness of the turbulent layer right after transition is found from 
Equation 4.7 assuming the layer to have started at the leading edge. 

8 = 0.37(0.0814)(3 x 10 5 ) _1/5 
= 2.418 x 10“ 3 m 
= 0.0952 in. 

At the trailing edge, the thickness of the turbulent layer will be 
5 = 0.37(0.762)(2.81 x 10 6 r ,/5 
= 0.0145 m 
= 0.5696 in. 

The displacement thickness at the trailing edge is thus only 0.0018 m 
(0.071 in.). 

Before leaving the topic of skin friction drag, the importance of surface 
roughness should be discussed. Surface roughness can have either a beneficial 
or adverse effect on drag. If it causes premature transition, it can result in a 
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reduced form drag by delaying separation. This is explained more fully in the 
next section. Adversely, surface roughness increases the skin friction 
coefficient. First, by causing premature transition, the resulting turbulent Q is 
higher than C f for laminar flow, in accordance with Figure 4.1. Second, for a 
given type of flow laminar or turbulent, Q increases as the surface is 
roughened. 

It is difficult to quantify the increment in Q as a function of roughness, 
since roughness comes in many forms. For some information on this, refer to 
the outstanding collection of drag data noted previously (e.g.. Ref. 4.4). 
Generally, if a roughness lies well within the boundary layer thickness, say of 
the order of the displacement thickness, then its effect on Q will be minimal. 
Thus, for the preceding example of the horizontal tail for the Cherokee, the 
use of flush riveting near the trailing edge is probably not justified. 

An approximate estimate of the effect of roughness, at least on stream¬ 
lined bodies, can be obtained by examining the airfoil data of Reference 3.1. 
Results are presented for airfoils having both smooth and rough surfaces. The 
NACA “standard” roughness for 0.61-m (2-ft) chords consisted of 0.028-cm 
(0.011-in.) carborundum grains applied to the model surface starting at the 
leading edge and extending 8% of the chord back on both the upper and lower 
surfaces. The grains were spread thinly to cover 5 to 10% of the area. 

An examination of the drag data with and without the standard roughness 
discloses a 50 to 60% increase in airfoil drag resulting from the roughness. It 
is difficult to say how applicable these results are to production aircraft. 
Probably the NACA standard roughness is too severe for high-speed aircraft 
employing extensive flush riveting with particular attention to the surface 
finish. In the case of a production light aircraft for general aviation usage, the 
standard roughness could be quite appropriate. 


FORM DRAG 

In addition to skin friction drag, a body generally experiences some form 
drag. Unlike the skin friction drag that results from viscous shearing forces 
tangential to a body's surface, form drag results from the distribution of 
pressure normal to the body’s surface. The extreme case of a flat plate normal 
to the flow is pictured in Figure 4.3. Here the drag is totally the result of an 
unbalance in the normal pressure distribution. There is no skin friction drag 
present in this case. 

Generally, form drag is difficult to predict. For that matter, so is skin 
friction drag except for the simplest cases. Thus, in general cases, such as that 
pictured in Figure 4.4, where the total drag results from both normal and 
tangential stresses (or pressures) one must usually resort to experimental data 
to estimate the drag. 



Figure 4.3 Flat plate normal to flow. 



Figure 4.4 Body having both skin friction and form drag. 

As with skin friction drag, form drag is generally dependent on Reynolds 
number. To see why, consider the flow around the circular cylinder pictured 
in Figure 4.5. In Figure 4.5a flow is pictured at a low Reynolds number. Here, 
beginning at the stagnation point, a laminar boundary layer develops. On the 
surface of the cylinder, the static pressure (normal) is highest at the stag¬ 
nation point and decreases to a minimum at the top and bottom. Moving 
around toward the rear, beyond these points, the static pressure increases, 
tending toward the stagnation pressure at the very rear. In the absence of 
viscosity the normal pressure distribution would be symmetrical (Equation 
2.78) and there would be no drag. This is a clear example of D’Alembert’s 
paradox, which states that a body in a inviscid fluid will experience no drag. 
As the slower moving fluid in the laminar boundary layer moves beyond the 
minimum pressure point on the cylinder, its momentum is insufficient to move 
against the positive pressure gradient, referred to as an adverse gradient, and 
thus the flow separates just past the top and bottom locations on the cylinder. 
In the separated region over most of the rear portion of the cylinder the static 
pressure is constant and equal to the low pressure at the top and bottom. Thus 
the high pressure acting over the front and the low pressure over the rear 
result in a high form drag. 
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<b) 


Figure 4.5 Flow over a circular cylinder, (a) Low Reynolds number. Separation 
occurs before transition. Large wake, (b) High Reynolds number. Transition 
occurs before separation. Small wake. 


The high-Reynolds number case is shown in Figure 4.5b. Here the 
laminar boundary layer undergoes transition to a turbulent boundary layer 
before separating. The subsequent turbulent mixing increases the momentum 
and energy of the boundary layer so that it remains attached toward the rear 
of the cylinder, well beyond the separation point of the laminar layer. Thus, in 
this case, the separation region is much smaller and the static pressure is 
higher on the rear of the cylinder than for the laminar case. Therefore, 
because of reduced form drag, the drag coefficient of a cylinder is lower at 
higher Reynolds numbers. 

Cd as a function of Reynolds number is presented in Figure 4.6 for both 
spheres and two-dimensional circular cylinders. Here, C d is based on the 
projected frontal area. Note the rapid drop in C d above an R value of 
approximately 2 x 10 6 . This is the so-called critical Reynolds number, where 
the transition point is nearly coincident with the separation point. “Sub- 
critical” refers to flow at Reynolds numbers that are less than critical; 
“supercritical” denotes R values that are higher than critical. A body shape 
having a well-defined separation point will not exhibit a critical Reynolds 
number; neither will streamlined shapes. 

Although not concerned with drag per se, Figure 4.6a also includes the 
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(b) 

Figure 4.6 Drag coefficients of cylinders and spheres versus Reynolds number, 
(a) Two-dimensional circular cylinders, (b) Spheres. 


quantity fD/V , known as the Strouhal number, S. S characterizes an interes¬ 
ting behavior of bluff bodies with rounded trailing edges. As such a body first 
begins to move through a fluid, the vorticity in the boundary layer is shed 
symmetrically from the upper and lower surfaces to form two vortices of 
opposite rotation. However, the symmetrical placement of the vortex pair is 
unstable, so that succeeding vortices are then shed alternately from the upper 
and lower surfaces. The resulting flow pattern of periodically spaced vortices 
downstream of the body is known as a Karman vortex street. 

In the definition of Strouhal number, / is the frequency at which the 
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vortices are shed. As a vortex is shed from one surface of the cylinder, it 
produces a momentary circulation around the cylinder opposite in direction to 
the vortex. From the Kutta-Joukowski law, a force on the cylinder normal to 
V results. As the next vortex is shed, the force reverses its direction, resulting 
in an alternating force on the cylinder. This particular phenomenon is the 
cause for the “singing” of telephone wires in the wind. 

As an example of the use of Figure 4.6 a, consider a wire 2 cm in diameter 
in a wind blowing at a speed of 8 m/s. Assuming standard sea level conditions. 


R = 


Vd 


8 ( 0 . 02 ) 
1.456 x 10'* 
= 1.099 x 10 4 


From Figure 4.6a, for this Reynolds number, 


C d = 1.1 


S = 0.18 


Thus the drag per unit length on the wire will be 
D = qC d d 

= ^ (1.226)(8) z ( l. 1 )(0.02) 
- 0.863 N/m 


The frequency of the alternating lift force on the wire will be. 



8(0.18) 
0.02 
= 72 Hz 


Let us now consider the extreme case of form drag illustrated in Figure 
4.3, where the point of flow separation is well defined and not dependent on 
Reynolds number. It is not too surprising to find that drag coefficients for such 
shapes are nearly constant over a wide range of Reynolds number values. A 
number of such shapes are pictured in Figure 4.7a. 

This figure presents values for both two-dimensional and three-dimen¬ 
sional shapes. Three-dimensional shapes are all bodies of revolution. Observe 
that for the same profile shape, 


CK2-D) 

Cd(3-£>) 


= 1.8 



Figure 4.7a Examples of shapes having C a values nearly independent of Rey¬ 
nolds number. 



Figure 4.7b Transition from three-dimensional to two-dimensional drag for 
cylinders at supercritical Reynolds numbers. 
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If the ratio of the span to the height (or diameter) of a flat plate (or 
cylinder) normal to the flow is approximately 5 or less, C d is nearly constant 
and equal to the 3-D value. For aspect ratios greater than 5, Cj varies 
approximately in the manner given by the normalized curve of Figure 4.7 b. 
This curve is based on data from several sources, including Reference 4.4. 

A qualitative evaluation of the drag coefficient for a given shape can be 
made using some “educated intuition.” Referring to Figure 4.8, the drag 
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coefficient of a bluff shape depends on the width of the wake behind the body 
(before viscosity dissipates it). Beginning with the top figure and working 
downward, one would expect, for the same projected frontal area, that the 
widths of the wakes would diminish progressively. Intuitively, such a pro¬ 
gression is visualized by picturing the flow as separating tangent to the 
surface and then being turned gradually in the direction of the main flow. 

With regard to drag, the trailing edge shape of a body is usually more 
important than the leading edge shape. For example, the drag of the top shape 
in Figure 4.8 can be reduced significantly by providing a body behind it to 
which the flow can reattach. This is illustrated in Figure 4.9. As opposed to 
Figure 4.8, in this case, the low pressure in the separated region between the 
front and the afterbody reacts on both parts, contributing little or nothing to 
the drag. 

In order to provide an additional basis for estimating the drag of two- 
dimensional sections, the data in Figure 4.10 are provided (Ref. 4.5). This 
figure shows that for a shape with sharp corners, a rounding of the corners 
will reduce the drag coefficient as well as the critical Reynolds number. 


DRAG OF STREAMLINED SHAPES 

The drag of shapes such as airfoils, fuselages, nacelles, torpedoes, 
submarines, and airships is composed of both form drag and skin friction 
drag. As the fineness ratio (length/maximum thickness) of a streamlined shape 
increases, more and more of its drag is attributable to skin friction. Con¬ 
versely, at low fineness ratios, the drag is principally form drag. 

Data on the drag of two-dimensional and three-dimensional streamlined 
shapes are presented in Figures 4.11, 4.12, and 4.13. Based on the projected 
frontal area, Q is given as a function of fineness ratio in Figure 4.11 at a 
high Reynolds number. In this figure, a fineness ratio of 1.0 corresponds to a 
circular cylinder and sphere for two-dimensional and three-dimensional 
shapes, respectively. Notice that the minimum drag occurs at a fineness ratio 
of approximately 2 for a three-dimensional shape, and at a value of ap¬ 
proximately 3 for a two-dimensional shape. However, in view of the sharp 
rise in both curves at the lower fineness ratios, it might be well in either case 
to use fineness ratios higher than these, say around 4, if one wishes to fair a 
blunt shape of a given frontal area. 

The crossover of the two curves in Figure 4.11 is to be expected. At low 
values of fineness ratio, Cj for the two-dimensional shapes is higher than that 
for the three-dimensional bodies, based on the data of Figures 4.6 and 4.7. At 
the other extreme, as the fineness ratio becomes large, the skin friction drag 
predominates. If C/ is assumed to be the same for either the two-dimensional 
or three-dimensional shapes, the ratio of the Q values, based on the pro- 



Flgure 4.11 Drag coefficients for streamlined shapes as a function of fineness 
ratio. Gf based on frontal area; ft-10 7 based on length. 
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Figure 4.12 Minimum C a for four and five-digit airfoils of varying camber ratios 
as a function of thickness ratio. C<* based on chord; smooth surface. 
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Figure 4.13 Drag of fuselages and similar shapes. 


jected frontal areas becomes 

C d (3-D) 5 k (3-D) 4 

Q(2-D) S*(2-D) t tD 

where D is the maximum three-dimensional body diameter or the maximum 
thickness of the two-dimensional shape. For an elliptical two-dimensional 
shape compared to an ellipsoid, this becomes 


Cg (3-D) -it 
Cd (2-D) 2 


(4.11) 


This is close to the ratio from Figure 4.11 for a finess ratio of 8 and only 
slightly lower than the corresponding ratio given earlier for the form drag. 

Minimum profile drag coefficients for NACA four- and five-digit airfoils 
are presented in Figure 4.12 as a function of thickness ratio at a Reynolds 
number of 6 x 10 6 . Here, as is usual for airfoils, C d is based on the chord 
length. The several data points at each thickness ratio result from airfoils of 
different camber ratios. Note that does not vary significantly with 
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camber. C dmm appears to vary almost linearly with tic and extrapolates to a 
value of 0.004 for a tic of zero. This corresponds to a C { value of 0.002. 
According to Figure 4.1, this would require laminar flow over these sections 
more extensive than one would expect. Probably, transition is delayed until 
approximately the 25% chord point, the location of maximum thickness. One 
would then expect a O mj[i value of about 0.005. 

Figure 4.13 presents three-dimensional drag data directly comparable to 
Figure 4.11, but with more detail. Data representing practical fuselage and 
nacelle construction are included in Figure 4.13 together with C d results from 
torpedo-shaped bodies. Assuming a reasonable relationship between the 
frontal and wetted areas of such bodies, expected Q values for various 
values of Q are also included on the figure. For a given Q value, the 
experimental results should approach one of these lines as the fineness ratio 
gets large. 

For fully turbulent flow at an R of 25 x 10 6 , Q for a flat plate would be 
0.0026, whereas the data appears to be approaching a Q of 0.0032 to 0.0034. 
The higher skin friction drag on the bodies is probably the result of surface 
roughness. 

It is interesting to examine the data of Figure 4.13 in terms of minimum 
drag for a given body volume. This is particularly important for airship and 
underwater applications. It is also of interest to the design of tip tanks, where 
minimum drag for a given volume of fuel is desirable. Denoting the volume by 
V m , we will define another drag coefficient. 


Q v is related to C d in Figure 4.13 by 


Obviously, the ratio of the frontal area, A, to the 2/3 power of the volume 
depends on the particular body shape. We will assume the body to be 
composed approximately of a hemispherical nose, a cylindrical midbody 
extending to the middle of the body, and a tail cone. For this particular shape, 


Using this relationship and Figure 4.13, the graphs presented in Figure 4.14 
were obtained. From this figure it can be seen that to enclose a given volume 
with a minimum drag body, its fineness ratio should be higher than the 
optimum values from Figure 4.13. Indeed, for fuselages, the drag for a given 
volume is nearly constant for lid values from 4 to 10. 


